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of Some lron-Thionine Photogalvanic Cells 


By A. E. Potter, Jr., and L. H. Thaller 


Lewis Research Center, National Aeronautics and Space Administration, Cleveland, Ohio 


Rabinowitch has suggested that the photogal- 
vanic effect might be used to convert sunlight 
into electricity. In order to explore this sugges- 
tion, some photogalvanic cells using the iron- 
thionine system as the photosensitive fluid were 
built and tested. The observed maximum power 
conversion efficiency was 3 X 10“ per cent, 
based on light absorbed. The principal reason 
for the low efficiency was shown to be polariza- 
tion of the polished platinum electrodes. Coating 
the electrodes with platinum black reduced 
polarization sufficiently, so that it appeared 
possible to achieve an efficiency of 6 
per cent, although this value was not actually 
observed. In principle, it appeared possible to 
make further increases in efliciency by increasing 
electrode area and decreasing the electrolyte 
resistance. 

Aside from electrode polarization, the major 
problem encountered was the rapid loss of the 
photochemical activity of the dye. No further 
progress can be made until this problem is 


solved. 


INTRODUCTION 

When certain dyes are dissolved in a solution con- 
taining a suitable reducing agent, the solution becomes 
photochemically active, in the sense that the dissolved 
dye can be reduced to a colorless form by exposure to 
light." ~ When illumination is terminated, the reaction 
reverses, and the dye reappears in its original form. A 
change in the potential of the solution accompanies the 
photoreduction of the dye. This latter phenomenon has 
been called the photogalvanic effect by Rabinowitch”’ ° 
who studied the iron-thionine system. This system can 
be represented by the equilibrium 

Fe'' + thionine @ Fe'"' + reduced thionine 


In the dark, this equilibrium lies far to the left as a 
result of the relatively large difference between the 
oxidation-reduction potentials of the iron and dye 
systems. Under strong illumination, the equilibrium is 
shifted to the right, and the solution bleaches. The 
basic reason for the photosensitivity of the equilibrium 


appears to be the formation of a long-lived electroni- 
cally excited state by thionine upon absorption of 
light.’ This state is very reactive and rapidly oxidizes 
ferrous ion. In the bleached solution, the major reduced 
dye species appears to be semithionine.'’ With intense 
illumination, the potential of the bleached solution can 
be as much as 0.25 v more negative than that of the 
unbleached solution. The relation between the amount 
of bleaching and the potential of the solution is very 
complex, and poses an unusual problem which has not 
vet been satisfactorily solved. This research concerns 
the possible use of the potential difference between il- 
luminated and dark iron-thionine solutions for con- 
version of sunlight into electricity. This use has been 
suggested by Rabinowitch, who states that an effi- 
ciency approaching | per cent based on light absorbed 
Was attained in preliminary experiments.'® re- 
ported efficiency was thought to be sufficiently promis- 
ing to justify an attempt to uncover the factors limiting 
efficiency in this cell, with the hope that the efficiency 
could be improved. For this reason, some iron-thionine 
photogalvanic cells were constructed, and their effi- 
ciency measured. 
EXPERIMENTAL DETAILS 

The Photogalvanic Cells 

Basically, a photogalvanic cell consists of a pair of 
inert electrodes immersed in a solution of dye and re- 
ducing agent. One electrode and the solution surround- 
ing it is illuminated. The other electrode is not il- 
luminated. A potential difference between the two 
electrodes is produced as a result of concentration 
differences at these electrodes, and may be used to 
perform electrical work. The potential of the illumi- 
nated solution is such that electrons flow through an 
external circuit from the illuminated electrode to the 
dark electrode. 

Two phenomena connected with the transfer of 
electrons from the light to the dark solution must be 
considered in the design and operation of photo- 
galvanic cells. They are: (1) concentration polarization 
of the bulk of the solution, and (2) competition be- 
tween the electrode reaction and the solution reaction. 

Let us first consider concentration polarization of the 
bulk of the solution. In the opinion of Rabinowitch,’ 
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the principal species which delivers electrons to the 
illuminated electrode is some form of the reduced dye, 
probably semithionine, since it is the predominant dye 
species in the illuminated solution.” The ionic species 
which accepts electrons from the electrode in the dark 
solution is ferric ion. Thus, when electrons are trans- 
ferred from illuminated to dark solution through an 
external circuit, the following electrode reactions occur: 


Illuminated electrode 
semithionine — thionine + electron 
Darkened electrode 
+ electron Fe" 


In the illuminated solution, the thionine formed at the 
electrode will react with Fe"' in the presence of light to 
reform semithionine and to produce Fe’. The net 
result of the electron transfer is thus an accumulation 

Fe’ in the illuminated solution and of Fe™ in 
the dark solution. The resultant potential difference 
counters that produced by photoreduction of the dye. 
Thus, if sufficient power is withdrawn from the cell, 
it becomes polarized because of concentration dif- 
ferences between the solutions. Two approaches to a 
solution of this problem were made. The simplest 
way to describe them is to discuss the particular 
cells which were built. These are shown in Figs. 1(a) 
and (b). Cell A, shown in Fig. 1(a), consists of two 
polished plates of platinum cemented to a glass back- 
ing plate and covered with a thin (l-mm) layer of 
photogalvanic solution. One plate or electrode is 


illuminated, the other is not. The intended method of 
operation was as follows: power would be withdrawn 
from the cell until volume polarization reduced the 
output appreciably. Then, the illuminated electrode 
would be darkened, and the previously darkened 
electrode would be illuminated. Current would then 
flow in the opposite direction than previously. When 
volume polarization again reduced power output, 
the operation would be repeated. By see-sawing the 
cell in this manner, the effect of volume polarization on 
power production can be minimized. The second 
approach to the volume polarization problem is seen 
in cell B, shown in Fig. 1(b). This cell consists of a 
sintered glass disk whose front and back are covered 
with a thin porous layer of gold. The gold layers form 
the two electrodes of the cell. The disk is immersed in a 
photogalvanic solution. One face, or electrode, is 
illuminated, the other is not. When power is withdrawn 
from the cell and volume polarization sets in, a con- 
centration gradient is produced across the sintered 
glass separating the two electrodes and the illuminated 
and dark parts of the solution. As a result of this 
gradient, interdiffusion between the two solutions will 
occur which will tend to counter the volume polar- 
ization. 

The second and perhaps most important factor 
affecting cell design is the competition in the illumi- 
nated solution between the electrode reaction and the 
back reaction in solution. The competing reactions 
are as follows: 


ELECTRODE REACTION 


A 


— dye + electron 
reduced dye 
+ oxidizing — dye + reducing 


SOL U TION REACT ION” 
Unfortunately there does not seem to be any really 
effective means of gaining a favorable balance in this 
competition. This difficulty is one of the major factors 
which limit cell efficiency. All that can be done is to 
make the ratio of electrode area to solution volume 
as large as possible. For this reason, the depth of 
solution in the cell A was made as small as possible 
(if the depth is too small, almost no light is absorbed 
by the cell), and the electrode in cell B was of a rough, 

porous nature, thus having a large surface area 


The Photogalvanic Solution 

It seemed reasonable to use a solution giving the 
highest possible photogalvanic potential with the 
light source used. The chemical variables which con- 
trol the photogalvanic effect are the concentrations of 
Fe", Fe'", and thionine. The concentrations of 
these species were varied in a systematic manner to 
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find the best combination. The procedure was as 
follows. A glass half cell was filled with a known 
volume of an air-free solution of a particular H*, Fe", 
and Fe’ concentration. The anion of all the solutions 
studied was Cl’. The solution was then illuminated 
and its potential measured with a polished platinum 
electrode while a concentrated thionine solution was 
added from a burette. The solution was constantly 
stirred by nitrogen bubbled through it. When addition 
of thionine caused no further decrease in the potential 
of the solution, the titration was stopped. This poten- 
tial L, of the solution and the amount of thionine 
added to achieve £, were recorded. The illumination 
was then turned off, and after the solution returned to 
equilibrium the dark potential 2» was measured. 
The H*, Fe", and Fe’ mixtures studied by this 
technique were all the combinations of the following 
concentrations (in moles/liter ) : 


[Fe] = 0.5, 1, and 2X 10° 
[Fe] = 0.2, 1, and 5 xX 10° 
(H*]=1, 3, and 6X 10” 
The solution which yielded the largest photogalvanic 
potential of all (2, — Ep = 0.185 v) had the follow- 
ing composition: 
[Fe"] = 1.0 x 10° 
[Fe] = 5.0 x 10° 


[H*] = 3.0 x 10° 


(Thionine] = 1.5 X 

The composition of this solution is almost identical 
with that found by Rabinowitch to give a large poten- 
tial change. This solution and one having double the 
thionine concentration, but otherwise identical, were 
used in all the tests reported. It was found that the 
solution deteriorated rather rapidly, so that it was 
necessary to prepare it under nitrogen immediately 
before use from air-free stock solutions stored under 
nitrogen. The stock solutions also deteriorated rapidly, 
and were never more than 48 hours old. The photo- 
galvanic potential of the solution (2, — Hp) was 
always measured just prior to use in the photogalvanic 
cells. 


The Light Source and Radiant Energy Measure- 

ments 

The light source used in all the experiments was a 
750-watt photoflood lamp. The light was _ filtered 
through 3 em of running water and 8 cm of a saturated 
solution of ferrous sulfate. These filters removed the 
portions of the spectrum which would have been ab- 
sorbed by water and the ferrous ions in the photo- 
galvanic solution. The energy from the light was 
measured with a calibrated Epply thermopile in con- 


junction with a microvoltmeter. The light beam from 
the source was only approximately parallel, so that the 
radiant intensity varied with distance from the lamp. 
For this reason, when measurements of radiant energy 
incident on the photogalvanic cells were made, care 
was taken to make the measurement at the same dis- 
tance from the lamp as the cell was placed. 

The amount of radiant energy absorbed by the 
photogalvanic cell A during operation was estimated 
as follows: the cell was turned from its usual position 
normal to the incident light beam, so that it formed a 
slight angle from the normal to the incident beam. 
The beam reflected from the polished electrode could 
then be observed, and its intensity measured. The 
intensity of the reflected beam was measured with cell 
empty, and with the cell filled with the photogalvanic 
solution. Calling the energy reflected from the empty 
cell J, and the energy reflected from the full cell J, , 
it is evident that 7, — J; is a measure of the energy 
absorbed by the solution, and that (J, — J;)/J, 
is approximately the fraction of incident light absorbed 
by the solution. Therefore, the energy absorbed by the 
cell is approximately the energy of incident light times 
the fraction (J, — J;)/7, . A similar procedure could not 
be applied to cell B, because the rough surface of the 
gold electrode reflected light diffusely rather than 
specularly. For cell B, a crude estimate of the light 
absorbed was made as follows: it was assumed that 
light absorption in the cells obeyed Lambert’s law. 
Then an empirical absorption coefficient was calculated 
from the light absorbed in cell A, and used to estimate 
the light absorbed in cell B. The amount of absorbed 
radiant energy calculated by this procedure is certainly 
accurate to within an order of magnitude, and is 
probably better than that. It will be seen that the 
accuracy is sufficient for the purposes of this research. 

For the power measurements reported here, the 
amount of radiant energy incident on the cells was 
about 3 X 10° watts per sq cm. Between 10 and 20 
per cent of this energy was absorbed by the cells. It 
was found that sunlight of about ten times the intensity 
of the filtered lamplight was needed to produce the 
same cell voltage as the lamplight. Thus, the cells 
can effectively use only about one per cent of the energy 
of sunlight. 

POWER AND EFFICIENCY OF THE CELLS 

In order to eliminate spurious voltages arising from 
dirty electrodes, it was necessary to soak the cells in 
concentrated HCl for a day or two before use. The 
voltage of the cells when exposed to light from the 
lamp and filled with the photogalvanic solution con- 
taining 1.5 X 10° moles per liter thionine was only 
about 80 mv, less than half the expected value of 185 
mv. Possibly this was because of leakage of light from 
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taining 3 X 10° moles per liter thionine was used in 
cell A, a voltage of 126 mv was observed. Possibly 
less light leaked to the dark side of the cell when the 
more concentrated dye solution was used, so that a 
higher voltage was obtained than before. The cell 
with this solution vielded a voltage of 170 mv in 
direct sunlight of 0.068 watt per sq cm. It was esti- 
mated that a voltage of 182 mv would have been 
produced in the most intense sunlight available at 
the earth’s surface. 

The power output of the cells was determined by 
measuring the voltage drop across a resistance con- 
nected across the cell terminals. Resistances ranging 
from 10° to 10° ohms were used. Current was withdrawn 
only for very short periods of time to minimize con- 
centration polarization in the bulk of the solution. 

It was found that the cell would never return to its 
initial voltage after power had been withdrawn from 
it. This was the result of irreversible changes in the 
dve. Thus, in one case, the initial cell voltage was 
126 mv. After about three minutes of operation at 
maximum current, the cell voltage dropped to 105 
my, and could not be increased to the original value in 
spite of the fact that concentration polarization was 
eliminated by back-electrolysis of the cell solution. 
Paper chromatographs of the solution from the cell 
showed that a new, nearly colorless compound had 
appeared in the solution. Thus, side reactions occurred 
which converted the dye into a substance having no 
photochemical activity. Because of this discouraging 
result, it was not thought worthwhile to operate the 
cells for the half-hour or so required to produce suffi- 
clent concentration polarization in the bulk of the 
solution to test the means for combating it. 

The efficiency of the cells was calculated by dividing 
the power output by the radiant energy absorbed by 
the cell. Since so little power could be withdrawn from 
the cells, it seemed reasonable to assume that the 
withdrawal of power did not appreciably affect the 
absorption of radiant energy This considerably simpli- 
fied the measurements, since the radiant energy ab- 
sorbed by the cells could be assumed constant and 
equal to that absorbed when no power was withdrawn. 
This assumption also has the result that maximum 
efficiency is observed when maximum power is de- 
veloped by the cell. This is in contrast to the common 
galvanic cell where maximum efficiency is achieved 
as the power approaches zero. The cell efficiencies are 
shown in Fig. 2. Cell A was tested with photogalvanic 
solutions containing 1.5 and 3.0 K 10° moles per 
liter of thionine. Cell B was tested with only the solu- 
tion containing 1.5 & 10 ° moles per liter. The highest 
efficiency is observed for cell A with the high concen- 


one side of the cell to the other. When a solution con- 
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Fic. 2—Power efficiency as a function of load resistance for two 
photogalvanie cells. 


tration of thionine. Even for this case, the maximum 
efficiency is very low, hardly reaching 3 X 10 per 
cent. 

It can be seen in Fig. 2 that the load resistance at 
which the cells reach maximum power and efficiency 
is in the range of 10° to 10° ohms. It can be shown that 
the power delivered to a load resistance by a simple 
battery is a maximum when the load resistance equals 
the internal resistance of the battery. Consequently, 
the apparent internal resistance of the cells is 10° to 
10° ohms. However, the resistance of the electrolyte 
in the cells is caleulated (from cell dimensions and the 
conductance of the solution) to be of the order of 10° 
ohms. Thus, some factor other than just the resistance 
of the electrolyte limits the power output of the cells. 
This factor will be shown later to be chemical polariza- 
tion of the electrodes. 

The experimental results pose two questions: (1) 
Why is the cell efficiency so small? (2) What measures 
are required to increase the cell efficiency? 

In order to discuss these questions, it is helpful to 
have some idea of the maximum efficiency that could 
in theory be attained by a photogalvanic cell. It can 
be shown that the maximum power that can be de- 
veloped in a circuit consisting of a battery, the internal 
resistance of the battery, and an external load resist- 
ance is £.°7,, where E,’ is the zero-current voltage, 
and /, is the short-circuit current yielded by the bat- 
tery. At most, one-half of this power can be developed 
in the external load resistance. The other half of the 
power is lost in the internal resistance of the battery. 
The maximum power output is thus }/£,'/, . The power 
efficiency 7 is defined as the ratio of power input to 
power output. Since we have assumed a constant 
power input, the maximum power efficiency for the 
photogalvanie cell is 

_ LE, 
where ¢ is the energy per absorbed photon, and N 


is the number of photons absorbed per second. As 
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written, the equation is strictly correct only for mono- 
chromatic light. This is a convenient approximation 
for the iron-thionine system, since most of the light 
is absorbed near 5900 A. Equation [1] shows that even 
if the conversion of light into electric current is per- 
fectly efficient, so that E£,°7, = eN, the power effi- 
ciency of the cell is 50 per cent. This is because the 
cell cannot deliver more than half the absorbed power 
to an external load. 

Equation [1] can also be used to estimate the maxi- 
mum efficiency that is possible for the cells used in 
this research. For this purpose, it is convenient to 
consider the power efficiency as the product of two 
other efficiencies, one /./N, the other E.’/e. The 
ratio J./N can be considered as a photon efficiency; 
the number of electrons produced in the circuit per 
photon absorbed. If one assumes that each molecule 
of reduced dye formed by the photochemical reaction 
yields one electron at the electrode, the maximum 
value of the photon efficiency must be the quantum 
yield of the photoreduction. The data of Hardwick" 
may be used to calculate the quantum yield for the 
solution used in the cells. The value obtained is 0.7, 
corresponding to a maximum photon efficiency of 
70 per cent. The ratio 2,’ may be regarded as a volt- 
age efficiency. Its value may be calculated, since E,” 
was measured, and the energy per photon absorbed 
by thionine is about 2.1 ev. As mentioned previously, 
cell A in maximum sunlight will yield an open-circuit 
voltage EF.” of 0.182 v, which gives a maximum voltage 
efficiency of 8.7 per cent. When these maximum theo- 
retical values of the photon and voltage efficiency are 
used in Equation [1], it is found that the maximum 
power efficiency is 3 per cent. The voltage efficiency 
is the major factor making the theoretical power effi- 
ciency so small. The voltage of the cell is small because 
the back reaction between reduced dye and oxidizing 
agent keeps the concentration of reduced dye low. 
It does not seem reasonable to expect voltages much 
larger than about 0.5 v in sunlight, corresponding to a 
possible power efficiency of about 12 per cent, assuming 
a quantum yield of unity. 

Returning to the original question, why is the ob- 
served efficiency a factor of about 10° less than the 
theoretical maximum? For the answer, compare the 
observed photon efficiency with the theoretical maxi- 
mum of 70 per cent. The largest observed photon 
efficiencies are of the order of 10° per cent. There- 
fore, the poor efficiency of the cells is primarily the 
result of a very small photon efficiency, which is to 
say, a very small current from the cell. 

A possible cause of the smallness of the current from 
the cell is electrode polarization. Electrode polarization 
will cause the cell potential to decrease as the current 


withdrawn from the cell increases, and thereby will 
limit the maximum current that can be withdrawn. 
In order to measure this effect, the potential of a small 
platinum electrode immersed in the photogalvanic 
solution was measured as a function of current density. 
The experimental arrangement was similar to that 
described on page 436 of Ref. 12. The overvoltage, 
the difference between the reversible electrode poten- 
tial and the observed electrode potential, is plotted 
against current density in Fig. 3. Results are given for 
both anodic and cathodic polarization. In the photo- 
galvanic cell, the electrode in the dark solution is the 
anode, and the electrode in the illuminated solution is 
the cathode. Accurate measurements of the cathodic 
polarization in illuminated solutions could not be 
obtained because the passage of current slowly de- 
stroyed the photochemical activity of the dye, so 
that the reversible electrode potential changed as 
current flowed. However inaccurate they may be, the 
results that were obtained for cathodic polarization in 
the illuminated solution are seen in Fig. 3 to be very 
similar to those obtained for cathodic polarization in 
the dark solution. This is a somewhat surprising result, 
since the cathode reaction should be different in the 
dark and illuminated solutions. Both anodic and 
cathodic overvoltages are directly proportional to the 
current density at low current densities. At higher 
current densities, the overvoltages become propor- 
tional to the logarithm of the current density. This 
behavior is typical of the iron group metals and is 
displayed when a reaction at the electrode is the rate- 
determining step in the discharge of an ion at the 
electrode (Ref. 12, p. 462). 

The experimental results shown in Fig. 3 may be 
used to calculate the variation of voltage of the iron- 
thionine photogalvanic cell with current density. 
At low current densities (0.2 wa per sq cm or less), 
the result is 

E. = — 10° I, [2] 


where F is the cell voltage, /, is the current density 
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voltage when no current is withdrawn. This result may 
be used to show both why the current from the cell 
is so small and why the load resistance at maximum 
efficiency does not equal the internal resistance of the 
cell as mentioned previously. Equation [2], which 
applies to the iron-thionine photogalvanic solution 


used here, can be written as 


where J is the current from the cell, A is the electrode 
area, and b is the proportionality factor between the 
current density and extent of electrode polarization 
(equal to 10° ohm per sq cm for the case at hand). 
The power developed by the cell in an external re- 
sistance R, is 
P=I'R. 

The current in the circuit 1s 

E. 

+R 
where FR; is the internal resistance of the cell. Com- 
t|, and [5], we find that 


E. R. | 6] 


b \- 
G +R; + 


This equation may be used to show that the maxi- 


bining Equations [5}, 


mum power is developed by the cell in the external 
resistance when 

The current when maximum power is being trans- 
ferred is therefore 

[8] 

Then it follows from Equation [4] that the maximum 
power developed by the cell in an external resistance is 


l E. 0 


9 
+ (b A)| 9) 


Prax = 

The value of 6 as found from the electrode polariza- 
tion measurements is 10° ohm per sq cm. The exposed 
area of each electrode of cell A was 6.25 sq em. The 
value of b/A is therefore 1.6 X 10° ohms. Since the 
internal resistance of the cell was very much smaller 
than this, of the order of 10° ohms, Equation [7] indi- 
cates that the power, and hence the efficiency for 
constant power input, of cell A should be a maximum 
at a load resistance of 1.6 X 10° ohms. Examination 
of Fig. 2 shows this to be approximately true. Cell 
A with the thionine concentration at 1.5 X 10° 
molar displays a maximum efficiency at 2.5 X 10° 
ohms, and with the thionine concentration at 3.0 X 


10° molar, the efficiency maximizes at about 1 X 10° 
ohms. The reason for the difference between the two 
is not clear. The reason why cell B has maximum 
efficiency at a lower load resistance, 5 X 10° ohms, is 
probably because of the greater electrode area of cell B, 
which makes b/A smaller. (The electrode area of B 
is large because the electrode surface is very rough. 
The actual frontal area of cells A and B is almost the 
same. ) 

Examination of Equation [8] shows that the current 
at maximum cell power is limited by b/A for a fixed 
value of 2," and for negligible internal resistance of the 
cell. Thus, the current and therefore the photon 
efficiency is very low for the cell because the elec- 
trodes are so easily polarized, that is, 6 is very large. 
The major factor which limits efficiency of the iron- 
thionine photogalvanic cell is seen to be the easily 
polarized electrodes. 

It is standard physical chemical laboratory practice 
to reduce or eliminate the polarization of platinum 
electrodes by coating them with platinum black. This 
was done for cell A, and attempts were made to repeat 
the efficiency experiments with the following curious 
results: the cell voltage did not increase quickly to a 
steady value upon illumination as before. Rather, it 
increased very slowly. With the usual photogalvanic 
solution, it increased from zero to 30 or 40 mv over a 
period of an hour. When citric acid was used instead of 
hydrochloric acid to acidify the solution, the cell volt- 
age rose to 245 mv in 45 minutes under illumination 
of the 750-watt lamp. Obviously, efficiency measure- 
ments on such sluggish cells has little value unless 
one operates the cell at each load long enough to 
achieve equilibrium, about an hour in this case. This 
is not practical, since the dye deteriorates. Thus, a 
rapid series of power measurements were made to 
determine the effect of blacking the electrodes on cell 
polarization by measuring the load resistance at 
maximum power output. This resistance is equal to 
the internal resistance of the cell plus the apparent 
internal resistance caused by electrode polarization. 
The result was 5 X 10° ohms. The previous result for 
cell A with polished platinum electrodes was 1 X 10° 
ohms. Thus, coating with platinum black reduces the 
electrode polarization a great deal, by a factor of 200 
in this case. As a result, the efficiency could, in prin- 
ciple, be increased by a factor of 200, assuming the 
low voltages and peculiarly sluggish behavior of the 
cell observed with platinum black electrodes does not 
signal some new undesirable phenomenon. Unfor- 
tunately, a factor of 200 brings the efficiency up to 
only 6 X 10° per cent. This is still not satisfactory, 
and further measures are required to increase the 
photon efficiency by increasing the current. 

Equation [8] shows that the cell current is inversely 
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proportional to (R; + b/A), where FR; is the resistance 
of the electrolyte in the cell, b is the polarization con- 
stant, and A is the electrode area. The polarization 
constant b has been reduced to a minimum by blacking 
the electrodes. Any further increase in current must 
be the result of decreasing the electrolyte resistance 
R,; and increasing the electrode area A. The former 
could be done by adding an inert electrolyte such as 
KCl to the solution. The latter could be done in several 
ways, such as making the electrode porous and pump- 
ing the solution through it, or by making the electrode 
from a number of slightly separated thin plates placed 
parallel to the light. In order to increase the current 
to a value such that the 3 per cent theoretical efficiency 
could be achieved, the electrolyte resistance must be 
decreased by a factor of about 50, and the electrode 
area must be increased by a factor of about 250. In 


principle, such changes are possible, so that it should 


prove possible to increase the current enough to ap- 
proach the theoretical efficiency. 

It should be pointed out that the dye used in this 
work absorbs only a small fraction of the solar spec- 
trum. For a practical cell, other photochemically active 
dyes which absorb in different spectral regions would 
have to be found, so that a mixture of such dyes could 
be used to absorb a large fraction of the solar spectrum. 
This is a smaller problem than that of the stability of 
the dye, since many dyes of different colors are known 
to be photochemically active. 


CONCLUSIONS 


A photogalvanic cell was constructed which would 
yield a voltage of 0.182 v in the brightest sunlight 


available at the surface of the earth. The theoretical 
maximum efficiency of the cell in sunlight was caleu- 
lated to be 3 per cent, based on light absorbed. The 
observed maximum efficiency was 3 X 10° per cent, 
based on light absorbed. The reason for the low ob- 
served efficiency was shown to be the ease of polariza- 
tion of the electrodes. Partial elimination of electrode 
polarization was achieved by coating the electrodes 
with platinum black. However, the reduction of polar- 
ization was insufficient to make the possible efficiency 
larger than 6 X 10° per cent. Further increases in 
efficiency appeared to be possible by increasing the 
electrode area and decreasing the electrolyte resist- 
ance, so that an efficiency of 1 per cent is not incon- 
ceivable. However, a major problem remains, that of 
chemical instability of the dye. No further progress is 
possible until this problem is solved. 
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The Photovoltaic Effect and its Utilization’ 


By Paul Rappaport 


RCA Laboratories, Princeton, N. J. 


The growth of the photovoltaic effect from an 
interesting scientific phenomenon into one of the 
most efficient methods yet devised for converting 
solar energy directly into electrical energy is de- 
scribed. Attention is given to considerations of 
the optimum semiconductor material for solar 
energy conversion and how it changes with 
operating temperature. Experimental results 
are presented showing operating characteristics 
and conversion efficiencies of such materials as 
silicon, gallium arsenide, indium phosphide, 
cadmium sulfide, and others. The high voltage 
(~100-y) photovoltaic effect from evaporated 
CdTe films is described. Included is a discussion 
of the major problems yet to be solved. 


INTRODUCTION 

It is the purpose of this paper to present a current 
review of the photovoltaic effect with special emphasis 
on the problem of solar energy conversion. A testi- 
monial to the importance of this field today is the 
Vanguard Earth Satellite that has been sending solar- 
powered signals back to Earth for a year and a half 
now as a direct result of the photovoltaic cells on the 
Vanguard’s surface. In the future as power needs 
become greater and conventional energy sources 
(fossil fuel and fissionable materials) are used up, the 
direct conversion of solar energy striking the waste 
lands of the world could conceivably become an im- 
portant source of power. 

A major step towards the solution of this difficult 
problem has been achieved by the solar converter or 
photovoltaic cell. This device is sometimes called a 
solar battery; however, strictly speaking, it is incorrect 
to use the term battery, since the primary energy 
source is not self contained. Such cells have already 
been made that can convert the sun’s radiated energy 
directly into electrical energy with an efhciency of up 
to 14 per cent, the best experimental over-all efficiency 
reported to date’ for any artificial conversion method. 


* Much of the work discussed in this paper was supported 
by the United States Army Signal Research and Development 
Laboratories. Also published in the RCA Review, September, 
1959. 


BRIEF HISTORY 


It has been over 100 years since Becquerel in 1839 
first discovered that a photovoltage was developed 
when light was directed onto one of the electrodes in 
an electrolyte solution.” Adams and Day* were the 
first to observe the effect in a solid, selenium, about 
forty years later. A number of other early solid state 
workers including Lange,’ Grondahl,’ and Schottky’ 
did pioneering work on selenium and cuprous-oxide 
photovoltaic cells. This work eventually resulted in the 
photoelectric exposure meter that has become so im- 
portant in photography. However, it was not until 
1954 that workers became interested in energy sources 
utilizing the photovoltaic effect. First there was the 
work of the RCA group’ who demonstrated practical 
efficiencies in converting radioactive radiation into 
electrical energy using a silicon p-n junction photo- 
voltaic cell. Then using a similar method Chapin, 
Fuller, and Pearson’ reported a solar conversion effi- 
ciency of about 6 per cent. A third paper appeared in 
1954 by Reynolds, Leis, Antes, and Marburger’ who 
reported about 6 per cent solar conversion efficiency in 
cadmium sulfide p-n junctions. Since 1954 the silicon 
and CdS work has been pursued with vigor. Other 
materials and techniques were also developed as will 
become evident later in the discussion. 


DESCRIPTION OF PHOTOVOLTAIC 
EFFECT 


The photovoltaic effect can be defined as the genera- 
tion of a potential when radiation ionizes the region in 
or near the built-in potential barrier of a semiconductor. 
It is characterized by a self-generated emf and current 
that can deliver power to a load; the primary power 
coming from the ionizing radiation. A potential barrier 
can be formed in a semiconductor by several means. 
The two techniques of most interest to photovoltaic 
cells are (1) the metal-semiconductor barrier formed 
by the deposition of a suitable transparent metallic 
film by evaporation or sputtering onto a semicon- 
ductor as in the selenium cell, or (2) the p-n junction 
formed by the introduction of a suitable type im- 
purity into the opposite impurity type semiconductor. 

We will now attempt to show how the photovoltaic 
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Fic. 1—Electron-hole production by photons. 


effect comes about when light strikes a semiconductor 
p-n junction. The left side of Fig. 1 shows how carriers 
are generated when light photons enter the semicon- 
ductor crystal. The circles are the atoms in this two- 
dimensional crystal lattice, and the double lines repre- 
sent the two electrons in the covalent bonds as in the 
case for silicon. Photons with energy above that re- 
quired to break these bonds, i.e., 1.1 ev for silicon, will 
create electron-hole pairs on a unit quantum efficiency 
basis as shown on the diagram. The right side of the 
figure shows the same process using the energy band 
model for a semiconductor. Here are seen both an n- 
and p-type semiconductor. In one case, a column V 
element contributes extra electrons, and in the other 
case a column III element does not contribute enough 
electrons and so contributes holes. 

The n-type material has its Fermi level (denoted by 


f in Fig. 1) near the top of the forbidden gap with 


many electrons, n, in the conduction band and few 
holes, p, in the valence band. The opposite is true in 
the p-type material. At a given temperature the 
product np is a constant (np ~ 10” for silicon at room 
temperature). For n-type highly conducting Si, n can 
be 10" per cu cm and p would therefore be 10° per 
cu cm. We therefore develop the concept of the 
majority and minority carriers, the electrons and holes 
respectively in the cited case. When sunlight strikes 
this semiconductor, those photons with energy greater 
than the forbidden gap energy produce both type 
carriers in equal numbers. The net effect is that an 
intense light source can increase the minority carrier 
density by many orders of magnitude while the ma- 
jority carrier density is barely affected. These photo- 
induced carriers are in excess of the thermal equilibrium 
number and they will diffuse randomly about the 
semiconductor and recombine in times of the order of 
tenths of microseconds. 

If the n and p-type semiconductors are brought 
together, a p-n junction with a potential barrier is 
formed since thermodynamics requires that the average 
energy of the carriers (the Fermi level) be the same in 
the two materials. Now, excess carriers that are within 


a diffusion length (the average distance minority 
carriers diffuse before they recombine) of the potential 
barrier will be “trapped” by the barrier and caused to 
flow across it in an attempt to reduce their energy. 
This is shown in Fig. 2 at the bottom where the energy 
level diagram of a p-n junction is illustrated. The 
excess electrons will flow to the right and excess holes 
to the left. These constitute an electrical current, and 
when suitable electrical connections are made such a 
device will convert radiation into electricity. The cur- 
rent produced will be proportional to the number of 
photons absorbed, and the voltage will depend on the 
height of the barrier, which is always less than the 
width of the energy gap in the semiconductor depend- 
ing on how heavily doped the p and n regions are. 

At the top of Fig. 2 is seen a typical commercially 
available silicon solar cell. Boron is diffused into an 
originally n-type wafer at about 1100°C for 20 minutes, 
causing a p-type skin a few microns thick to form on 
all exposed surfaces. All surfaces are etched or lapped 
down to the n region except the top face. Ohmic con- 
tact is made to the top and bottom by first nickel 
plating and then soldering. Leads are connected to the 
load from the ohmic contacts, and light striking the 
top face will deliver a current to the load. 

In order to achieve high efficiency of conversion it is 
desirable to produce an electron-hole pair within a 
diffusion length of the junction for each photon ab- 
sorbed. For short diffusion length material (which is 
the usual result after the high temperature firing re- 
quired to form the p-n junction) it is therefore de- 
sirable to absorb most of the photons in the region 
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where the field exists. The field can be tailored to some 
extent to match the absorption of the light by varying 
the time and temperature of the boron diffusion process; 
however, because of the variation of absorption coeffi- 
cient with \ a compromise must eventually be ac- 
cepted. Fig. 3, kindly supplied by Dr. Prince of Hoff- 
man Semiconductor Division, shows how the genera- 
tion of carriers varies with wavelengths within the cell. 
Notice that the shorter wavelengths create carriers in 
the p region and the longer wavelengths create carriers 
in the n region. Since quite a large number are created 
in the n region, the diffusion length must be between 
10° and 10 * em for good efficiency. 


ELECTRICAL CHARACTERISTICS 


The electrical characteristic of a photovoltaic cell 
can be understood from Fig. 5. Above, in Fig. 4, is 
the simplified equivalent circuit showing a constant 
current generator, a nonlinear junction impedance, and 
a load. Neglected is a shunting capacitance and _ re- 
sistance and a series resistance. In operation, light 
causes a current, 7, to flow in the load which is the 
difference between the generated short circuit current, 
T,, and what flows in the non-linear junction, J; . 

T; is the usual junction current given by 


IT; = — 1) {1] 
where /, is the dark or reverse saturation current’ and 
\ = e kT, k being the Stefan Boltzmann constant and 
T the absolute temperature. V is the voltage applied 
to the junction. Since 

=I, 
I = I], — — 1) 


The maximum voltage, Vis is given when / 
From Equation [3] one obtains 


Fig. 5 shows how the current voltage characteristic 
of a solar cell comes about. The dark current charac- 
teristic is given by Equation [1]. In the light the whole 
curve is given by Equation [3] as shown. Note that /, 
can be many times the dark junction reverse current. 
The portion of the curve in the fourth quadrant is the 
region of power generation. The maximum power, 
Pox , delivered is represented by the largest rectangle 
that can be fitted into this area, and here can be seen 
one of the advantages of this type of generator. It can 
deliver almost 90 per cent of the power generated to a 
load, whereas in normal-type batteries 50 per cent of 
the power is so available. This was first pointed out by 
Lehovec in 1948." Therefore, the squarer the 7-V 
characteristic, the higher the efficiency will be. It is 
interesting to note that in a photovoltaic junction the 
current flows in the reverse direction and the voltage 
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Fig. 6—Vimax and J, as a function of light intensity. 


is in the forward direction for the usual p-n junction. 
Under conditions of low light level where /,/7, < 10 or 
where the series resistance becomes appreciable, the 
?-V characteristic can become linear. Examples of this 
will be shown later. 

Fig. 6 shows how the voltage and current varies as a 
function of the flux density in silicon cells. Note that 
the current is directly proportional to the flux, while 
the voltage is dependent on the log of the flux accord- 
ing to Equation [4]. Note that the voltage characteristic 
shows a saturation value of voltage at high flux densi- 
ties, which is quite convenient for charging storage 
batteries. Since the voltage is insensitive to the light 
level, charging can be maintained even though the 
light fluctuates considerably. To the right is shown 
the power as a function of voltage on the cell. 


THEORETICAL CONSIDERATIONS 


The expression for the maximum solar conversion 
efficiency of a solar photovoltaic cell is as follows: 


max AV ap N pn 


where K is a constant depending on the reflection and 
transmission coefficients and the collection efficiency, 
Vip is the voltage delivered at maximum power, 
Np» (Ee) is the number of photons that generate 
electron-hole pairs in the semiconductor of energy gap 
E,, and N,, £., is the input power where N,,, is the 
number of incident photons and £,, is their average 
energy in electron volts. 

A number of interesting facts can be obtained from 


Equation [5], since AV,,, is usually much greater than 
1, and K ~ 1, the equation can be written as follows: 
(6] 

LV ph 

An estimate of the maximum theoretical efficiency 
for silicon can be obtained directly, since n,, is about 
and is about yielding about 22 per 
cent. Where monochromatic light is used with energy 
equal to the band gap, the efficiency can be well over 
75 per cent since then n»,, = Np», and Vp», ~ 0.75 E,, . 
Also, since < Vmax « InJ,, 

max In J, . (7] 
Therefore at higher light intensities the maximum 
efficiency should increase. 

Equation [6] also shows how max changes as the 
semiconductor energy gap is varied. Two effects take 
place: first the number of photons absorbed n,,(2¢) 
decreases as the band gap is increased. This is seen in 
Fig. 7, where the number of photons absorbed is 
plotted against the semiconductor energy gap. Second, 
the voltage increases as the band gap increases due to 
the reduction in J, the dark current.'® Fig. 8 shows the 
result of a calculation including these effects which 
shows how nmax Varies with EF, ."* The upper curve is 
calculated using standard p-n junction theory, and the 
lower curve shows the result if the junction reverse 
saturation departs from its theoretically predicted 
dependence on Hg. Note that these curves are for the 
solar spectrum that exists above the atmosphere. The 
optimum band gap is seen to be 1.6 ev; at that point 
the nmax is 20 per cent better than what it is for silicon. 
Note that in the lower curve there is little difference 
between materials. 
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Fic. 7—Absorbtion of solar photons as a function of energy 
gap. 
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Fic. 8—Theoretical efficiency vs. semiconductor band gap. 
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9—Comparison of silicon current voltage characteristics 

in sunlight and in artificial light. 
EFFICIENCY MEASUREMENTS, CALIBRATION 

AND SILICON RESULTS 

Fig. 9 shows the 7-V characteristic for a_ silicon 
photovoltaic cell. Note that for convenience the first 
quadrant is used for this plot instead of the fourth. 
This cell has an efficiency of 7.8 per cent, as measured 
in sunlight with a Vax of 0.52 v and an J, of 12 ma. 
Note the rectangular characteristic of this curve. Also 
shown is the comparison of the cell’s operation in sun- 
light and in artificial light. The purpose of this com- 
parison is to show how efficiency is measured with an 
artificial light source. The major problem here is that 
one cannot easily duplicate the solar spectrum with an 
artificial light source and that the efficiency of a solar 
cell depends on the spectrum used. However, if one 


reproduces with an artificial source the same number of 
carriers or short circuit current in a solar cell, as the 
sun would produce, one can then say that the artificial 
source is effectively at the same solar intensity that 
one would measure with an absolute power measuring 
device such as a pyrheliometer. In the figure, the sun- 
light intensity was measured by the pyrheliometer to 
be 72.6 mw. The artificial light source was then ad- 
justed to produce the same short circuit current, and 
as can be seen from the figure, the 7-V characteristic 
was duplicated within a few per cent of the charac- 
teristic obtained in sunlight. Small deviations might be 
expected from cell to cell if the diffusion length varies 
by very much. This serves as a convenient laboratory 
method of measuring efficiency. The same technique 
can be used for other semiconductors; however, a 
separate light calibration is required for each material. 

Presently available silicon solar energy converters 
have vielded efficiencies as high as 14 per cent; how- 
ever, in production quantitites, efficiencies be- 
tween 7 per cent and 10 per cent. Such cells can be 
connected in series or parallel to increase the voltage 
or current. These cells are the easiest of all semicon- 
ductor devices to make, the major cost being that of 
producing the single crystal silicon wafers. Assuming a 
10 per cent efficiency and a solar radiation flux of 
100 mw per sq em (the average sunlight intensity in 
the northeastern part of the United States during 
the summer), a 1l-kw solar power supply would be 
10 sq m in area, cost $200,000, and weigh 20 Ib. 


EXPERIMENTAL RESULTS ON GaAs, InP 
AND CdTe 

Shown on Fig. 8 are three semiconductors, which, 
from the viewpoint of availability and stability, were 
selected to verify experimentally the theory described 
above. These materials are the III-V compound semi- 
conductors GaAs and InP and the II-VI compound 
CdTe. 

Two techniques were used to make p-n junctions on 
these new materials, namely, the diffusion process and 
the alloy process. These are similar to techniques 
employed with germanium and silicon. However, the 
scheme for doping the materials n or p-type is some- 
what different. In general, a Column II element in the 
periodic table converts GaAs or InP into a p-type 
semiconductor, while a Column VI element converts 
it to n-type. For CdTe a Column I or VI element 
yields p-type and a Column VII or II element yields 
n-type. Fig. 10 shows the construction of a typical 
diffusion cell. 

The major problems encountered in making efficient 
solar converters from these new materials has been 
the following: 

(1) Making suitable ohmic contact to the n and p 
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Fic. 10—Experimental photovoltaic cell. 
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Fic. 11—KEffect of contacts on 7-V characteristics of GaAs 
solar cell. 


regions. Rectifying or high resistance contacts lower 
the output voltage and current appreciably. In general, 
metallic plating and soldering or an indium alloy dot 
solves the problem. Fig. 11 shows the difference in 7-V 
characteristics when poor contacts are improved in 
GaAs. 

(2) Low resistivity in the p- and n-type regions. 
This is required to obtain the maximum voltage per 
junction and a low series resistance. This problem is 
solved by using highly doped original material and by 
controlling the alloy or diffusion process to get a 
maximum impurity concentration. 

(3) The fact that the p-n junctions do not behave 
according to theory. This can be caused by poor crystal 
perfection, current leakage across the junction, or 
poor ohmic contact. The solution to this problem lies 
in better material, proper etching, lower diffusion 
temperature, and improved ohmic contacts. 

Solar cells were successfully made on GaAs, InP 
and CdTe. Fig. 12 shows typical 7-V plots in an artificial 
light equivalent to sunlight of about 100 mw per sq cm 
for a GaAs cells. The curves were obtained by varying 
the load resistance and measuring the voltage and 
current while the units were illuminated. The units 
were made by baking n-type GaAs in cadmium vapor 
for four hours (G34 and G41) and 40 hours (G49). 


Note that the open circuit voltage is as high as 0.9 v. 
This is a consequence of the large energy gap of GaAs 
and should be compared to a maximum voltage of 
0.55 v in silicon. The impedance of these units at 
maximum power is about 100 ohms. The best efficiency 
yet obtained is 6.5 per cent in GaAs in small areas.” 
For cells about 1 sq cm in area efficiencies of 5 per 
cent have been obtained. 

Fig. 13 shows the results obtained on InP. These 
units were made by the alloy process where a zinc 
indium dot was fired on to the InP at 500°C for 10 
minutes. The zine converts the n-type InP to p-type 
where it recrystallizes after having alloyed with the 
InP. The ohmie contact was alloyed indium. Because 
of the presence of the dot the exposed face was masked 
somewhat so that the effective area was about ;'5 sq 


MILLIAMPS 
REPRESENTATIVE I-V CHARACTERISTICS OF GaAs DIFFUSION 
UNITS IN BRIGHT LIGHT (FML). 
12—Current-voltage characteristics of GaAs diffusion 
cells. 
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em. The maximum efficiency of IA-18 was 2 per cent. 
Note that the curves do not depart much from lin- 
earity. This implies the existence of a large internal 
resistance, about 1200 ohms for IA-11 and 700 ohms 
for LA-18. This is probably due to poor ohmic contact 
as supported by the low short circuit current, which 
is about 1 ma compared to about 10 ma for the GaAs. 

The results for CdTe are shown in Figs. 14 and 15. 
With this material we were hampered by the lack of 
low resistivity crystals. In general, the technology of 
CdTe is behind that for GaAs and InP. Fig. 14 shows 
three current voltage curves for diffusion type june- 
tions in CdTe. C-5 and C-6 were baked in iodine vapor 
at 800°C for 4 hours and C-3 in sodium vapor at 
800°C for 4 hours. Indium was used as the ohmic con- 
tact to the n-type base wafer and silver paste to the 
diffused region. The output voltage is rather low con- 
sidering the large band gap, i.e., 1.45 ev, and the short 
circuit current is quite low—about ;,j, of that in the 
other materials—thus reflecting a very high series 
resistance. This is also evident from the straight line 
characteristic of the curves. The resistance is given by 
the slopes of these lines for C-5, 10° ohms, C-6, 2 < 10° 
ohms, and C-3, 3 X 10° ohms. 


Fig. 15 shows the results on the alloy-type CdTe 
units. The junction was prepared by alloying indium 
into p-type CdTe at 600°C for 10 minutes. Note that 
here the curves are rectangular, thus indicating con- 
siderably lower series resistance. The maximum con- 
version efficiency in CdTe has been achieved with the 
alloy-type junction and is 2.1 per cent. 
RESULTS ON CdS AT WADC 
The work on CdS solar cells under the general direc- 
tion of D. C. Reynolds at WADC is quite interesting.” 
According to the theory CdS should be low on the 
efficiency scale. Because of its large band gap, 2.4 ev, 
only about 4 of the solar photons should be absorbed. 
However, spectral response measurements show that a 
large absorption takes place in the region considerably 
beyond the cut-off at 0.8 ev where the crystal should 
be transparent, thus accounting for the high efficiencies 
which Reynolds reports to be up to 7.5 per cent.” 
This device does not obey the usual theory and to 
date is not understood. It may be possible that some 
ternary compound is formed near one of the contacts 
which has a considerably smaller band gap than Cds, 
thus accounting for the response at lower energy. 
teynolds finds extremely long ‘diffusion lengths’’ in 
Cds. He finds the optimum thickness of a cell is about 
0.9 mm and suggests this is due to exciton migration 
or a luminescence effect. A third explanation may be in 
the possibility of the formation of impurity bands in 
the heavily doped Cds. 

A brief description of these cells is in order. CdS 
crystals grown by the vapor phase method and made 
n-type by indium doping to about 10° @ em are cut to 
about 1 x 1 x 0.09 em. The surfaces are roughened and 
copper is plated onto one face and diffused into the 
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I-V CURVES FOR CdS AT VARIOUS TEMPERATURES 


Fig. 16—Current-voltage curves for CdS at various tem- 
peratures. 
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interior at 460°C for less than one minute, making ¢ 
p-type layer between 1 and 10 yw thick. Contact to this 
layer is made with a Hg-loaded expoxy resin. Contact 
to the n-region is made by soldering after zine and 
copper plating. Note that this is a back-wall-type cell 
as contrasted to the front-wall cells made of silicon 
and GaAs, where the barrier is at the front surface. 

Fig. 16 shows a series of 7-V curves as a function of 
temperature on a CdS cell kindly supplied by Rey- 
nolds. This cell had an over-all nmuax of 4 per cent as 
measured in sunlight. Note the excellent square shape 
of the 7-V curve. The Vines is 490 mv and J, is 4.5 ma. 
It could be expected that Vinax should be well over 1 v 
in aé/4 per cent efficient CdS cell according to the usual 
theory. One disappointing characteristic here is the 
rapid fall off of 7, with temperature compared to the 
almost constant /, of GaAs and silicon. Because of 
this and the usual decrease in Vynax of about 2 mv 
per °C, the CdS max falls faster with temperature than 
either Si or GaAs. In fact, it appears that CdS behaves 
like a material with a band gap of 0.8 to 1.0 ev rather 
than 2.4 ev. At present there does not seem to be any 
high temperature advantage in using currently avail- 
able Cds cells. 

EFFECT OF INCREASED TEMPERATURE 

One expects that any device utilizing semicon- 
ductors would be sensitive to temperature change, and 
photovoltaic cells are no exception. Calculations were 
made to see how the maximum efficiency as a function 
of band gap varies with temperature. It was assumed 
that the short circuit current was relatively independent 
of temperature, an assumption generally borne out 
experimentally and that the main variation that one 
observes is due to the fact that Via, is a function of 
/,, a quantity quite sensitive to temperature. 

Fig. 17 shows the results of these calculations.’ 
Here, the efficiency is plotted as a function of the 
energy gap with the temperature as the parameter. 
There are two main points to be seen here, the first 


6 


one being that the optimum semiconductor energy gap 
increases as the temperature increases. Thus at room 
temperature cadmium telluride is the optimum ma- 
terial, while at 300°C the optimum material is Z with 
an energy gap of 1.8 ev. The materials of Y and Z are 
ones that may be ternary compounds; as, for example, 
combinations of GaAs and GaP. The second point to 
observe is that the relative efficiency increases at 
higher temperatures when comparing a higher band 
gap material to silicon. For example, consider the 
comparison between silicon and GaAs at 150°C where 
it can be seen that GaAs would be almost twice as 
efficient a converter as silicon. However, at room 
temperature, GaAs is not quite 20 per cent better than 
silicon. Note that on these curves, CdS, although 


somewhat better than silicon at 150°C, is still not as 
good as GaAs would be. As a result of these calcula- 
tions, GaAs is considered the most promising material 
as a high temperature (~200°C) solar converter. 
Fig. 18 shows some specific temperature measure- 
ments on silicon. This is one of the early silicon devices 
with a somewhat lower efficiency than those made 
today; however, the temperature results are typical of 
all silicon devices. Note that the maximum voltage 
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Fic. 17—Efficiency vs. energy gap for various temperatures. 
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Fig. 19—Variation of GaAs alloy junction solar cell parameters 
with temperature 
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APPLIED VOLTAGE — VOLTS 
Fig. 20—Current-voltage characteristics for CdTe films at 
varying light intensities. 


decreases with temperature at the rate of 2 mv per °C. 
The middle curve shows how the short circuit current 
remains relatively constant with temperature. The 


bottom curve shows how the maximum power de- 
creases with temperature. The slope is approximately 
0.02 mw per °C. 

Fig. 19 shows similar results for GaAs. This repre- 
sents one of the better cells. Many other cells have 
shown quite varied characteristics, many of which 
were poorer than the one shown. The voltage is seen 
to decrease at the rate of 1.7 mv per °C. The current 


is seen to be reasonably constant, and the power de- 
crease is somewhat less than that of the silicon unit. 


A NEW HIGH VOLTAGE PHOTOVOLTAIC 
EFFECT 

Since the major expense in making silicon solar 
converters is in the production of the single crystal, 
possibilities of preparing cells by deposition methods 
such as evaporation, sputtering, etc., were studied. 
During this work a very interesting effect was dis- 
covered by Pensak and later worked on by Goldstein.” 
They found that, when cadmium telluride was evap- 
orated at an angle onto a pyrex substrate at about 
100-200°C to a thickness of about one micron and 
when suitable electrical contacts were made to the 
film, under bright light illumination such films could 
produce as much as 1500 v per cm. 

Fig. 20 shows the 7-V characteristic of the CdTe 
evaporated films. Note the high impedance and high 
voltages. Decreasing the temperature increases the 
voltage to about 300 v at liquid nitrogen temperatures. 
Such films exhibit high resistance on the order of 10° 
ohms per square in the dark, and, as a result, the cur- 
rents produced by such cells are quite low. As a matter 
of fact, at present, these cells have very low energy 


FILM 


TIC COATING ELECTRODES 
ON PYREX GLASS 


TYPICAL PHOTOVOLTAIC CELL 


CELL #2 


+ 


a) 


b) 
| FERMI 


CELL #1 


Fic. 21—Representations of series-connected cells. 
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conversion efficiencies, being on the order of 10° per 
cent. However, this is a breakthrough into a new and 
interesting type of photovoltaic converter which 
warrants further research. 

Fig. 21 shows schematically what such films look 
like. At the top you see the film on glass. Contacts are 
made by using a transparent conducting coating and 
silver paste. In order to explain the extremely high 
voltages, one must assume that there are many p-n 
semiconductor barriers, possibly grain boundaries, 
that are somehow connected in a series chain. This 
would, however, not be enough to explain the results, 
since such a series connection would result in a voltage 
typical of just one barrier. The remaining barriers 
would cancel one another. In order to explain the 
present result, one must assume that the barriers are 
asymmetric. This is shown schematically at the bottom 
of the figure where a steep junction is shown between 
the n-p transition and a rather broad junction between 
the p-n transition. With light shining normally on 
such a configuration, assuming the lifetime to be small, 
it can be seen that carriers may be collected in the 
n-p region, but not from the p-n region, since those 
‘arriers generated in the p-n region would very easily 
recombine before being collected. These asymmetries 
may also be caused by impurity gradients, lifetime 
inhomogeneities, or stacking faults which are hexagon 
to cubic crystal structure transitions. Such stacking 
faults have been observed in II-VI compounds. Work 
in ZnS has shown correlations between a high voltage 
photovoltaic effect (produced by ultraviolet radiation ) 
and stacking fault structures. Voltages of 50 to 100 v 
have been observed at room temperature, and 350-400 
v at liquid nitrogen temperatures. 


PROBLEM OF IMPROVED EFFICIENCY 


There are five areas in which to seek higher efficiency 


in the single-crystal-type solar converters. Not all of 
the following areas can be optimized independently of 
the others. A good compromise is the best that can be 


expected. 

(1) Spectrum matching. Here we have the problem 
of utilizing the maximum energy of the photons ab- 
sorbed. Only for photons of band edge energy is this 
realized. At higher photon energies the excess energy 
beyond the band gap is imparted to the electron-hole 
pairs and is eventually dissipated as heat. A con- 
siderable improvement in efficiency could be realized 
if this problem were solved. Jackson has suggested” 
the possibility of using sandwich structures of differing 
band gap semiconductors. The sunlight is allowed to 
strike first the widest gap material, which absorbs part 
of the spectrum and transmits the rest to the next 
widest gap semiconductor. This material absorbs part 
of the spectrum and transmits the rest in turn to the 


next material, and so on. He calculates that for three 
semiconductors of 0.94, 1.34, and 1.91 ev band gap 
respectively a solar conversion efficiency of 69 per cent 
is possible and where ten different materials are used 
86 per cent efficiency is theoretically possible. 

Such devices would be difficult to construct mechani- 
‘ally, and each semiconductor would require a separate 
electrical connection. Series or parallel connection of 
the different semiconductors would either degrade the 
individual Vinax’s or J,’s. The optical reflection from 
ach interface would have to be reduced by special 
interference films. It has been suggested that such 
cells would be more expensive than the most valuable 
real estate, so that it would be much more economical 
to use more area and lower efficiency for most applica- 
tions. 

(2) Reflection and surface-recombination problems. 
These two problems are grouped together since it 
appears that the optimum condition exists when a 
compromise between the two is made. A surface with 
low optical reflection shows high surface recombination 
(loss of generated carriers by recombination at the 
surface) and vice versa. Prince has measured a 4 per 
cent reflection on 10 per cent efficient silicon solar 
cells whereas normal clean silicon surfaces have a 30 
per cent reflection coefficient.’ If a steep doping gradient 
or band gap gradient (i.e., GaP on GaAs) is built in 
very close to the surface even high surface recombina- 
tion or a high absorption coefficient do little harm, 
since the carriers would be accelerated inward towards 
the junction before recombination can take place. 

(3) Bulk recombination. This is a question of 
placing the junction at the optimum depth from the 
surface so that the greatest number of generated pairs 
are collected by the junction. A large diffusion length 
of about 10°* em is required to collect carriers generated 
by near band gap radiation. Such diffusion lengths are 
difficult to obtain where a semiconductor is highly 
doped and where it has first to be subjected to high 
temperatures in the fabrication procedure. In general, 
where a junction is shallow (limited by high series re- 
sistance ), the carriers produced by high energy photons 
that do not get trapped by the surface are collected 
while the ones generated by longer wavelengths are 
lost due to bulk recombination. The opposite becomes 
true as the junction is made deeper. 

(4) Vinax limited to less than the band gap. The 
Vinax for good silicon cells is about > the value the- 
oretically possible for a junction formed by highly 
doped n- and p-type regions. In GaAs this voltage has 
been as high as $ of the gap. While it is not possible to 
realize the complete band gap voltage higher doping 
should yield higher efficiencies. 

(5) Series and shunt resistance. These parameters 
‘an reduce efficiency seriously, although in practice it 
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is not difficult to eliminate them entirely. Once the 
proper doping has been achieved, optimum contact 
arrangement can be employed to negate these two 
parameters. A shallow junction will add to series re- 
sistance, so that here again junction depth must be 
optimized. 

CONCLUSIONS 

In these conclusions the advantages and disad- 
vantages of present solar cells as energy conversion 
devices are summarized. 

The advantages are: 

1) Photovoltaic cells have given the highest over-all 
conversion efficiency from sunlight to electricity yet 
measured. The best value reported is 14 per cent. 

2) Such cells are easy to fabricate being one of 
the simplest of semiconductor devices in this category. 

3) Unlimited life. 

(4) High power-to-weight ratio for satellite applica- 
tions yielding 20 lbs per kw on the ground where 10 
per cent cells are used in sunlight of 100 mw per sq 
em and about 80 lbs per kw in satellite applications 
where 10 per cent cells are used in sunlight of 135 mw 
per sq cm but where tumbling, rotation, cell matching, 
higher temperature, light losses in cover, etec., are 
taken into account. A cell 15 mils thick with contacts 
is considered in the above. 

5) Photovoltaic cells are area devices that don’t 
require high temperatures for high efficiencies, thus 
eliminating lenses and concentrators in all but very 
special applications. 

6) Where large scale application will one day be 
achieved the problems of power distribution by wires 
will be eliminated by using solar cells at the site where 
the power is required. 

The two main disadvantages to date are the high 


cost and the need for a storage device. 


In the near future silicon solar cells should be avail- 
able with efficiencies in the 15 per cent range. Cells 
using other materials or combination of materials 
should be capable of 20 per cent conversion efficiency 
and large area polycrystalline films in the square foot 
range should yield 5 per cent efficiency. Such devices 
will represent very strong competition in the solar 
energy conversion field for the thermoelectric or 
thermionic devices which are showing promise. 
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Selective Radiation Coatings. 
Preparation and High Temperature Stability 


By Panos Kokoropoulos, Ehab Salam, and Farrington Daniels 


Solar Energy Laboratory and Department of Chemistry, University of Wisconsin 


Selective radiation coatings suitable for opera- 
tion at high temperatures with focused sunlight 
are described. Coatings of CuO and Co,0, on 
polished nickel, silver, and platinum have high 
absorptivity for solar radiation and low emissivity 
in the infrared when heated. They are stable at 
at high temperatures. 


Selective radiation coatings are useful for heating by 
solar radiation. The development of efficient thermo- 
electric converters and proposals for utilizing solar 
energy in space vehicles call for increased research 
efforts in the development of surfaces which will ab- 
sorb most of the sun’s radiation (from 0.3 to 2.5 w) and 
which will emit only a small fraction of the farther 
infrared radiation emitted by a heated black body. It 
is necessary that these surfaces withstand extended 
operation at high temperatures, above 500°C in boilers 
and thermoelectric generators, as well as operation in 
the neighborhood of 100°C. 

Tabor! has shown that selective radiation surfaces 
can be made by placing a thin black coating on a bright 
polished metallic surface. Hottel and Unger®* > have 
described an excellent selective surface obtained by 
spraying an aqueous solution of cupric nitrate on 
polished aluminum under controlled conditions and then 
heating to give a film of cupric oxide. The preparation, 
emissivity in the infrared and solar absorptivity of 
several selective radiation surfaces has been reported,’ 
including an oxidized and chlorinated surface on copper, 
manganese dioxide on stainless steel, molydenum oxide 
compounds on aluminum, and other special coatings. 

The object of the present investigation’ was to pre- 
pare selective coatings by electroplating copper, cobalt, 
and other suitable metals on polished, reflecting metallic 
surfaces and oxidizing them by heating in air; to 
measure their emissivities and their absorptivities for 
solar radiation, and to investigate their stability at 
high temperatures. 

It is planned to continue this work in an attempt to 
contribute to the theories involved and to find still 
better selective coatings. 


EXPERIMENTAL PROCEDURES 


Circular plates of metal 7.5 cm in diameter and 0.1 
to 1 mm thick were treated with acid, highly polished 
with rouge and cleaned with carbon tetrachloride to 
remove grease until water gave a uniform film on the 
surface. They were then electroplated to give layers 
10-° to 10-* em in thickness. The thickness of the 
plated film was determined in two ways—by measuring 
the electric current and the time of plating, and by 
direct weighing on an analytical balance before and 
after plating. It was assumed that the density of the 
oxide film is the same as the bulk density of the oxide 
and the roughness of the surfaces has been neglected 
in these calculations. 

Copper was electroplated from a solution of 2.25 per 
cent CuCN, 3.4 per cent NaCN and 1.5 per cent NasCO; 
at a current density of about 0.003-0.015 amperes per 
sq em. Cobalt was electroplated from a hot solution of 
CoSO,;:7H.O 50 per cent, NaCl 1.5 per cent and 
H;BO; 4.5 per cent at 40°C and at a current density of 
about 0.05 amperes per sq cm. 

Extensive experiments were carried out with black 
coatings produced by dipping copper into a concen- 
trated solution of 500 gm of sodium chlorite and 250 
gm of sodium hydroxide per liter of water at 100°C. 
However, the work was discontinued because although 
the coating is adherent it is actually a mixture of oxides 
and chloride compounds. 

The absorptivity of solar radiation was determined 
with an apparatus similar to the one described by 
Tabor.® It consists of a tube 9 cm in diameter and 16 em 
long on a heliometer base, which can be rotated quickly 
so that either end may be pointed directly at the sun. 
The metal disc at the center has on one side the selec- 
tive radiation coating and on the other side a coating of 
lamp black (camphor smoke) assumed to have an 
absorptivity of 0.95. The tube is insulated and each end 
is covered with a quartz window and a removable cap. 
A sensitive thermocouple is attached to the edge of 
the disc so that it is not directly radiated by the sun, 
and is connected to a recording potentiometer. 

Absorptivity measurements were made only on 
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Fic. 1—Emissivity measurement apparatus. 
A.—heating thermostat 
B.—radiating sample 
C.—radiometer 
1).—cooling svstem 


cloudless days, and a succession of cycles 150 seconds 
long was recorded, alternating between the lamp black 
coating and the selective radiation coating. The absorp- 
tivities were calculated from the temperature differ- 
ences as recorded for the two sides using formulas 
described by Tabor® and Guild.® 

The relative emissivities were determined on the same 
coated discs that were used for the absorptivity ex- 
periments. They were measured with an instrument 
similar to that described by Taylor et al,’ as shown in 
Fig. 1. The coated metal disc is placed at one end of an 
insulated aluminum cylinder and thermostated at 110°, 
140°C, or other temperature. A radiometer of 110 copper 
constantan junctions was placed 5 em away at the 
other end of the cylinder and protected with radiation 
screens. It was made under the direction of Professor 
VY. E. Suomi of the Department of Meteorology of the 
University of Wisconsin from a design described by 
Heilman.* The recorded radiometer readings were con- 
verted into emissivities by a calibration procedure in 
which a radiation emitter of bright aluminum foil, 
assumed to have an emissivity of 0.03, was covered 
tightly with sectors of black paper of various areas, 
assumed to have an emissivity of 0.95. The ratio of 
black and polished aluminum areas permitted a setting 
which agreed roughly with the emissivity recorded for 


the selective radiation coating, and the relative radia- 


tion received was compared with that from a known 
radiator. The instrument was further checked with the 


TaBLe I—Emissivity of Selective Radiation Coatings 


Per cent After 
CuO thickness thickness U nderly black body | second 
metal 
140°C heating 


2.4 10-5 cm 
2:3 


9 


6.8 K em 
12 Ag .29 
19 


(Flat black paint, 0.95) 


emissivity of polished iron at 0.52 and with the emis- 
sivity of polished copper, nickel, and brass. 


RESULTS 


Imissivity data are summarized in Table I and Fig. 2. 
It is seen in Fig. 2 that the emissivity of the surface 
coated with thin layers of CuO or Co;O, is directly 
proportional to the logarithm of the thickness of the 
layer. 

The higher emissivity of cupric oxide on nickel than 
on silver at high temperatures is due to the fact that the 
nickel surface becomes oxidized and roughened by 
oxygen which diffuses from the air through the oxide 
coating. This explanation was confirmed as indicated in 
Table I by finding that a second period of heating 
increased the emissivity still more. 

The measurements of solar absorptivity were less 
accurate. A coating of CuO (2.3 X 10-° cm) on polished 
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Fig. 2—Emissivity data. 
A.—emissivity vs. thickness 
B.— emissivity vs. log thickness 
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silver gave an absorptivity of 0.76, as much as that of a 
perfect blackbody absorber, a coating of CuO (2.4 X 
10-° em) on nickel gave an absorptivity of 0.81 and a 
coating of Co;0,4 (6.8 & 10-° em) on silver gave 0.90. 
In general the absorptivity increases with the thickness. 
With very thin coatings, the reflectivity of the under- 
lying metal is an important factor in the absorptivity 
of solar radiation, but with thick coatings it is less 
important. 

At temperatures in the neighborhood of 100°C, at- 
tainable with flat-plate collectors, many black coatings 
are stable and suitable for extended operation. At tem- 
peratures in the neighborhood of 500°C and above, only 
a few coatings are permanent. Black paint, organic 
matter, and sulfides will all be removed by oxidation. 
The black oxides are best, but at very high tempera- 
tures they will dissociate, and if there is a chance of 
removing the metal or the gaseous oxygen, they will 
not be permanent. Again, at high temperatures some 
metals will reerystallize and change the character of 
the surface and at still higher temperatures they will 
melt. 

The results on the stabilities of coatings are sum- 
marized in Table II. 


Tasie I]—Stability of Selective Coatings in Air 


Under- 
Coating Thickness lying Temperature °C 
metal 


Properties 


Above 600 
Above 800° 


Ag crystallizes 


CuO 3X Ag 
Cu alloys with 


Ag 
CuO 3 Pt 600 Excellent (at 
least 3 hr.) 
800° Cu alloys with 
Pt 
Co,0; | 6 Ag 600° Excellent coat- 
19 Ag 800° - ing but Ag 
900°} crystallizes 
Co;0,4 | 7 Pt 1000° Excellent 


The emissivities for cobalt oxide films were consider- 
ably higher than those for cupric oxide, but part of the 
difference is due to the fact that the cobalt oxide coat- 
ings are thicker. 

The nature of the surface of the underlying metal is 
extremely important in keeping the emissivities low. 
It was found essential to give the metals a high polish. 
In the case of chemically produced black coatings on 
copper or brass it was found possible to increase the 
emissivity of the selective radiation surface from 20 
per cent to more than 40 per cent by sand-blasting 
the polished surface. 

In the case of cupric oxide on silver the coating dis- 
appeared after 17 hours at 700°C. Similar results were 
determined with cupric oxide on platinum. The dis- 
sociation of CuO at this temperature is slight, but the 


metallic copper released by the dissociation appears to 
become alloyed with the metal and removed, so that 
the equilibrium is upset and the dissociation becomes 
complete. Cobalt oxide gave superior stability and did 
not alloy with either silver or platinum. After 12 hours 
at 900°C, the Co;0, on silver was in excellent condition, 
but the underlying silver contained many small crys- 
tals. Cobalt oxide on platinum was heated at 1100°C 
for 26 hours and showed no alloying nor loss of black- 
ness. 

Various combinations of metals were tried. Copper is 
very unsatisfactory as an underlying metal at elevated 
temperatures because the oxygen of the atmosphere 
diffuses through the oxide film and attacks the bright 
copper surface, making the oxide coating thicker, so 
that it loses much of its selective radiation property. 
Eventually it becomes so thick that it peels off me- 
chanically. Nickel oxidizes at still higher temperatures, 
but silver and platinum are inert toward oxygen, 
although oxygen dissolves in hot silver. Satisfactory 
results were obtained at 500° by electroplating nickel 
with silver, electroplating a thin coating of copper on 
the silver, and then oxidizing the outer copper layer. 
Unsatisfactory results were obtained when copper was 
plated with silver and then with a thin coating of 
copper, because atmospheric oxygen diffused through 
both the copper oxide and the silver layer and formed 
a coating of oxide on the underlying copper which then 
peeled off. 

Thin coatings of copper oxide and cobalt oxide made 
from the electroplated films as described here were 
remarkably adherent and could not be easily scratched 
off. If the coatings were too thick, however, they were 
not mechanically satisfactory. Microscopic photo- 
graphs showed the coatings to be finely crystalline, but 
extremely uniform. If the oxide exists in different crys- 
talline forms there is a chance for instability over a tem- 
perature range. No such difficulty was detected with 
cobalt oxide, but nickel oxide with its green and black 
oxides was not studied for this reason. 

Qualitative experiments were carried out with the 
oxide-chloride coating on copper. Reflectivity measure- 
ments were made at different wavelengths throughout 
the infrared with a Perkin and Elmer infrared spectrom- 
eter provided with a 45° prism reflector which inter- 
cepted the light beam from the source and reflected it 
at right angles to the coated copper surface, which in 
turn reflected it back into the spectrometer. Compari- 
son of the intensity of radiation at each wavelength for 
the coated copper and a front-silvered mirror showed 
that the reflectivity is less with thicker coatings and 
that the specular reflectivity of polished copper is 
reduced (i.e., the absorption is increased) by the copper 
oxide coatings much more for the short wavelengths, 
around 3 uw, than for the longer wavelengths, around 9 u. 
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These measurements support the direct measurements 
of absorptivity and emissivity previously recorded. 

In another experiment, coated copper tubes, 18 cm 
long and 2.5 em in diameter, heated to 250°C, were 
used as sources of radiation, and the radiation intensity 
was measured at different wavelengths with the infrared 
spectrometer, using a wide (2-mm) slit. The total area 
under the recorded emission curve was taken as a 
relative measure of the heat radiated by the tubes. 
Taking the copper surface coated with lamp black as 
unity, highly polished copper gave 0.10 and coated 
copper gave from 0.23 to 0.68, depending on the thick- 
ness of the black coating. 


DISCUSSION 


Many combinations of thin coatings of semicon- 
ductors on highly reflecting metals can give surfaces 
which absorb more than 85 per cent of the solar radia- 
tion at temperatures somewhat above 100°C and emit 
less than 20 per cent of the infrared radiation emitted 
by an ideal blackbody radiator. 

For good selective radiation properties (i.e., high 
absorption in the solar region (0.3 to 2.5 w) and high 
transparency and low emissivity in infrared in the 
general region of 5 to 10 w) the coating must have a 
thickness (10-° — 10-* em) which is about equal to the 
wavelength of the radiation to be absorved and about 
fs or less of the infrared radiation which is emitted by 
the heated surface. 

The thin coating must be a black semiconductor. 
There is a close relation between the electrical conduct- 
ance and the optical properties which in turn are deter- 
mined by the atomic number and the unfilled electronic 
shells. The oxides of the transition elements are suitable 
for such selective radiation coatings. The greater the 
electrical conductance, the lower is the emissivity. A 
thin film of cupric oxide has a lower emissivity than a 
thin film of cobalt oxide. The electrical properties of 
many of the semiconductors are well known, but the 
optical properties and the emissivity of thin films are 
less well known. The semiconductors, which increase in 
conductance as the temperature is raised, should give 
lower emissivities (in comparison with that of a perfect 
blackbody radiator) at higher temperatures. 

Another factor which needs further study is the role 
of extra dissolved oxygen in the conductance and the 
emissivity of the semiconductor. 

Another aspect of the problem which calls for further 
research is the practical production and engineering 
studies of these selective radiation surfaces. Small units 
can be placed in electrical ovens at 400°C and oxidized 
in air in a few hours. If they are too large or if they con- 
tain parts which cannot be heated to 400°C, other 
methods are needed. The heat conductance of large 


metal sheets is so good that it is difficult to bring the 
temperature high enough with small heating units or 
blasts of hot air. It has been found satisfactory to 
oxidize a surface electroplated with copper in air at high 
temperatures using focused sunlight from an army 
searchlight. 


CONCLUSIONS 


Selective radiation surfaces for operation in the 
neighborhood of 100°C can be made from thin coatings 
of semiconductors of black oxides on several different 
polished metals. For operation at higher temperatures, 
the number is quite limited on account of chemical or 
physical changes. Copper oxide on silverplated nickel 
is probably satisfactory at temperatures up to about 
500°C. Cobalt oxide on platinum appears to be chemi- 
cally and physically stable at 1000°C and is the most 
reliable selective radiation surface yet studied at these 
high temperatures. It has a somewhat higher emissivity 
than cupric oxide at 140° and experiments to measure 
its emissivity at 1000° are now underway. The actual 
temperature reached under equilibrium conditions with 
sharply focused sunlight will of course be the most 
significant measurement. 

If the selective radiation surfaces are to be used in 
space vehicles, a further limitation may be placed on 
the stability of the oxide coatings. When the surfaces 
are heated to a high temperature in an exceedingly 
high vacuum such as exists in outer space, the few 
molecules of oxygen released by slight dissociation will 
be quickly removed from the surface and the equi- 
librium will be upset, so that the dissociation into 
metal and oxygen gas runs to completion, leaving only 
the metal and destroying the surface which absorbs the 
solar radiation. This difficulty can be overcome by 
keeping the temperature low, by using black oxides 
with very low dissociation pressures, or very slow rates 
of dissociation, or by enclosing the surface with extra 
oxygen in a gas-tight envelope of quartz or other trans- 
parent material. It is possible, too, that in a high 
vacuum the escape of oxygen, enclosed in the oxide as 
an impurity, may decrease the conductivity and in- 
crease the emissivity. 

Further investigations of the absorptivities and 
emissivities of many oxides and metals at high tem- 
peratures are needed, together with a more complete 
understanding of the theories involved. The work is 
being continued in this laboratory. 
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Optimum Storage of Heat with a Solar House 


By Edward Speyer 


Central Research Laboratory, American Machine & Foundry Company, Stamford, Connecticut 


This paper reports calculations of the economic 
feasibility of solar house heating (and cooling) 
for different localities in the United States, and 
presents results in terms of maximum cost per 
square foot of solar collector which can be afforded 
if the solar system is to be competitive. Sunny 
mountainous regions, with cold winters (large 
fuel bills to be saved) are the most favorable areas. 

Particular attention is given to the efficiency and 
capacity of energy storage and to the possibilities 
of saving summer heat for winter use. It is found 
that long-term storage is less economical than 
short-term storage. Storage costs of about a 
dollar per therm (10° Btu), or solar collectors 
built and installed for $1.50 per sq ft, would make 
solar houses commercially competitive. 


INTRODUCTION 


The two main parts of a solar house heating system 
are the solar energy collector (generally on the roof) 
and the energy storage facility. Collection and storage 
should each be considered in terms of capacity and 
efficiency. Collector capacity is the area perpendicular 
to the sun; collector efficiency is the percentage of 
incident solar energy which is stored or usefully em- 
ployed. Storage capacity is the number of Btu’s 
which can be released at night or on cloudy days; 
storage efficiency is the percentage of energy previously 
put in storage which is available after a certain period 
of storage time. 

Solar house heating systems with short-term energy 
storage (for night use, and for two or three cloudy days 
in succession) use solar energy chiefly when it is least 
available, i.e. in winter. Long-term energy storage, 
which will substantially carry over summer heat for 
winter use, means that at least 50 per cent of the energy 
put into storage must be available three months later, 
and that the quantity of stored energy must be suffi- 
cient to reduce significantly the collector area needed. 
Long-term storage has generally been rejected as too 
expensive, the reasoning being that the extra storage 
capacity is needed only a small proportion of the time. 
However, long-term storage involves not only larger 
capacity but also improved efficiency; the solar heating 
system thus benefits every day from the reduction in 


energy leakage and loss. This report investigates energy 
storage as to both capacity and duration (efficiency) 
as functions of collector area and efficiency, climatic 
conditions (suitable geographical areas in the U. 8.), 
auxiliary heating, air conditioning, and cost. The pur- 
pose is six-fold: 

(1) To determine the optimum combination of 
collector area, heat storage capacity, and storage 
efficiency for solar houses in different parts of the U.S. 

(2) To determine the effect on the above combina- 
tion when auxiliary (non-solar) heat is furnished. 

(3) To determine the cost per square foot of collector 
which would make solar house heating economically 
competitive with fossil fuels in the United States today. 

(4) To determine the effect of possible future im- 
proved energy storage methods on the economic feasi- 
bility of solar house heating. 

(5) To determine the feasibility of combining solar- 
powered air conditioning with solar house heating. 

(6) To determine where, and by how much, break- 
throughs in technology would make solar house heating 
economically feasible. 

A number of investigations of solar house heating 
have appeared in the literature.’* Four distinctive 
features of the present investigation may be noted: 

(1) The optimum storage efficiency, and in particular 
when it should be long-term or short-term, is not 
decided in advance, but calculated. 

(2) Storage capacity is taken as one of the variables 
in the calculation. 

(3) In the solar house investigated in this paper, 
auxiliary (non-solar) heat is occasionally furnished, but 
only the hot water heater is used. The conventional 
furnace is eliminated. 

(4) The collection efficiency is considered as varying 
with collection temperature and the collection tem- 
perature as varying with the amount of heat in storage. 


PRINCIPAL CONCLUSIONS 


(1) There is no general area in the United States 
where solar house heating is economically feasible today. 
To become competitive, one of four things must happen: 

(a) Severe rationing of coal, fuel oil, and gas. 

(b) The cost of fuel for house heating must go above 
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Fic. 1—Maximum allowable cost for each square foot of solar 
collector, assuming conventional energy storage and 
fossil fuels at 18¢/therm. 


approximately 50 cents a therm. (This is somewhat less 
than the cost of electric heat today; a therm is 10° Btu.) 

(c) An 80 per cent to 90 per cent reduction in the 
present cost of long-term energy storage of between 
300 and 700 therms. 

(d) Development of methods of building and install- 
ing solar energy collectors of 50 per cent efficiency under 
average conditions (See Fig. 3) for $1.50 per sq. ft. 
(See Fig. 1.) 

(2) From an economic point of view, the most favor- 
able regions for solar house heating, fuel costs being 
assumed equal, are mountainous regions and high 
plateaus (the Appalachian highlands and the Southwest 
plains). These regions have cold winter temperatures 
(large fuel bills to be saved) and large amounts of 
winter sunshine (relatively small collector area re- 
quired). The least favorable regions are the Pacific 
Northwest, the Mississippi Valley, and the Gulf Coast 
(See Fig. 1). 

(3) The relative feasibility of solar house heating in 
different parts of the United States varies by a factor 
of about three. That is, the maximum price that can be 
afforded per square foot of collector in the Southwest 
is three times that which can be afforded in the Pacific 
Northwest and in southern Florida (See Fig. 1). 

(4) The collector area required to provide sufficient 
energy to heat a one-family dwelling varies from 450 
sq ft in the Southwest to 2,000 sq ft in the North. 
This precludes the use of roof collectors for apartment 
buildings of more than two or three stories; thus, solar 
house heating is impractical in the heavily populated 
sections of cities. 

(5) The required collector area depends strongly on 
the amount of winter insolation. A cloudy winter with 
20 per cent less than average insolation, but average 


temperature, may increase the required collector area 
by more than a third. In order to provide for year-to- 
year weather variation without greatly increasing costs, 


an auxiliary heating unit should be used to make up 
the deficit during the less favorable winters. 

(6) Although the stored energy may not reach its 
minimum until January or February, November and 
December are the critical months because they repre- 
sent the greatest drain in stored energy. On the other 
hand, during the five or six warmest months of the year, 
the system will not be able to use, i.e., it will throw 
away, a considerable portion of the available solar 


Fic. 2—Maximum allowable cost for each square foot of solar 
collector, assuming long-term energy storage cost at one 
dollar per therm and fossil fuels at 18¢/therm. 


EFFICIENCY (%) 


TEMP. OF COLLECTION 


Fic. 3—Efficiency of collector and one glass plate vs. collection 
temperature for various ambient temperatures. 
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TEMPERATURE OF COLLECTION (°F) 
Fic. 4—Efficiency of collector and two glass plates vs. collec 
tion temperature for various ambient temperatures. 


> 
i= 
- 
. 
jo 
& 
a 
S59 
90 
FACE 
eo} TABOR SURFACE ~ ~~ 
20}— 
| 
| 
~ 
25 


energy, unless the storage efficiency is very high (and 


expensive). 

(7) The major determinant of the engineering require- 
ments for a solar house in a given locality is the effi- 
ciency of the collector. The collector is more efficient 
when it has only to heat water (or air) to a moderate 
temperature than when it must heat a fluid to nearly 
the boiling point of water. (See Fig. 3 and 4.) 

(8) If long-term energy storage is substituted for 
short-term storage, surprisingly little reduction in the 
required collector area is achieved. In some localities, 
the reduction is significant if storage capacity is in- 
creased, but nowhere is there much reduction as the 
efficiency is increased. (In some cases, because of a drop 
in collector efficiency, the required collector area may 
actually increase with storage efficiency.) In all locali- 
ties, solar houses with short-term storage are more 
economical, even when the daily saving of heat with 
long-term storage is taken into account. 

(9) If the energy storage is hot water, the most effi- 
cient maximum storage temperature is about 190°F. 
This assumes a Tabor surface collector similar to that 
of Fig. 3. 

(10) The major cost of a solar heating system is the 
collector; however, even if there is no significant im- 
provement in collector design or manufacture, a major 
breakthrough in the cost of high efficiency energy 
storage might make solar house heating economically 
feasible. Something like a 90 per cent reduction im 
energy storage costs for storage capacities between 
300 and 700 therms would be required. Solar houses 
would then be designed with long-term storage instead 
of short-term storage. 

(11) A major breakthrough in the cost of long-term 
energy storage would only slightly reduce the relative 
feasibility of solar house heating in different parts of 
the country; i.e., the ratio remains slightly more than 
double from the most to the least favorable areas 
Moreover, the most favorable areas remain approxi- 
mately the same, i.e., the Southwest and the Appa- 
lachian highlands. (Compare Fig. 1 and 2.) 

(12) The maximum price that can be afforded per 
square foot of solar collector varies linearly with the 
cost of fossil fuel; that is, if the cost of fuel oil doubles, 
twice as much can be paid for the solar collector. 
(See Fig. 12.) 

(13) If solar house heating as a whole is not eco- 
nomical, calculations will indicate that the greater the 
proportion of auxiliary heat furnished, the cheaper the 
over-all system. (This assumes, of course, that the 
solar equipment is scaled down as the demands on it 
are decreased.) The cheapest system will consist 
entirely of auxiliary heat, i.e., no solar heat at all. The 
optimum proportion of auxiliary heat, therefore, is 


dependent primarily on the point of view adopted. 


(In this paper most of the calculations assume that the 
auxiliary heat is to be furnished only to the domestic 
hot water drawn from faucets, not to the house itself.) 

(14) Solar air conditioning by itself does not appear 
economical; solar heating and air conditioning com- 
bined in the same house may be more nearly econo- 
mically feasible than solar heating alone, but the over- 
all prospects for solar air conditioning seem to stand or 
fall with the prospects for solar heating. (See section at 
the end of this paper.) 

Further detailed conclusions are given in a later sec- 
tion, following explanation of the calculations on which 
they are based. 


Orientation of Future Work 

The idea of utilizing solar energy for space heating 
and cooling is not new. Practically any competently 
designed system will work, but none, so far, has been 
cheap enough to market. There are, broadly speaking, 
two ways of approaching the design problem: the hand- 
book approach and the Einstein approach, i.e., intelli- 
gent use of materials, techniques, and principles which 
represent the most advanced state of the art, and an 
innovation based on a new and different concept. The 
past ten years have seen several competently designed 
solar houses built by the handbook approach and one 
new development—the silicon solar battery—embody- 
ing the Einstein approach. It does not seem likely that 
economical application of solar energy to heating and 
cooling will come about within the next 5 or 10 years 
by “successive approximations” of small improve- 
ments (and increases in the cost of fossil fuels), but only 
by an original “breakthrough”? such as a different 
method of collecting solar energy (perhaps with a new 
type of fuel cell), or a practical means of dispelling 
winter overcast. 

The difficulties encountered in applying solar energy 
economically indicate that research might more fruit- 
fully be directed toward the basic physical processes 
than toward engineering projects. Applications of solar 
energy to date have taken the energy as it comes from 
the sun, and sought to use it efficiently. Perhaps prog- 
ress can be made by modifying solar energy so that a 
transducer, like the silicon battery or an antenna, can 
use it more efficiently. 


DESCRIPTION OF HEATING SYSTEM 


The calculations are made for a house whose heating 
system is similar to that in the Massachusetts Insti- 
tute of Technology solar house.': *: * (See Appendix A.) 


A. Collection 

Solar energy is collected by fixed, flat-plate collectors 
tilted from the horizontal at 15° + the latitude, using 
water as the heat transfer medium. The absorbing 


Fe 
VOL. 
1959 
5 
‘ 
4 
26 : 
; 
q 


(blackened) surface of the collector has _ selective 
properties*: * absorbing 92 per cent of the solar energy 
incident on it and re-emitting only 14 per cent as much 
energy (infrared) as would a Planckian radiator at the 
same temperature. The collector is covered with a single 
pane of glass. The efficiency to be expected from such a 
collector is shown in Fig. 3; the calculations follow the 
methods given by Hottel and Woertz."® Their exact 
simultaneous equations, rather than their approxima- 
tion, must be used because of the selectivity of the 
absorber. Fig. 4 shows the results of similar calculations 
carried out for two glass cover plates. Further details 
are given in Reference 15. 

The collector operates only when the sun heats it to a 
suitable temperature. When the sun sets, the pump 
stops and all the water runs back down into the storage 
tank. This prevents loss of heat and/or freezing the 
pipes on winter nights. The pump works faster when the 
sun is high than when the sun is low, so that the 
quantity of water heated per unit time is the variable, 
rather than the temperature to which the water is 
heated during a single day. 


B. Storage 


Energy is stored as hot, but not boiling, water. Water 
has 10 per cent higher thermal storage capacity per unit 
volume than iron, and at least 25 per cent higher than 
rocks. The superiority of water’ is even greater when 
account is taken of the passages which must be pro- 
vided for heat transfer to and from solid materials. 
Since the storage tank must be insulated, it must be at 
least fairly tight, so that not much additional expense 
is incurred in making it water-tight. 

The water (or other liquid) should not be kept in 
vapor phase or under pressure. Pressure tanks are much 
more expensive, and the storage capacity per unit 
volume usually decreases rather than increases as the 
liquid passes to the vapor phase. Moreover, the boiling 
point must be kept low in order not to lower the collec- 
tion efficiency, or less importantly, the storage effi- 
ciency. 


C. Tank 

A. cylindrical steel tank in the house or basement 
stores several thousand gallons of water. (See Appen- 
dix B.) 


D. Insulation 


Glass wool or other conventional insulation is used 
around the storage tank. This is preferable to burying 
the tank in the ground, which can be calculated as an 
infinitely thick insulation which absorbs the heat that 
leaks out of the tank. The mathematical solution for an 
infinitely buried tank has been given.'' The temperature 
versus time curve is quite different from that for a 


conventionally insulated tank; the buried tank cools 
off much faster at first, but over a long period will hold 
the residual heat longer than the conventionally insu- 
lated tank. Poorly conducting clay, loam, or trap rock 
has a conductivity of about 4 Btu per sq ft per hr per 
°F for 1-in. thickness. A spherical 3,000-gal tank, filled 
with water at 207°F and buried in dry loam at 44°F, 
will cool down to 82°F in five weeks; the same initial 
and final temperatures, for a five week period, corre- 
spond to a mere 1.8 in. of glass wool with K factor of 
0.3 Btu per sq ft per hr per °F for 1-in. Other methods 
of insulation are discussed in Reference 15. 


E. Air Heat 

The rooms are heated by convected hot air. People 
who live in one-family houses are divided in their prefer- 
ences for hot water heat furnished via baseboard 
radiators and air heat furnished via ducts and gratings. 
In a solar-heated house, however, there is an over- 
riding consideration in favor of air heat. If a room is 
to be satisfactorily heated by a radiator of normal size, 
the water in the radiator must not fall much below 
140°F. On the other hand, if air is circulated through 
the house after being blown over hot-water coils, the 
water in the coils (circulated from the storage tank) 
‘an fall as low as 95°F and still heat the air sufficiently 
in a normal-sized heat-exchange device. If auxiliary 
heat is used for the air also, the storage tank water may 
be allowed to fall as low as 85°F, i.e., serve as first-stage 
heating of very cold air. Moving the minimum tem- 
perature from 140°F down to 85°F (using 207°F as 
maximum) increases the energy storage capacity of the 
tank by 82 per cent. Even more important, the efficiency 
of solar energy collection is improved during the most 
critical months of the year. 

A person who is willing to go to the extra expense in- 
volved in solar heating is quite likely to want air con- 
ditioning in his house too. If a house is to have air 
conditioning, it should be designed for air heat, because 
the same ducts and blowers can be used for both 
systems. 


F. Design Factors 


The heating load is assumed to be 25,000 Btu per 
degree day for a house built for an outdoor design tem- 
perature of O°F. The design temperatures used for the 
present calculations, based on data given in the ASHAE 
Annual Guide, are given in Table II. 

The correction factors for the Btu per degree day 
rating of a house with design temperature different 
from O°F are taken from the ASHAE Annual Guide 
and given in Table III. Its engineering significance is 
this: in localities with milder climates houses are gen- 
erally built with less insulation and radiator area than 
the same size and style house would have if built in a 
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TaBLe I—Boston Water Temperatures 
Oct. Nov De Jan. Feb. Mar. Apr May June July Aug. Sept. 


55 50°| 45 40° 40° 45 50° 55 60 65° 65° 60° 


Fu Annual Annual 
I ty nace 10use water 
st eat neat 
Therms Therms 
Blue Hill, Mass +4°F 42.0 S800 1693.91 319.10 
Central Mass. 0 $2.5 SOO 1685.75 320.17 
Ithaea, N. ¥ —10 $2.5 S00 1512.47) 320.17 
Madison, Wis. —11 13.1 SOO =1600.21 321.32 
Rhode Island +4 $1.7 SOO =1621.60 318.38 
Columbus, Ohio 0 39.9 SOO 1403.75 315.09 
Medford, Ore +2) 42.2 600 1591.45 319.46 
Lander, Wyo 13 $2.9 SOO 1746.74 320.96 
tapid City, 8. Dak.... —15 14.0 S00 1544.71 322.75 
N. Y. City and Long 
Island +5 $1.8 SOO =-1342.03 318.74 
Columbia, Mo —§ 38.9 S00 1193.89 313.23 
Grand Jetn.. Colo -10 39.1 800 «1267.89 313.66 
Nashville, Tenn. +7 36.1 S00 976.17 308.15 
Atlanta, Ga +20 33.8 600 989.1 300.85 
Fresno. Calif +30 36.6 600 =1064.69 309,22 
Albuquerque, N. M +S 35.0 800 1238.80 306.36 
Los Angeles, Calif +35 34.0 600 922.35 304.50 
Charleston, 8. C +20 32.8 600 690.55 302.28 
kK] Paso, Texas +18 31.8 600 890.69 300.51 
Ft. Worth, Texas +10 32.75 600 688.82. 302.28 
Phoenix, Ariz +-27 33.5 600 677.08 303.43 
TABLE III—Correction Factor (3) for Outside 
Design Temperature* 
Design Temperature —~20° F —10° F 0° +10° F'+20° F 
Correction Factor 0.778 | O.875 1.00 1.167 1.40 


* Taken from ASHAE Guide 


colder area. Thus, the number of Btu needed to heat 
the house per degree day is greater in the South than in 
the North. As can be seen in Table XXVI, fuel costs 
(or savings) may be greater in a milder locality, like 
Rhode Island or southwestern Oregon, than in a locality 
with severer winters, like South Dakota. 
G. Hot Water 

Domestic hot water load is figured at 120 gal per day, 
heated to 140°F. To compute the heat load for this 
water, however, one must know the city water tem- 
perature. City water (underground) temperatures are 
not readily available in the published literature, so 
Whillier’s figures for Boston are used with a correction 
made for difference of latitude. Boston and Blue Hill 
are taken as latitude 42° and the correction, plus or 
minus according to whether south or north, is one-half 
a degree Fahrenheit for each degree of latitude different 
from 42°. Whillier’s figures for Boston are given in 
Table I. 
H. Auxiliary Heat 

Hottel' and Whillier,? who studied houses with short- 
term storage, have shown that the annual dollar return 


in fuel savings per square foot of collector is lower if the 
collector must provide all the house heat (and hot 
water) than if some auxiliary (non-solar) heat is pro- 
vided in winter. This conclusion holds also for solar 
houses with long-term storage. 

There are additional reasons for providing auxiliary 
heat in a solar house. First, a stand-by source of heat 
eliminates the need for designing the solar heating 
system against the worst winter on record. That is, 
auxiliary heat can be used during worse-than-average 
winters. Secondly, even if the house is unoccupied dur- 
ing the fall, the solar collecting and storage system 
would have to be kept operating, if there were no 
auxiliary source of heat which could be tapped in 
December and January to make up for the lack of 
stored heat. The same consideration applies if people 
want to move into a new solar house in winter. 

The greatest dollar saving of conventional heat 
equipment is accomplished by eliminating the furnace. 
Since the auxiliary heating equipment is not planned 
to carry the full house heating load, i.e., it will provide 
hot water in winter and sometimes part of the house 
heat, a slightly oversized water heater using gas or oil 
or off-peak electricity, can suffice. Thus, two hot water 
tanks are to be provided in the solar house: a large 
energy storage tank, and a smaller tank (Whillier sug- 
gests 275 gal) for domestic hot water. Since the water in 
the solar system is circulated and stored for long periods 
it is kept completely separate from the domestic hot 
water, which should be of drinkable purity. A valve 
arrangement, however, permits filling the storage tank 
from the domestic hot water supply and another set of 
valves permits heat transfer from one water supply to 
the other. (Detailed design considerations are given by 
Whillier. ) 

The system for which the calculations have been 
made assumes that, with average weather conditions, 
the auxiliary heater provides the domestic hot water in 
December and January and half of February. Thus, 
the solar system provides all the house heat and 93 
months of domestic hot water, or about 96 per cent of 
the total annual heat load. If a calendar is prepared 
showing the temperature to be expected in the storage 
tank for each 3-day period during the first half of the 
winter, the auxiliary heater can be run whenever the 
storage tank temperature falls more than a few degrees 
below its average temperature for that date. The hot 
water heater would then be used for 12 or more hours 
a day during a cloudy winter month instead of only 
when needed to heat up the domestic hot water supply. 

On the other hand, when more winter solar energy 
than usual is collected and/or the number of degree 
days is less than usual, the auxiliary heater may be 
relieved of the necessity of furnishing any heat at all. 
The variation in winters from year to year is different 
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for different parts of the country, but as a rule of 
thumb, one can assume a 50-50 chance that the monthly 
difference in collected energy, or in needed house heat, 
essentially, degree days) will vary plus or minus at 
least 10 per cent from the average. If the period taken 
is a year instead of a month, the percentage deviation 
with a 50-50 chance drops to about 5 per cent.!- 


I. Storage Temperature Limits 

The temperature of the water used in heating the house 
is permitted to vary from a minimum of 85°F to a maxi- 
mum of 207°F. (The minimum temperature would 
probably have been better chosen at 90°F, since 85°F 
water requires either a very large or very slow heat 
transfer to the air which is to be convected into the 
rooms.) The maximum temperature was determined 
by the stipulations that steam was not wanted and 
that at least a 5° differential is needed between the col- 
lection temperature and the temperature of the water 
in the storage tank. 


SYMBOLS USED MORE THAN ONCE 


Area of solar energy collector (sq ft) 

Area of solar energy collector required when storage 
efficiency is perfect 

.Auxiliary (non-solar) heat furnished to the house (Btu) 
House heating load per month (Btu) 

Average monthly insolation (sunshine) falling on unit 
area of a horizontal collector (Btu/ft?/month) 

Subseript for 7** month (¢ = 1 is October) 

Heat conductance per unit area through insulation d in. 
thick (Btu/ft?/hr/°F) 

Storage tank exponential heat loss rate factor (see 
Equations 7 and 8) (reciprocal minutes) 

Number of days in 7*® month 

Money (dollars) 

Incremental capital investment in solar energy equip- 
ment 

Annual fuel savings 

Heat storage capacity (Btu) 

Area of storage tank (sq ft) 

Temperature (degrees F). Subscripts: a@ = ambient; 
b = beginning of month; ¢ = collection; e = end of 
month 

Volume of energy storage tank (gal) 

Heat load of domestic hot water per month (Btu) 
Amortization period in which solar equipment pays for 
itself in fuel savings (vears) 

Horizontal-to-tilt correction factor: ratio of insolation 
on south-facing surface tilted 15° plus latitude to insola- 
tion on horizontal surface, dimensionless (see Table IV) 


Correction factor to heat load/degree day for outdoor . 


design temperature different from O°F, (See Table III) 
dimensionless 

Efficiency of collection of solar energy, dimensionless 
(See Figs. 3 and 4) 

Shape factor of storage tank. Ratio of area to two-thirds 
power of the volume, dimensionless. (See Fig. 15) 


SIMPLIFIED CALCULATION 


The problem of designing a solar heating system with 
long-term storage is somewhat like that of calculating 


TaBLeE [V—Horizontal-to-Tilt Insolation Ratio 
(Based On Data From Reference 14) 


Insolation on south-facing surface tilted 15° + latitude 
a = . . 
Insolation on horizontal surface 


Lati- 
tude in| Jan. | Feb. 
Degrees 


April| May | July Aug. |Sept.| Oct. |Nov. |Dec. 


1.601.38 1.140.920. 29/1 .53)1.6 
38 .37/1.66)1.§ 
.64/2.16)2. 


30 
35 |1.77/1.48)1.21] .96) . 

40 |2.04)1.621.28) .99) .81) .73 
45 .84 .74 


a bank balance: each month certain withdrawals and 
certain deposits are made, and the balance at the end of 
the month is carried over for the calculation of the next 
month’s transactions. The difference between the heat 
collected and the heat furnished must be made up by 
(or, if it is positive, added to) the heat in the storage 
tank. 

The monthly heat collected consists of the product 
of four factors: the average monthly insolation on a 
horizontal surface (J), the horizontal-to-tilt correction 
(a) for the fact that the collector is not horizontal (see 
Table IV), the efficiency of collection 7 (see Figs. 3 and 
4), and the area of the collector (A). The monthly 
heat furnished by the solar system consists of three 
parts: the house heating load (#7), the domestic hot 
water heat (W), the heat leakage through the insulation 
of the storage tank. The house heat 7 is 25,000 times 
the number of degree days times the correction factor 8. 
The degree days are reprinted in the ASHAE Guide 
from Weather Bureau publications. The correction 
factor can be read from a graph made from Table IIT. 
The difference between the heat collected and the heat 
furnished must be divided by the heat storage capacity 
(Q) per degree to give the change in temperature 
between the beginning and end of the month (7, — 7) 
in the stored water. 

The insolation data (J) was obtained privately'® from 
the Weather Bureau. Similar data are given in Refer- 
ence 16. It may be remarked that to get its monthly 
5-year averages, the Weather Bureau averages each 
month separately, with no regard for possible correla- 
tion between months. It might turn out, for example, 
that years with unusually high summer insolation 
might generally have unusually low autumn values, 
but such correlations are lost if one simply averages 
the values for each month in successive years. It would 
be better to plot for each station a curve for each year’s 
values, month by month, and calculate the center of 
gravity or the maximum density of the curves. 

In the simplified calculation, we assume that an as- 
yet-undiscovered method of energy storage provides 
perfect storage efhciency and permits energy collection 
at a constant efficiency of 50 per cent. Since there is 
perfect storage efficiency, there is no heat leakage. 
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Moreover, since the energy can be stored for as long as 
we please, no special arrangement need be made for 
the peak load in winter. Thus there is no need to supply 
auxiliary heat (except during colder-than-average 
years), and the total annual heat collected equals the 
total annual heat furnished. For the simplified calcu- 

lation we can write, 

207 — 

The calculation is begun by assuming that the tem- 
perature of the water in the storage tank on October 
Ist is 207° = T,. We substitute into Equation [I] the 
appropriate values for October, and calculate T,, 
which then becomes the 7, for November. We continue 
on through the year, making sure that the temperature 
in the storage tank remains always between 85° and 
207°. If the temperature falls below 85°, we halt the 
calculation, assume a larger value for A, and begin with 
October again. The problem is to find the minimum 
collector area which will satisfactorily heat the house. 
The minimum value of collector area is found for each 
size of tank (heat storage capacity). As the storage 
capacity increases, a smaller collector area is required 
and less heat is thrown away in summer, until finally < 
capacity is reached beyond which no further decrease 
in collector area is possible. That is, all the available 
heat is used. The result of the calculation for 21 differ- 
ent localities in the United States is shown in Fig. 5; 
the sharp elbow in each curve represents the combina- 
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Fic. 5—Area of collector vs. storage capacity, assuming per- 
fect storage insulation and 50 per cent constant col- 
lection efficiency. 
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Fic. 6—Collector area vs. relative storage capacity for various 
efficiencies and localities. 


tion of minimum collector area with maximum useful 
storage capacity which will satisfactorily heat the 
standard house. The horizontal portion of each curve 
beyond its elbow indicates that further increase of 
storage capacity permits no concomittant reduction in 
collector area. 

Fig. 5 indicates that larger storage capacity can be 
utilized in the colder cities. To put it in terms of insola- 
tion instead of degree days, the less insolation a locality 
receives, the greater the saving in collector area that 
can be achieved with long-term storage. The curves 
with the steepest slopes—especially Medford, Lander, 
and Fresno—indicate that the ratio of the insolation 
per square foot of collector to the heat load per square 
foot of collector, falls unusually far in the late autumn 
and early winter. Medford, for example, although 
having mild winter temperatures, is in the Pacific 
Northwest where winter overcast typically blocks out 
more than 60 per cent of the possible sunshine. 

Using 85°F as the base level, 1,000 gal of water at 
207° represents a storage capacity of approximately 
1,000,000 Btu. Thus the maximum useful storage 
capacity of the most even-climated city shown in Fig. 
5—Los Angeles—requires a tank holding more than 
26,000 gal. This would fit in one-half of a full basement, 
but it would be expensive. The practicality of a design 
based on Fig. 5, therefore, depends on the economics of 
our imaginary breakthrough in energy storage methods 
(which might not, for example, require a water-tight 
tank). 

The figure of 50 per cent collection efficiency is a pure 
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guess, since we do not know the conditions under which 
the energy will be stored. If we assume other efficiencies, 
such as 30 per cent or 70 per cent, the results will be 
as shown in lig. 6. It will be noted that the maximum 
useful storage capacity is independent of the collection 
efficiency; this is due to the constant value of the 
product of collection area times collection efficiency. 
In other words, graphs of collector area versus collection 
efficiency would be a family of hyperbolas, each hyper- 
bola corresponding to a different storage capacity. 


ECONOMIC FEASIBILITY 


The solar collector and associated energy storage 
facility cost considerably more than the conventional 
heating equipment which is replaced. AM is defined 
as the difference, i.e., the incremental capital invest- 
ment in solar equipment. Interest lost on the money 
invested, plus depreciation, maintenance, and repair 
of the solar equipment, are assumed to total 10 per 
cent. Then 


AM 


where Y is the period required for the solar system to 
pay for itself, and M, is the annual fuel saving, in 
dollars. 

This formula is quite general. It applies to house 
heating or house cooling, or both, and holds whether or 
not auxiliary heat is supplied. Fig. 7 shows the amorti- 
zation period Y plotted against the incremental in- 
vestment AJ/ for different dollar values for annual 
fuel saving. Fig. 7 provides sufficient reason for reject- 
ing many plans for solar houses, such as those for pro- 
viding heat on occasional cold nights in Florida. If we 
take 10 years as the maximum acceptable amortization 
period (mortgages customarily run longer, but money- 
saving equipment usually claims less), we see from 
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Fig. 7—Amortization period vs. incremental capital invest- 
ment required for solar system, for various annual 
amounts of fuel saving. 


Equation [2] that to be economically feasible on other 
than a novelty-luxury basis, the annual fuel savings 
M; must exceed one-fifth of the incremental invest- 
ment AM, 

In fact, the equation 


AM =5M, [3] 


is a good rule-of-thumb criterion for estimating eco- 
nomic feasibility. 

The amortization period Y is not a suitable index for 
economic feasibility. We choose instead, the maximum 
allowable cost of the collector per square foot 0M/dA 
when the amortization period is 10 years. 

lig. 2 indicates the value of 0.7/0A in different 
regions of the United States, using the simplified basis 
of calculation, i.e., the data shown in Fig. 5. The 
annual fuel saving .W; is calculated as the sum of the 
monthly heats furnished times the assumed fuel heat cost 
of 18 cents a therm. Fuel costs vary considerably, of 
course, with location, being cheaper near gas fields and 
coal mines and more expensive where transportation 
costs are higher. No attempt has been made to take 
this variation into account; the ratio of 18 to the actual 
number of cents/therm should be taken as a correction 
factor in using Fig. 2 (and also Fig. 1). 

The cost of the furnace which is saved is taken as 
either S800 or $600 according to the location. (See 
Table II.) The cost of the solar heating system con- 
sists of the cost of the collector plus the cost of energy 
storage: 

aM aM {800 


AM=A +Q- 


aA aQ 600 


To calculate the data for Fig. 2, the cost of energy 
storage 0. /dQ was assumed to be one dollar per therm, 
which is only about 10 per cent of the actual present-day 
cost. We substitute Equation [4] into Equation [2]. 
read off the “elbow” values of Q and A from Fig. 5, 
and compute 0.//dQ. The data used in computing 
Fig. 2 are given in the last column of Table XXVII. 

The “elbow” figures, rather than a smaller storage 
capacity and larger collector area, represent the maxi- 
mum 0M /dA because of the extremely low storage cost 
assumed. This can be checked readily by substituting 
values from Fig. 5 into Equation [4] with AJ/7 assumed 
constant. 

No account has been taken of the cost of piping to 
and from the collector, or of the various valves and 
thermostats required. The cost of such solar equipment 
is certainly greater than the cost of the conventional 
equipment of this type, which is eliminated. This extra 
cost, which is independent of locality, is thus implicitly 
considered as part of the cost of the collectors. There is 
a slight inaccuracy introduced, because lump costs are 
considered as distributed over the collector area. The 
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result of the inaccuracy is to make solar heating in cold 
localities seem a little more nearly feasible (with large 
collector areas), and in warm localities a little less 
feasible, relative to each other, than would be the case 
if the lumped and distributed costs were accurately 
divided. 


FULL CALCULATION 


We now proceed to calculate the design data for our 
standard solar house using conventional methods of 
heat storage. We follow the same general procedure as 
in the simplified calculation except that we must also 
take into account: 

(1) The heat leakage through the insulation of the 
storage tank. 

(2) The thickness (and cost) of the insulation of the 
storage tank. 

(3) The effect of furnishing auxiliary heat. 

(4) The effect of the above three factors and the 
average outdoor daytime temperature on the collector 
efficiency. 

The basic iterative equation, giving the change in 
storage tank temperature for the 7th month is 


lan A—(H W — B). 


— 24 — Ta) 
d [5] 
(T, — T,); = K 

8.337V + 12 ao V23m, 


Equation [5| is essentially Equation [1], but includes 
an auxiliary heat term (B) and a heat leakage term, 


24K /daV™ (7: +f r,) 


The iteration starts with 7 = 1 (October) and T, 
207°. T, is the ambient temperature in the basement 
(around the storage tank), not on the roof. The tem- 
perature inside the storage tank is taken as the average 
temperature for the month (JT, + T,)/2. The factor 
8.337 is the number of pounds of water in one gallon: 
the specific heat of water, being unity, does not appear 
explicitly. 

Long-term storage was defined in the opening para- 
graph as retaining at least 50 per cent of the energy for 
a period of three months. The three-month period is 
somewhat arbitrary, but it permits heat collected in 
the sunny months of August, September, and October 
to be used during the worst deficit months. Before 
carrying out the full calculations, we inquire as to the 
practicality of reducing the conductance through the 
walls of the storage tank so as to achieve long-term 
storage. 

We first define a shape factor o for the storage tank 
as the ratio of the area of the tank to the two-thirds 
power of the volume. (See Appendix B.) 


¢ = [6] 


sol 


‘ SS FACTOR 
Fic. S—Residual storage capacity vs. heat loss factor, for 
various storage periods. 


We next define a decay factor L by the equation 
a 1 
( ‘i 7.) = ( To de \7] 


where 7 is the temperature at time ¢, and To is the 
temperature at time ¢ = 0. L is in reciprocal minutes. 

Equation {7| is the equation of a standard cooling 
curve when the heat absorbed by the insulation is neg- 
lected and the conductivity of the medium on each side 
of the insulation (water and steel inside and air outside) 
is infinite. These two simplifications tend to cancel each 
other, although for most of the situations dealt with 
here the latter is larger. This means that Equation [7] 
tends to show the tank as cooling more quickly than it 
actually would, i.e., the error is in the conservative 
direction. To calculate the cooling per month, we as- 
sume a standard-length month, i.e., 43,800 minutes. 

It is not hard to show that 

Kid 
= [8] 
19148 
where 1914 is a dimensional constant. 

The insulation required for various percentages of 
energy retention and various periods of time is most 
conveniently expressed in terms of L. Fig. 8 shows this 
relationship, assuming 7) as 206.6° and T, as 44.6°. 
To interpret the Z values in terms of tank capacity and 
insulation thickness, Table XXVIII is provided. The 
shape factor o is taken as 5.54, which applies to a 
“square cylindrical tank,” 1.e., a right circular cylinder 
whose height equals its diameter. (See Appendix B.) 
From Fig. 8 it can be seen that if 50 per cent of the 
heat stored is to be available three months after 
storage, L must be 3.75 & 10-° per minute. Table 
XXVIII shows that for 6,000-gal storage capacity this 
corresponds to a conductance of slightly more than 
0.02 Btu per square ft per hour per °F. 

Glass wool, which is a fairly typical and cheap insula- 
tion, has a conductivity of about 0.3 Btu per sq ft per 
hour per °F per in. thickness, at 190°F. This means that 
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for long-term storage the tank should be enclosed by 
at least 15 in. of insulation. 

Equation [5] cannot be solved in a straightforward 
manner, because we do not know the proper value of 
A, nor of n;. We start with a value of A obtained by 
approximation methods and from experience, and raise 
or lower it according as the calculation shows the 
minimum temperature in the storage tank goes below 
85°, or does not ever go down to 85°. Five square feet is 
taken as the minimum significant increment of collector 
area. 

Collection efficiency can be read off specially pre- 
pared graphs similar to Fig. 3, but with smaller incre- 
ments of ambient temperature. This ambient tempera- 
ture is out-doors at the collector, not in the basement. 
Since the collector is operative only during the warmest 
part of the 24-hour day, the ambient temperature is 
computed, as is done in the calculations for the M.1.T. 
solar houses,® as 79 of the mean ambient temperature 
plus ;°5 of the maximum temperature. All of these are 
daily temperatures averaged over the month. Unfor- 
tunately, the collection efficiency is different for each 
month, and in addition depends on the value chosen 
for A, since A controls the monthly change of the stor- 
age tank temperature. The collection temperature, it 
will be recalled, is 5° above the storage tank tempera- 
ture. We use the average storage tank temperature for 
the month, which requires knowing the temperature at 
the end of the month, which is just what the iterative 
formula is supposed to find out. 

The problem was programmed and run on an IBM 
650. It was necessary to provide for a calculation of 
collection efficiency, rather than a graphical determina- 
tion. This was accomplished to a satisfactorily close 
approximation (within a few tenths of a per cent) by 
regarding the curves of Fig. 3 as portions of the upper 
branches of a family of parabolas with axes parallel to 
the x-axis. 

Some difficulty was encountered in running the pro- 
gram because successive iterations did not always form 


Biue Hil! , Mass. 


Be {10,000 gais)___ — 


8, (10,000ga8) 


COLLECTOR AREA | SQUARE FEET) 


go's 


Fic. 9—Collector area vs. thermal conductivity of storage 
tank insulation, for various schedules of auxiliary heat 
(see Table XXIX) and 3,000-gal and 10,000-gal tank ca- 
pacity. 
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Fic. 10—Collector area vs. thermal conductivity of storage 
tank insulation, for 3,000-gal and 10,000-gal tank ca- 
pacity. 
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Fig. 11—Collector area vs. thermal conductivity of storage 
tank insulation for 3,000-gal and 10,000-gal tank ca- 
pacity. 


a sequence converging to vanishing terms. The trouble 
was traced to the unexpected sensitivity of the calcula- 
tion to the value of the collection efficiency, and was 
overcome by having the computer determine if the 
differences between two successive pairs of computed 
collection temperatures have the same sign; if not, the 
terms may not be approaching a limit, and it is advis- 
able to replace the third trial temperature by a tem- 
perature closer to the first temperature. In symbols, if 


<0, then a new 

TCs (9) 
Te;’) = Te; + — 


is formed. The former Tc; and Tc, become Tc,’ and 
Tc’, respectively, and a new Tc,’ is computed. 

The results of the calculations, together with the 
input data, are given in Tables IV to XXYV. Figs. 9-11 
are plotted from Tables V, XII and XX, respectively. 
These results are discussed in a later section. 

Calculations were run for one place (Blue Hill) using 
different schedules of auxiliary heat. The “regular” 
schedule (B,), it will be recalled, provides auxiliary 
heat during an average winter only for the domestic 
hot water, and then only for December, January, and 
half of February. The second schedule (B.) assumed no 


» 
— 
| 
2000 
| 
= 
“48 
ae 10,000 gais 
1600 
1500 
| 
: ° 02 0 2 2 
« 
Albuquerque , N.M. 
| 
{720 
a | 
no 
{700 
| 0,000 
= = 
680 
3.000908 
660 
| 
640 
59 
630 
| | 
620 | | 
04 06 2 4 6 8 2 
¥ 
4 
. 
1800 
}790 
200 
=) 


Insolatic 

ne 
House heat 
Water 
Degree 
Aver 

ft 
Efficiency 7 


T+ for gal; K 


gal: 


Efficiency 
T, for 


1410 ft 


000 gal: Kod 


Efficiency 
for 


ft 


10.000 gal; K 


Efficiency 7 


000 gal: 


000 gal 
1390 ft 
Efficiency 
T, for 10,000 gal 
1520 ft 


Efficiency 


therms ft 


1000 gal: Kd 


Efficiency 7 


for 000 gal; K 


ambient temp 
3000 gal; K/d 
ft? 
Efficiency 7 
for 
1680 ft? 


Efficiency 7 


3000 gal; K/d = 0.2; 


1400 


TABLE V 


Dex 


58.10% 


122.88 


60.71% 


102.54 


63.60% 


119.02 


6§1 


TaBLe VI 


De 


0. 466 


108.00 
26.44 

432 
53.33 


207 


TaBLeE VII 


0.445 


66 


Blue Hill, Massachusetts 


Jan. Feb. March April 


0.415 


5.29 
7.00 


15. 
2 


586 
47.06 
193. 45° 


34.487 


190.38 


35.07° 
180.58 


168.33 


41.68% 


150.52 


46§.90° 


140.65 


Central Massachusetts 


Feb March April 


143.00 
27.09 

572 
47.01 

191.72° 


59.32% 45.66% 34.77 
85.40 128.05 


58.63% 


Ithaca, New York 


Feb March 


48.14°7 38.55% 
101.41 161.57 


51.48% 43.99% 


34 


174. 98° 


June 


0.424 


16.70 
24.00 
63 
67.78° 


207° 


38.58% 


207° 


38.58% 


207° 


38.58% 


207 


38.587 


June 


0.430 


7.50 
.09 


June 


0.439 


18.16 
24.09 
83 

70.98 


207° 


39.54% 


207° 


39.54% 


July 


0.481 


10.04% 


207 


10.04% 


207° 


40.04% 


40.04% 


July 


0.463 


3.00 
23.34 
12 
73.20° 


207° 


40.19% 


207° 


40.19% 


91% 


Aug. 


0.445 


6.89 
23.25 
26 
70.77° 


207° 


39.47% 


207 


39.47% 


207° 


39.47% 


207° 


39.47% 


207° 


39.47% 


207 


39.47% 


207 


207 


39.62% 


0.491 


8.75 
23.34 
40 
73.57° 


207° 


40.30% 


207° 


40.30% 


Sept. 


0.472 


32.33 

24.00 
122 

64.50° 


207° 


37.59% 


207° 


37.59% 


207 


37.59% 


207° 


37.59% 


207° 


37.59% 


207° 


37.59% 


207° 


37.59% 


207° 


37.59% 


Sept. 


0.498 


36.00 
24.09 
144 

63. 98° 


207 


37.43% 


207° 


37.43% 


Sept. 


0.493 


34.13 
24.09 
156 

66. 52° 


207° 


38.21% 


207° 


38.21% 


Oct 
‘ 
n therms ft? (tilt correct 0.417 0.336 0.339 0.376 0.365 0.425 0.437 3 
therms 103.62 183.65 286.47 316.68 284.08 234.79 in 71.29 2.12 
therms 26.35 27.00 29.45 31.00 28.00 29.45 26.35 23.25 ee ee 
39] §93 1081 1195 1072 S86 269 
nt tem} roof 53.89 $3.16 31.50 28.18 28.31 36.77 58.85 72.71 
WO gal; K/d = 0; A = 207 207 145.64 86.27 111.01 108.49 207° 207° 
34.28% 41.52% | 57.23% | 61.34% | 58.21% 47.98% | | 35.85% | | 40.01% | 
02; A = | 207 207 140.57 85. 60 108.89 106.89 (07° 207° = = 
1400 ft 
Efficiency 34.28% 42.39% | 61.75% | 58.77% 18.70% | 35.859 | 04> 
fox t/a = | 207 07 85.03 89 2.53 203.2 207 id 
T) for V = Sm = 0.1; A 207 207 5.03 101.88 102.53 3.24 7 = 
1410 ft? a 
34. 30° 45.28% 62.79% 60.40% 50.97% 36.48% 10.04% a 
= = 0.2; A 207 207 85.13 92.07 97. 48° | 190.35° | 207° = 
34. 28° 18.54% 64.17% | 62.49% | 53.63% | | 38.66% | 40.01% | 
- Sd = 0; A = 207 207 158.71 96.98 85.81 91.02 zz 203 . 87° 207° 207° iz || ss 
Ee 34.30% 39.31% 53.73% 63.37% 64.23% 57.25% 42.36% 36.35% 38.58% a a2 
T, for V = = 0.02; A 207 207 152.49 93.42 85.38 91.16 148.92° 199.09 207 | 
1310 ft 
Efficiency 7 34. 30° 40.38% 55.17% 63.93% 64.26% 57.46% 43.44% 37.18% 38.58% a 
Kd =0.1; A 207 203.54 131.57 85.56 85.63 91.94 143.34 183.35° | 207 | 
34. 89° 14.43% | 59.40% | 64.97% | 64.13% | 58.14% | | 39.83% | 38.58% | | 39.47% | : 
K/d = 0.2; A = 207 180.35 85.26 89.89 95.20 171.37 207 207° 207° 
7 
RRS 50. 12° . 22° 64.41% 63.085 58.07% 49 41.83°; 38.58% 39.47% 3 
1959 
Oct Nov = Jan May Aug. 
Insolation in (tilt correct 0.383 0.384 0.351 0.364 0.433 0.415 0.437 0.458 
incl 
House heat (therms 193.50 284.75 310.00 284.25 235.00 63.50 1 | 7.25 i. 
Water heat (therms 27.09 29.54 31.09 28.08 29.54 26.44 2 23.34 a 
Degree days 774 1139 1240 1137 940 254 70 29 a 
‘ Aver ambient temp roof 41.94 29.56 26.76 26.56 35.66 59.12 68.58 71.28 7 
T, for V = 0; 4 = 1425 207 162.13° | 115.24 85.50° 123.03 207 207 207° = 
ft 
0.2; A 207 207 137.55 121.28 207° 207 
Efficiency 34.11 12.50% 52.69° 59.61% | MM, | 47.53% | 37.74% | 35.93% | 38.82% | | | 
2 
Oct Nov De Jan April May uly Aug. 
Insolation in therms/ft? (tilt correct — 0.292 0.259 0.283 0.369 0.402 0.398 0.455 0.514 
House heat (therms 168.44 246.97 270.38 252.88 213.94 132.56 63.88 3.72 
Water heat (therms 26.44 27.09 29.54 31.09 28.08 29.54 27.09 26.44 23.34 
Degree days 451 770 1129 1236 1156 978 606 292 17 us lee aie 
vot 55.37 43.30 31.14 28.26 27.38 36.27 48.88 60.70° 75.68 
0; A 1590) 207 207 158.03 85.04 111.21 173.31 184.88 207 207 | 
34.75 39. 487 55.54% 61.54% 36.54% 36.4% 40.91% 
A= | 207 207 121.15 85. 80 164.20° | 207° 207° 
34.75% 15.59% 60.78% 62.77% | 40.08% | 36.43% | = 


TasB_e VIII—Madison, Wisconsin 


Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept. 


Insolation in therms/ft? (tilt correct 0.467 0.321 0.302 0.374 0.369 0.491 0.437 0.461 0.431 0.481 0.502 0.523 
incl.) 
House heat (therms) 99.03 192.23 280.91 307.01 260.41 217.48 124.92 58.68 17.69 2.80 6.69 32.36 
Water heat (therms) 26.54 27.18 29.64 31.19 28.17 29.64 27.18 26.54 24.18 23.44 23.44 24.18 
Degree days 459 891 1302 1423 1207 1008 579 272 82 13 31 150 
Aver. ambient temp. (roof) 53. 88° 37.31° 25.73° 21.77° 24. 48° 35.74° 48.73° 60.71° 70.85 76.57° 74.36 65.94° 
T, for V = 3000 gal; K/d = 0; A = 1570 | 207° 207° 142.65° 85. 63° 145.05° 140.60° 207° 207° 207° 207° 207° 207° 
ft? 
Efficiency 7 34.28% 40.34% 56.24% 54.82% 47.19% 40.24% 32.62% 36.43% 39.50% 41.16% 40.52% 38.03% 3 
T;, for V = 3000 gal; K/d = 0.2; A = | 207° 207° 104.02° 85. 29° 124.39° 126.66 207° 207° 207° 207° 207° 207° 
1600 ft2 
Efficiency » 34.28% 46.67% 61.83% 57.84% 52.50% 42.57% 32.62% 36.43% 39.50% 41.16% 40.52% 38.03% 
TasBLe IX—Rhode Island 
Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept. : 
Insolation in therms/ft? (tilt correct 0.481 0.405 0.382 0.375 0.407 0.481 0.436 0.457 0.439 0.489 0.480 0.528 
incl.) 
House heat (therms) 100.87 177.91 274.02 297.84 269.78 231.39 150.91 68.31 15.36 0 6.88 28.33 
Water heat (therms) 26.29 26.94 29.39 30.94 27.94 29.39 26.94 26.29 23.94 23.19 23.19 23.94 
24 Degree days 672 1035 1125 1019 874 570 258 58 0 26 107 
Aver. ambient temp. (roof) 45.15° 33. 88° 31.49° 31.03° 39. 23° 48. 85° 59. 86° 68.75° 74.18° 72. 46° 65.76° 
T, for V = 3000 gal; K/d = 0; A = 207° 161.38° 88.06° 85.77° 127.92° 183.41° 207° 207° 207° 207° 207° 
1210 ft? 
Efficiency 7 34.87% 39.47% 55.31% 65.51% 59.76% 47.19% 36.84% 36. 16% 38.87% 40.47% 39.97% 37.98% 
for V = 3000 gal; K/d = 0.2; A = 207° 207° 130.09 94.48° 85. 15° 126. 22° 166. 13° 192. 98° 207° 207° 207° 207° 
1310 ft? 
Efficiency n 34.87% 44.66% 58.99% 64. 68% 60.12% 50. 16% 42.14% 38.53% 38.87% 40.47% 39.97% 37.98% be. 
X—Columbus, Ohio 
‘ Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept. 
“ Insolation in therms/ft? (tilt correct 0.477 0.367 0.326 0.303 0.354 0.438 0.449 0.458 0.435 0.493 0.509 0.542 
incl. 
House heat (therms) 84.25 173.25 258.00 273.50 236.50 195.25 111.00 45.00 7.75 0 2.00 17.25 
Water heat (therms) 26.01 26.67 29.11 30.66 27.69 29.11 26.67 26.01 23.67 22.91 22.91 23.67 
Degree days 337 693 1032 1094 946 781 444 180 31 0 s 69 
Aver. Ambient Temp. (roof) 57.65° 44.78° 34. 16° 31.86° 33. 75° 42. 89° 53.17° 64.49° 74.21° 78.27 76.21° 69.98 
T, for V = 3000 gal; K/d = 0; A = 207° 207° 170.91 93. 96° 85. 16° 164.11 204.41° 207° 207 207° 207° 207 4 
1380 ft? 
Efficiency 7 35.47% 37.69% 53.07% 64. 88°7 55.29% 34.48% 37.59% 40. 41.65% 41.06% 39.24% 
T, for V = 3000 gal; K/d = 2; A = 207° 207° 140.04 97.10° 85. 05° 183.77° | 207° 207° 207° 207° 207° 
1500 ft2 
Efficiency 7 35.47% 42.92% 57.21% 64.44% 56.52% 43.65% 38.05% 37.59% 40.48% 41.65% 41.06% 39.24% 
TaB_Le XI—Medford, Oregon 
Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept. 
Insolation in therms/ft? (tilt correct 0.537 0.347 0.262 0.306 0.413 0.539 0.573 0.583 0.551 0.561 0.660 0.630 
incl.) 
House heat (therms) 114.10 218.40 287.70 301.70 219.45 193.20 133.35 72.45 24.15 0 0 26.95 
Water heat (therms) 26.38 27.03 29.48 31.03 28.03 29.48 27.03 26.38 24.03 23.28 23.28 24.03 
Degree days 326 624 822 862 627 552 381 207 69 0 0 77 
Aver. ambient temp. (roof) 57.56° 47.29° 39.94° 39.63° 45.84° 49.66° 56.47° 62. 69° 69.61° 76.81° 76.35° 69.68° 
T, for V = 3000 gal; K/d = 0; A = 207° 207° 177.97° 85. 16° 162.87° 207° 207° 207° 207° 207° 207° 207° 
1840 ft? 
Efficiency n 35. 45% 37.30% 54. 86° 57.04% 39.44% 32.96% 35.10% 37.04% 39.13% 41.23% 41.10% 39. 15% 5 
T;, for V = 3000 gal; K/d = 0.2; A = 207° 207° 150.27° 85.52° 151.51° 207° 207° 207° 207° 207° 207° 207° 4 
1970 ft 
Efficiency » 35. 45% 42.00% 58.94% 58.63% | 41.36% 32.96% 35.10% 37.04% | 39.13% 41.23% 41.10% 39.15% 3 
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Tas_e XII—Lander, Wyoming 


Oct Nov Dec. Jan. Feb. March April May June July Aug. Sept. 


Insolation in therms/ft? (tilt correct 0.673 0.576 0.297 0.602 0.559 0.654 0.640 0.576 0.573 0.615 0.551 0.651 
incl 

House heat (therms 132.96 220.89 290.95 314.30 248.03 219.84 144.53 83.31 34.29 1.47 4.84 51.33 

Water heat (therms 26.51 27.15 29.61 31.16 28.14 29.61 27.15 26.51 24.15 23.41 23.41 24.15 

Degree days 632 1050 1383 1494 1179 1045 687 396 163 7 23 244 

Aver. ambient temp. (roof $8.14 33.03 24.03 20. 46° 26. 65° 35. 20° 45. 40° 55. 28° 65.37° 74.84° 72.63° 61.64° 

T, for V = 3000 gal; K/d = 0; A = 207 207 207 85. 54° 207° 207° 207° 207° 207° 207° 207° 207° 
1910 ft 

Efficiency » 32.45% 27.46% 45.93% 44.97% 25.31% 28.19% 31.54% 34.73% 37.86% 40. 66% 40.02% 36.72% 

T» for 3 3000 gal; K/d = 0.02; A = 207 207 207 85.38 207° 207° 207° 207° 207° 207° 207° 207° 
1940 ft? 

Efficiency » 32.45% 27.46% 46.01% 45.01% 25.31% 28.19% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 

= 3000 gal; Kd = 0.1; A = 207 207° 207 85.38 207 205. 78° 207° 207° 207° 207° 207° 207° 


Efficiency 7 32.45% 27.46% 45.97°7 5 25.48% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 

T;, for V = 3000 gal: K/d = 0.2: A = 207 207 207 85. 32° 207 202.55° 207° 207° 207° 207° 207° 207° 
2240 ft? 

Efficiency 7 32.45% 27.46% 46.01% 45.02% 26.19% 29.02% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 

7, for V = 10,000 gal; K.d = 0; A = 207 207 195.67 85.20 190. 10° 183.90 207° 207° 207° 207° 207° 207° 
1400 ft 

Efficiency 7 32.46% 29.59% 47.81% 17.68% 32.83% 32.46% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 

T, for JV 10,000 gal; Kd = 0.02; A = 207 207 194.60 85.35 194.56 182.35 207° 207° 207° 207° 207° 207° 
1470 ft 

Efficiency 7 32. 46% 29.80% 47.96% 16.95% 32.29% 32.70% 31.54% 34.73% 37. 86% 40.66% 40.02% 36.72% 

T, for V = 10,000 gal; K d=0.1; A = 207 207 189.98 85.66 207 175.07° 207 207° 207° 207° 207° 207° 
1735 ft 

Efficiency 7 32.43% 30.62% 48.63% 44.95% 31.34% 34.02% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 

T, for V = 10,000 gal; K/d = 0.2; A = 207 207 184.04 85.13 207° 169.87 207 207° 207° 207° 207° 207° 
2025 ft 

Efficiency 7 32. 46% 31.74% 49.64% 45.05% 32.31% 34.95% 31.54% 34.73% 37.86% 40.66% 40.02% 36.72% 


TasB_e XIII—Rapid City, South Dakota 


Oct. Nov Dec Jan Feb. March April May June July Aug. Sept. 


Insolation in therms ft*® (tilt correct 0.601 0.491 0.429 0.449 0.484 0.625 0.565 0.494 0.449 0.518 0.576 0.613 


inel 
House heat (therms 102.50 182.66 249.69 279.01 235.96 214.23 126.08 73.19 30.34 6.56 4.92 39.57 
Water heat (therms 26.66 27.30 29.76 31.31 28.28 29.76 27.30 26.66 24.30 23.56 23.56 24.30 
Degree days 500 891 1218 1361 1151 1045 615 357 148 32 24 193 
Aver. Ambient Temp. (roof 51.95 38.48 29.63 24.46 27.62 34.69 47.50 57.40° 67.09° 76.71° 74.62° 64.64° 
T, for V = 3000 gal; K/d = 0; A = 207 207 143.44 90.19 85.01 146 .57° 189.62 207° 207° 207° 207° 207° 


63.44% 56.26% 41.81% 35.29% 35.40% 38.38% 41.20% 40.60% 37.64% 


I ficiency l 
T,, for V = 3000 gal; K.d 0.2; A = 207 207 116.55 96.91 85.77 141.76 173.47 207° 207° 207° 207° 207° 


17 
1070 ft? 


Efficiency » 44.97% 59.457 62.3877 56.847 45.29% 38.14% 35.40% 38.38% 41.20% 40.60% 37.64% 


New York City and Long Island 


Feb March April May June July Aug. Sept 


Jan. 


Insolation in therms ft? (tilt correct 0.512 0.415 0.395 0.395 0.462 0.502 0). 468 0.477 0.461 0.509 0.499 0.533 


inc! 


House heat (therms 67.25 146.87 244.25 266.85 243.98 204.44 120.24 37.93 2.69 0 0 7.53 

Water heat (therms 26.32 26.97 29.42 30.97 27.97 29.42 26.97 26.32 23.97 23.22 23.22 23.97 

Degree days 250 546 GOS 992 907 760 447 141 10 0 0 28 

Aver. ambient temp. (roof 59.73 18.60 37.17 34.77 34.46 42.48 52.74 64.18° 73 .55° 78.42° 76.80° 70.71 

T, for J 3000 gal; K/d = 0; A 207 207 172.33 86.53 91.03 147.00 174.01° 207° 207° 207° 207° 207° 
1030 ft? 

Efficiency 7 36. 12° 38.647 54.76% 65.875 57.15% 46.54% 39.58% 37.50% 40.29% 41.69% 41.23% 39. 

T» for V = 3000 gal; K/d = 0.2; A = 207 207 131.44 87.70 85.72 137.69 152.56 207° 207° 207° 207° 207 
1095 ft? 

Efficiency 7 36.14% 45.47%, 60.62% 66.41% 59.29% 51.38% 43.20% 37.50% 40.29% 41.69% 41.23%, 39.46% 
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‘a 
— 
XV 
Oct. | Nov. | Dec. | | 
518) 


TasBLeE XV—Colombia, Missouri 


Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept. 


Insolation in therms/ft? (tilt correct 0.570 0.503 0.360 0.376 6.408 0.535 0.419 0.472 0.472 0.537 0.538 0.614 
incl.) 
House heat (therms) 61.18 152.71 230.93 254.75 204.55 162.98 76.12 31.52 3.27 0 1.40 14.48 
Water heat (therms) 25.85 26.52 28.95 30.50 27.55 28.95 26.52 25.85 23.52 22.75 22.75 23.52 
Degree days 262 654 | 989 1091 876 698 326 135 14 0 6 62 
Aver. ambient temp. (roof) 61.19° 45.93° 35.77° 32.41° 36.55° 45.74° 57.05° 66.78° 75.78° 80.92° 79.89° 72.39° 
T» for V = 3000 gal; K/d = 0; A = 207° 207° 207° 85.05° 116.80° 170.37° 207° 207° 207° 207° 207° 207° 
1130 ft? 
ide Efficiency 7 36.60% 31.75% 49.27% 61.83% 50.26% 38.08% 35.28% 38.28%, 40.93°% 42.39% 42.10% 39.95% 
at T» for V = 3000 gal; K/d = 0.2; A = 207° 207° | 180.56° 85.20° 123.83° 157 .87° 207° 207° 207° 207° 207° 207° 
1280 ft? 
Efficiency 36.58% 36.39% 53.33% 60.80% 51.13% 40.27% 35.28% 38.28% 40.93% 42.39% 42.10% 39.95% 
TasLe XVI—Grand Junction, Colorado 
Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept 
re Insolation in therms/ft? (tilt correct 0.689 0.545 0.520 0.512 0.494 0.622 0.546 0.519 0.539 0.542 0.624 0.691 
incl. 
\T House heat (therms) 2.84 173.25 247.63 278.03 202.13 161.44 87.94 31.72 5.03 0 0 7.88 
/lie Water heat (therms) 89 26.55 28.99 30.54 27.58 28.99 26.55 25.89 23.55 22.79 22.79 23.55 
Degree days , 792 1132 1271 924 738 402 145 23 0 0 36 
159 Aver. ambient temp. (roof) 57.31° 40.91° 31.35° 26.40° 34.40° 44.14° 55. 13° 65. 46° 75. 43° 82.07° 79.23° 71.74° 
& T» for V = 3000 gal; K/d = 0; A = 207° 207° 154.54° 110.39° 85.89° 166.77° 207° 207° 207° 207° 207° 207° 
870 ft? 
7: Efficiency 35.39% 39.37% 52.29% 61.04% 54.95% 38.26% 34.68% 37.88% 40.83% 42.71% 41.92% 39.76% 
is T» for V = 3000 gal; K/d = 0.2; A = 207° 207° 127 .06° 115.13° 85.39° 158 .47° 188.97° 207° 207° 207 207° 207° 
960 ft? 
: Efficiency 7 35.37% 43.94% 55.71% 60.41% 56.25% 42.72% 37.77% 37.88% 40.83% 42.71% 41.92% 39.76% 


Tas_e XVII—Atlanta, Georgia 


Nov. Feb. March May June July Aug. Sept. 


Oct. Dec. Jan. 


Insolating in therms, ft? (tilt correct 0.635 0.563 0.543 0.539 0.484 0.562 0.515 0.443 0.393 0.411 0.493 0.530 

incl.) 
House heat (therms) 38.50 137.55 214.90 221.20 179.20 141.40 16.55 7.00 0 0 0 2.80 
Water heat (therms) 24.80 25.50 27.90 29.45 26.60 27.90 25.50 24.80 22.50 21.70 21.70 22.50 
Degree days 110 393 614 632 512 404 133 20 0 0 0 8 
Aver. ambient temp. (roof) 65. 98° 54.39° 47.99° 47 .03° 49.19° 55.25° 64.91° 72.76° 80.46 82.50 81.75 76.71 
T, for V = 3000 gal; K/d = 0; A = 207° 207° 151.73° 87.43 94.80 104.31 181.67° 207 207 207 207 207 

605 ft? 
Efficiency 7 38.04% 43.82% 60.54% 68.35% 66.54% 55.47% 41.89% 40.06 42.26% 42.83% 42.62% 41.20% 
T, for V = 3000 gal; K/d = 0.01; A = | 207° 207° 148.44° 86.35° 93.84° 103.07° 179.40° 207° 207 207 207° 207 

605 ft? 
Efficiency 7 38.04% 44.36°% 61.17% 68.67% 66.88% 55.99° 42.28% 40.06° 42.26% 42.83% 42.62% 41.20% 
T, for V = 3000 gal; K/d = 0.02; A = | 207° 207° 144. 98° 85.34° 92.72° 101.73 177.14° 207° 207 207 207 207 

605 ft? 
Efficiency 38.04% 44.86% 61.78% 68.95% 67.22% 56 .50°% 42.65% 40.06% 42.26% 42.83% 42.62% $1. 20°% 
T, for V = 3000 gal; K/d = 0.05; A = | 207° 207° 138.82° 86.44 93.21 101.57° 173.77 207 207 207 207 207 

615 ft? 
Efficiency 7 38.07% 45. 86% 62.55% 68.71% 67.14% 57.03% 43.19% 40.06% 42.26% 42.83% 42.62% 41.20% 
T, for V = 3000 gal; K/d = 0.1; A = 207° 207° 126.32° 85.11° 89. 98° 98.03 166.07° 207° 207° 207° 207 207 

620 ft? 
Efficiency » 38.06% 47.87% 64.50% 69.36% 68.11% 44.43% 40.06% 42.26% 42.83% 42.62% 41.20% 
T, for V = 3000 gal; K/d = 0.2; A = 207° 207° 108.13° 88.18° 86.38 95.76° 155.45° 207° 207° 207 207 207 

640 ft? 3 


Efficiency » 38.06% 50.74% 66.64% 69.44% 68.91% 60.62% 46.14% 40.06% 42.26% 12.83% 42.62% 41.20% 


- 

‘ee 
od 


Oct 


Insolation in therms/ft? (tilt correct 0.5 
incl 

House heat (therms 42.7 

Water heat (therms 25 

Degree days 154 
Aver. ambient temp. (roof 64 

T, for V = 3000 gal; = 0; A = 207 
935 ft 

Efficiency » 37 

T» for V = 3000 gal; Kid 0.2; A 207 


1045 ft? 


Efficiency 7 


Oct 


incl 


House heat (therms 36.16 
Water heat (therms 25.51 
Degree days 86 
Aver. ambient temp. (roof 67 
T;, for V = 3000 gal; K/d = 0; A = 207 
1040 
Efficiency 38 61! 
T,, for J 3000 gal: K/d = 0.2: A = 207 


1255 ft? 


Efficiency » 


Oct 


House heat (therms 61 


Water heat erms 25 
Degree days 218 
Aver. ambient temp. (roof 61.3 


K/d = 


3000 gal: 


Efficiency » 36 
for 3000 gal: K/d = 0.02; A = 207 
650 ft 
Efficiency 7 
7: for V = 3000 gal; K/d = 0.1; A = 207 
675 ft 
Efficiency 7 36.64% 
T, for V = 3000 gal; K/d = 0.2; A = 207 


710 ft2 


Efficiency 7 


T,, for V = 10,000 gal; K/d = 0; A = 207 
625 ft? 
Efficiency » 36.64° 
T» for V = 10,000 gal; K/d = 0.02; A = 207° 
640 ft? 
Efficiency 7 36.64% 


for V = 10,000 gal; K/d = 0.1; A = 207° 
690 ft? 


Efficiency 36.64 


T, for V = 10,000 gal; K/d = 0.2; A = 207° 
745 ft? 


Efficiency 4 


Insolation in therms ft tilt correct 0.618 


Insolation in therms ft? (tilt correct 0.741 


36.647 


39.10% 


TABLE XVIII 


Dec 


TABLE XIX 


Dec 


TaBLe XX 


Dec. 


0.623 


45% 


2.09° 


70% 
42° 


Nashville, Tennessee 


Fresno, California 


Albuquerque, New Mexico 


July 


0.497 


0 
22.32 
0 


83 .27° 


207° 


43.04% 


207° 


43.04% 


July 


0.638 


0 
22.41 

0 
87.14° 


207° 


44.09% 


207° 


44.09% 


July 


0.610 


83.17° 


207° 


43.01% 


207° 


43.01% 
207° 


43.01% 


207° 


43.01% 


207° 


43.01% 


207° 


43.01% 
207° 

43.01% 
207° 


43.01% 


Nov Jan. Feb. March April May June Aug. Sept. 
51 0.418 0.383 0.345 0.409 0.482 0.489 0.509 0.458 7 0.506 0.570 be 
130.88 201.46 216.19 176.73 138.38 51.68 11.95 0 0 6.11 + 
f 2 26.10 28.52 30.07 27.16 28.52 26.10 25.42 23.10 22.32 23.10 a 
171 725 778 636 498 186 43 0 0 22 sa9-9 
2 §2.12 44.36 42.45 44. 88° §2.71° 62.97° 71.62° 79.81° §2.15° 75.75° 
207 179.80° | 112.72 85.67° 179.40° | 207 207° 207° | 207° 207° 
5% | 38.42% | 51.57% | 64.92% | 55.90% | 38.70% | 37.13% | 39.72% | 42.08% | MMMM, | 42.73% | 40.93% ae 
207 143.71 111.46 85.80 169.02° 179.92 207° 207° 207° 207° 
37.63° $4.46 57.23% 65.10% 57.40% 44.93% 41.63% 39.72% 42.08% 12.73% 40.93% 
= Nov = Jan. Feb March April May June [| Aug. Sept. vi 
0.578 0.389 0.424 0.522 0.650 0.608 0.633 0.598 0.599 0.623 
145.07 243.89 264.49 168.20 127.83 60.97 18.08 0 0 0 s : 
k 26.19 28.61 30.16 27.24 28.61 26.19 25.51 23.19 22.41 23.19 : 
: 345 580 629 400 304 145 43 0 0 0 
57.13 18.61 47.55 53.01 58.83 65.55° 72.47° 80.47° 84.78° 78.24° VOL. 
207 207 85.08 136.71 207 207° 207° 207° | 207° 207° 3 
35.34% 52.74% 62.91% | 45.85% | 35.88% | 37.91% | 39.97% | 42.27% | MMM | 43.46% | 41.64% 1959 
207 207 85.10° 158.28 207° 207 207° 207° = 207° 207° 
38.63% 35.34% §2.74% 59.80% 42.33% 35. 88% 37.91% 39.97% 42.27% 43.46% 41.64% 
=. Nov Jan Feb March April May June Aug. Sept. 
0.640 0.636 0.630 0.710 0.671 0.628 0.607 0.646 0.719 
in 
5.) 177.82 253.74 273.78 201.53 166.25 81.57 19.76 0 0 0 2.82 
7 25.95 28.37 29.92 27.02 28.37 25.95 25.27 22.95 22.17 22.17 22.95 - 
2 630 899 970 714 589 289 70 0 0 0 10 ar 
2 16.91 39.45 36.43 42.77 49.51° | 59.64° | 69.23° | 79.46° | MMM | 80.50° | 74.31° Shs rene 
for 207 207 132.10 91.35 85.58° 163.40 186 .85° 207° 207° 207° 207° 
645 ft 
44.827 60.61% 66.39% 57.71% 43.78% 39.49% 39.02% 41.42% 42.27% 40.51% 
207 128.42 91.82° 85.51° | 161.74 183.55° | 207° 207° = 207° 207° By eee 
15.45% | 61.08% | 66.33% | 57.96% | 44.56% | 40.03% | 39.02% | 41.42% P| 12.27% | 40.51% Pe eran, © 
: 207 116.41 95.99 86.08° | 157.55° | 173.33° | 207° 207° 207° 207° eh eee 
47.35% | 62.20% | 65.68% | 58.50% | 46.90% | 41.74% | 39.02% | 41.42% | MMMM | 42.27% | 40.51% Bena res hy 
207 105.68 102. 40° 86.10 154.89° 163 .53° 207° 207° 207° 207° 
es mms, 49.08% 62.81% | 64.77% | 58.89% | 48.88% | 43.35% | 39.02% | 41.42% | ME | 42.27% | 40.51% fe Ai 
207 157 . 56 107 . 52° 85.22° | 125.82° | 170.73° | 207° 207° 207° 207° 
40.63% | 54.MMM | 64.20% | 63.21% | 52.30% | 42.17% | 39.02% | 41.42% | MEM | 42.27% | 40.51% eae 
207 85. 25° 126.07° | 167.93° | 207° 207° 207° 207° 
41.54% “EH 64.56% | 63.17% | 52.68% | 42.62% | 39.02% | 41.42% 42.27% | 40.51% or 
207 132. 96 .64° 85.36° | 124.57° | 156.79° | 203.87° | 207° 207° 207° Bere 
59.77% 65.66% 63.38% | 54.62% | 44.96% | 39.50% | 41.42% 42.27% | 40.51% 
192.40 115. 86° 93.14° 85.55° | 121.55° | 145.84° | 184.98° | 206.26° 207° 207° 
49.78% 62.66% | 66.13% | 63.78% | 56.74% | 49.71% | 42.74% | 41.50% | | 42.27% | 40.51% 
38 


a. 


TasBLeE XXI—Los Angeles, California 


Nov. Feb. March 


Oct. Aug Sept. 


Dec. Jan. April May June July 


Insolation in therms/ft? (tilt correct 0.593 | 0.534 0.515 0.523 0.578 0.670 0.543 0.559 0.480 0.570 0.592 0.606 


incl. ) | 
House heat (therms) 39.82 91.55 137.78 173.03 139.61 124.97 84.68 55.39 25.63 14.19 10.07 25.63 
Water heat (therms) 25.11 | 25.80 28.21 29.76 26.88 28.21 25.80 25.11 22.80 22.01 22.01 22.80 
Degree days 87 200 301 378 305 273 185 121 56 31 22 56 
Aver. ambient temp. (roof) 65.74° | 61.97° 58.66° 55.33° 56.17° 58.93° 61.27° 64.14° 66.95° 69.91° 70.69° 69. 10° 
T, for V = 3000 gal; K/d = 0; A = 207° | 207° 171.49° 137.40° 86.65° 156.04 174.37° 191.81° 207° 207° 207° 207 
475 ft? 
Efficiency 7 38.00% | 42.76% 52.84% 64.54% 62.07% 49.57% 44.53% 40.01% 38.34% 39.22% 39.45% 38.98% 
T, for V = 3000 gal; K/d = 0.2; A = 207° 199.72° 130. 40° 125.79° 85. 87° 151.57° 150.73 162.48° 196.43 201.47° 207 207° 
550 ft? 
Efficiency 7 39.19% 50.45% 60.75% 66.32% 62.83% 53.94% 52.88% 46.51% 41.00% 40.10% 39.45% 38.98% 


TaBLE XXII—Charleston, South Carolina 


Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept 


Insolation in therms/ft? (tilt correct 0.562 0.529 | 0.461 0.547 0.484 0.559 0.543 0.530 0.476 0.489 0.513 0.483 


incl.) 
House heat (therms) 18.20 | 94.50 159.60 165.20 132.65 98.35 22.05 0 0 0 0 0 
Water heat (therms) 24.92 25.62 28.02 29.57 26.71 28.02 25.62 24.92 22.62 21.82 21.82 22.62 
Degree days 52 270 456 472 379 281 63 0 0 0 0 0 
Aver. ambient temp. (roof) 69.56° 59.54° 54.04° 53.57° 55.02° 60.11° 68.14° 74.22° 80. 48° 82.83° 82.63° 78.73' 
T; for V = 3000 gal; K/d = 0; A = 207° 207° 173.12° 85.89° 136.30° 126.91 196. 85° 207° 207° 207° 207° 207° 

Toe 505 ft? 

Efficiency 7 39.14% 41.79% 59.18% 64.37% 58.77% 50.97% 40.35% 40.48% 42.27% 42.92% 42.87% 41.78% 
T, for V = 3000 gal; K/d = 0.2; A = 207° | 207° 127.80° 85.15° 131.10° 111.37° 172.72° 207° 207 207° 207° 207 

575 ft? 
Efficiency 39.15% 49.06% 65. 80% 65.23% 61.84% 56. 98°7 44.287, 40.48% 42.27% 42.92% 42.87% 41.78% 


TaBLE XXIII—ElI Paso, Texas 


Nov. Dec. Jan. Feb. March April May June July Aug. Sept. 


Insolation in therms/ft? (tilt correct 0.725 0.670 0.649 0.661 0.654 0.745 0.679 0.648 0.587 0.592 0.658 0.706 


incl.) 
House heat (therms) 23.61 131.53 211.12 225.96 150.08 111.29 37.10 0 0 0 0 0 
Water heat (therms) 24.77 25.47 27.88 29.42 26.57 27.88 25.47 24.77 22.47 21.67 21.67 22.47 
Degree days 70 390 626 670 445 330 110 0 0 0 0 0 
Aver. ambient temp. (roof) 68. 08° 55.30° 48. 88° 46.91 52.82° 58.10 66. 88° 75.35° 84.13° 84.93 83.25 78.65° 
T, for V = 3000 gal; K/d = 0; A = 207° 207° 155.53° 86 .56° 89.07° 179.32° 207 207° 207 207 207° 207° 
495 ft? 
Efficiency 7 38.70% 43.46% 60.34% 69.25% 57.44% 40.29% 38.31% 40.81% 43.28% 3.50% 43.04% 41.75% 
T» for V = 3000 gal; K/d = 0.2; A = 207° 207° 114.07° 90.00° 86.51° 158.92 174.12 207° 207 207 207 207 
535 ft? 
Efficiency 38.69% 50.07% 65.77% 69.12% 60.85% 48.97% 43.70% 40.81% 43. 28% 43.50% 13.04% 41.75% 


TaBLE XXIV—Fort Worth, Texas 


Oct. Nov. Dec. Jan. Feb. March April May June July Aug. Sept 


Insolation in therms/ft? (tilt correct 0.674 0.577 0.542 0.483 0.541 0.617 0.573 0.549 0.555 0.581 0.623 0.669 
incl. ) 
House heat (therms) 16.92 87.23 155.50 181.47 130.12 89.86 26.26 1.46 0 0 0 0 
Water heat (therms) 24.92 25.62 28.02 29.57 26.71 28.02 25.62 24.92 22.62 21.82 21.82 22.62 
Degree days 58 299 533 622 446 308 90 5 0 0 0 0 
Aver. ambient temp. (roof) 71.41° 58.32° 51.49° 48.09 52.47° 58 .52° 68.42 75. 85° 84.34 88.28 88.67 81.47 
T;, for V = 3000 gal; K/d = 0; A = 207° 207° 207° 134. 87° 86.12° 188.94 207° 207° 207 207 207° 207 
555 ft? 
Efficiency 7 39.66% 35.72% 45.69% 63.14% 56.35% 38.86% 38.78% 40.95% 43.34% 44.40% 44.50% 42.54% 
T, for V = 3000 gal; K/d = 0.2; A = 207° 207° 162.35° 123.43 85.76° 175.35° 181.72° 207° 207 207° 207° 207 
635 ft? 
Efficiency 7 39.68% 43.15% 54.50% 64.86% 58.42% 45.25% 42.90% 40.95% 43.34% 44.40% 44.50% 42.54% 
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Oct Nov Dec Jan. 
Insolation in therms ft? (tilt correct 0.737 0.632 0.594 0.616 
incl 
House heat (therms 8.77 88.92 159.50 189.01 
Water heat (therms 25.02 25.71 28.12 29.67 
Degree days 22 223 400 474 
Aver. ambient temp. (roof 73.92 62.18 56.30° 52.40° 
T, for V = 3000 gal; K/d = 0; A = 207 207 191.86 112.78° 
445 ft® 
Efficiency 7 40.40% 39.41% 52.85% 66.85% 
T, for V = 3000 gal; K/d = 0.2; A = 207° 207 136.98 100.04° 
485 ft* 
Efficiency 40. 40° 48.34% 62.92% 69.28% 


auxiliary heat at all. Six different schedules were 
assumed and calculated. The schedules are given in 
Table XXIX and the results shown in Fig. 9. The 
third schedule is based on the figures given by 
Whillier? for the M.I.T. house, which—as explained in 
Appendix A—was built ‘tight’ so as to reduce dras- 
tically the convection losses generally found in a resi- 
dential house. If one regards the decrease in heat 
demand due to the tight construction as though it 
were auxiliary heat furnished to our standard house of 
average tightness, the third schedule (B;) is the result. 

The fourth, fifth, and sixth schedules are chosen to 
ascertain the effect of different amounts of auxiliary 
heat on the required collector area together with the 


other parameters. 


FURTHER ECONOMIC CONSIDERATIONS 
A. Cost of storage 


The cost of energy storage is conveniently separated 
into cost of the tank and cost of insulation. Steel tanks 
are more practical than wooden tanks, and, in the 
larger sizes, cheaper as well. (See Appendix B.) A 
3,000-gal tank of diameter 5 ft 4 in. and 18 ft long can 
be built of ;°;-in. steel, and costs $350. Larger diameter 
tanks must be built of }-in. steel. An 8 ft x 8 ft tank 
(3,000 gal) costs $500; a 10 ft x 103 ft tank (6,000 gal) 
costs $750; a 103 ft x 153 ft tank (10,000 gal) costs $982; 
a 12 ft x 12 ft tank (10,000 gal) costs $1,300. These 
prices are not for quantity orders, but the saving can 
be considered canceled out by the additional cost, not 
figured in the above, of rust prevention treatment. 
(Data from the Buffalo Tank Corp.) 

In the 3,000-gal tank size, the longer tank is more 
economical than the square tank. The shape factors 
are 6.36 and 5.54 respectively—a 13 per cent dis- 
advantage for the long tank, but 13 per cent increase 
in insulation cost is significantly less than the $150 
difference in tank price. The same conclusion holds 
a fortiori for the case of the two 10,000-gal tanks; the 
long tank has a shape factor of 5.625, an increase of 
only 13 per cent, while the insulation cost which is to 


TaBLE XXV—Phoenix, 


F 


0 


123 


26. 
309 

57. 

89. 


55. 
85. 


59. 


Arizona 


eb March April May June July Aug. Sept. 

674 0.727 0.680 0.668 0.612 0.596 0.628 0.717 

21 78.16 29.51 0 0 0 0 0 

80 28.12 25.71 25.02 22.71 21.92 21.92 22.71 
196 74 0 0 0 0 0 

78° 63 .56° 71.15° 79.43° 88.37° 93. 58° 91.60° 87.14° 

80° 203 .24° 207° 207° 207° 207° 207° 207° 

01% 37.97% 39.59% 41.97% 44.42% 45.79% 45.28% 44.09% 

22° 177.20° 179.40° 207° 207° 207° 207° 207° 


60% 46.70% 44.06°% 41.97% 44.42% 45.79% 45.287, 44.09% 


be increased by this amount is not quite double that 
for the 3,000-gal tank. 

A plot of the different tank capacities against cost 
does not result in a smooth curve; moreover, for lowest 
cost, tanks must be purchased in stock sizes. Thus, it is 
better to consider tank costs as quantized, rather than 
as so many dollars per gallon capacity. 

The cost of insulation is figured on the basis of data 
furnished by the Pittsburgh Plate Glass Co. on their 
Superfine Fiber Glass of 1 lb per cu ft density and K 
factor at 137°F of 0.27 Btu per sq ft per hr per °F 
for 1-in. thickness. It comes in batts either 3 in. thick 
or 5 in. thick and the price for quantity orders for either 
thickness figures out to 6.7 cents per sq ft per in. thick- 
ness. This means, for example, that the 18 ft long, 3,000- 
gal tank mentioned above can be swathed in 15 in. of 
insulation for about $400, 15 per cent allowance being 
made for the fact that the outer batts must cover a 
larger area than the innermost batt. 

In determining the most economical combination of 
collector area, tank capacity, and insulation, it is con- 
venient to express the insulation cost as 1.8 cents per 
sq ft for an insulating value of 1 Btu per sq ft per hour 
per °F. The change in A factor with temperature is 
neglected. 


B. Optimizing the Parameters 

On a strictly economic basis it is cheaper in most 
localities to save some $400 invested in insulating the 
house when it is constructed and pay perhaps $35 more 
per year for fuel oil. (See Fig. 7). However, such a house 
will be less comfortable in extremely cold or extremely 
hot weather, which is the reason people are willing to 
pay for the insulation. Similarly, instead of making a 
“tight”? house at an extra construction cost of $500 or 
$1,000, a solar heating system, since it requires an extra 
capital outlay in the best of cases, should aim at ways 
to reduce the construction costs of the house. 

In this report we are not concerned with optimizing 
the proportions spent on different items in house con- 
struction but only with optimizing the items, or param- 
eters, used in the present calculations, i.e., the cost 


= 
; VOL. 
1959 
: 
= 
2 
40 


of collectors and of storage. The calculations are made 
for the standard house described above. The measure 
of economic feasibility is again taken as the amount 
which can be afforded for each square foot of collector, 
assuming a 10-year amortization period, and Equations 
[2] and [3] are used to calculate 0M/dA. Equation [4], 
however, is no longer adequate; instead we use: 
1 


300 
ad + Vay + \600 


AM =A 

The third term on the right is the cost of insulation; 
the partial derivative is the cost of insulation per unit 
area for a thickness and material of 1 Btu per sq ft per 
hr per °F’. The value of K/d is thus one of the variables, 
or rather f is, since A is fixed at 0.27 Btu per sq ft per 
hr per °F in. thickness. 

The problem is to find the set of values for A, V, 
OM /dV,S, and d which will minimize AM, or rather 
which will maximize 0//dA. We thus need not only 
Equation [10] but also the relation 


12 
AM [il] 
1 


5 
where the factor of 0.18 is our assumed fossil fuel rate 
of 18 cents a therm, and B is the B, Schedule of Table 
XXIX. The first step is to assume a value for V of 
either 3,000, 6,000 or 10,000 gal. Then, from the storage 
tank cost and dimension information given in Part A 
of this Section, we get 0.1//0V and S. We next pick a 
value of A/d and interpolate the corresponding value 
of A from the calculated data, and, substituting into 
Equations [10] and [11], compute 0.17 /0A. The process 
is then repeated for different pairs of A /d and A values 
until a maximum value of 0.1/7 /0A is obtained. Then a 
different choice is made for V, and the whole process 
again repeated. 

This trial and error procedure was used for only three 
localities: Lander, Blue Hill, and Albuquerque. It then 
became obvious that the optimum storage tank size 
for all localities is 3,000 gal. We cannot consider tanks 
much smaller than 3,000 gal or the storage will not be 
able to supply house heat during several days of suc- 
cessive cloudy weather, i.e., there is a minimum even 
for short-term storage. 

In Albuquerque, 0.17 /0A is $2.52 for a 3,000-gal 
tank, and about 70 cents for a 10,000-gal tank. As can 
be seen from Fig. 11, the 10,000-gal tank must be well 
insulated (long-term storage) in order to reduce the 
required collector area. In Lander, which can be seen 
from Fig. 10 to be a much more favorable place for 
large storage tanks, 0.7 /0Amax is $1.06 for a 3,000-gal 
tank, and 86 cents for a 10,000-gal tank. Lander is 
probably the best place of the twenty-one for large- 
‘apacity storage; when the smallest tank turned out to 
be the cheapest there too, further calculations were 


limited to houses with 3,000-gal storage tanks. Accord- 


ingly, V €M/dV was taken as $350, and S as 400 sq ft. 


Thus, Equations [10] and [11] reduce to 


am if. 400 X 0.018 . {450 
aA A (s [12] 


Figs. 9-11 led to a further simplification in calculat- 
ing the optimum values of the parameters. It became 
clear that a sufficiently good approximation could be 
obtained by assuming a linear relationship between 


A and K/d. 


_ AA 


A= 
0.2 


K/d + Ao [13] 
where Ao = required collector area when K/d = 0 and 
AA = difference in required collector area for K/d = 0 
and K/d = 0.2. 

The optimum value of Kd can then be found directly 
by setting equal to zero the derivative of 0M/dA with 
respect to A/d. Thus we get 


44 + 4/ 2.0736 + 0.288 (a1, Ay/AA 


d opt. 50 


90 


The results of these calculations are tabulated in 
Table XXVI. Fig. 1 is drawn on the basis of the values 
(maxima) found for @M/dA. The insulation costs, 
S0M/dS are found simply by dividing the optimum 


Albuquerque 
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CENTS / THERM OF FOSSIL FUEL 


Fig. 12—Maximum allowable cost for each square foot of solar 
collector vs. cost of heating fuel. 
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TaBLeE XXVI 


V = 3000 gal 


Location 


Blue Hill, Mass 
Central Mass 
Ithaca, N. Y 
Madison, Wis 
Rhode Island 
Columbus, Ohio 
Medford, Ore 
Lander, Wyo 
Rapid City, 8. D 
New y ork 
Long Island 
Columbia, Mo 
(Grand Jetn Colo 
Nashville, Tenn 
Atlanta, Ga 
Fresno, Calif 
Albuquerque N M 
Los Angeles, Calif 
Charleston, 8. C 
kK] Paso, Texas 
Ft. Worth, Texas 
Phoenix, Ariz 


a 


10,000 gal 


sander, Wyoming 


TaBLeE XXVII 


Blue Hill, Mass 
Central Mass 
Ithaca, N. ¥ 
Madison, Wis 
Rhode Island 


Columbus, Ohio 
Medford, Ore 
Lander, Wyo 


Rapid City, 8. D 
N. Y. City and Long 
Island 
Columbia, Mo 
Grand Jetn Colo 
Nashville, Tenn 
Atlanta, Ga 
Fresno, Calif. 
Albuquerque, N. M 
Los Angeles, Calif 
Charleston, 8. C 
Kk] Paso, Texas 
Ft. Worth, Texas 
Phoenix, Ariz. 


City 


6392 
6743 
6914 


Optimum Parameters 


Collector gy 


Insula Insula 
tio 
Btu f 
SSS 
OS S7 
12 62 
19 39 
7a 
10 79 
67 
06 120 
SI 
1] 65 
OS SO 
O9 79 
09 76 
14 53 
OS 4D 
10 76 
OY S|] 
10 69 
12 58 
10 70 
12 61 
aM 
—| 


Fuel 


avings 


$362 
$361 
$330 
$546 
$349 
$309 
$344 
$372 


$336 


$299 
$27 

$285 
$231 
$232 
$247 
$278 
$221 
$179 
$214 
$178 
$176 


785 
760 
725 
715 
695 
640 
615 
575 
570 


565 
$95 
440 
$50 
400 


area A 


aA max 


it 

1401 $1.50 
1487 1.41 
1641 1.20 
1598 1.30 
1256 1.63 
1438 1.29 
1909 0.96 
2010 1.06 
1017 1.95 
1063 1.71 
1192 1.38 
910 1.90 
QS5 1.49 
626 2.07 
1121 1.18 
674 2.33 
505 2.39 
540 1.86 
517 2.32 
592 1.70 
166 2.16 
= $14.15 


1640 


Simplified Calculation 


$672 
S662 
$652 
$660 
$632 
$644 
$712 
$596 
$528 


$570 
$536 
$536 
S488 
S484 
$532 


$332 
$420 
$348 
$316 


te 
! 

' 


K d values into $7.20, i.e., 400 sq ft times 0.18 dollars 


per sq ft for 1 Btu per sq ft per hr per °F. 


C. Changes in Fuel Cost 

Waiting for conventional fuel to become expensive 
(while presumably the cost of building and installing 
solar equipment remains unchanged) is a sad method of 
“improving” solar house heating. However, it is worth- 
while investigating the effect of such an increase in 
fuel costs. Fig. 12 shows the permissible cost per square 
foot of collector plotted against the cost per therm for 
fossil fuel used in three different cities. Where the graph 
crosses 18 cents per therm corresponds to the data 
given in Table XXVI. 

The calculation is made the same way as that de- 
scribed in Part B of this section, but with the fuel sav- 
ings proportional to the value assumed for the cost of 
fossil fuel. It is not hard to show that the curves should 
be approximately linear. They can be represented with 
fair accuracy by the equation: 

aM. 5M; 450K /d — 7.2 

“A, * A)K/d 
where p is the ratio of the cents per therm paid for fossil 
fuel to the assumed value of 18 cents per therm. 


[15] 


DISCUSSION OF RESULTS 
QOne’s first instinct when confronted with the mass of 
data presented in the tables and figures is self-preserva- 
tion; i.e., reduction of the data to a few simple general 
equations. This can be done only on a limited basis, but 
it is not worthwhile, for example, to try to fit the annual 
insolation or heat load data with sinusoidal expressions. 
It can be seen from Figs. 9-11 that, for a given stor- 
age capacity, there is a nearly linear relationship be- 
tween area (A) and insulation (AK/d), as 
expressed in Equation [13]. But if we try to extend the 
idea to include storage capacity and find a relation such 
that the constants hold for all localities, we can expect 
disappointment. For example, a relation of the form 


collector 


4 f 


[16] 


i 
1 


where the C’s are constants, was investigated and found 


A=(C,V°K/d+C; 


not worthwhile. Instead it seems better to be content 
with observing: 

(1) The slope of A vs. A/d curves increases with V. 

(2) Ay with 10,000-gal storage is always less than 
Ay with 3,000-gal storage. 

(3) A increases with the severity of the winter (degree 
days) and decreases with the amount of winter insola- 
tion. Summer insolation makes no difference because 
heat is thrown away in the summer anyway unless the 
storage capacity is very large and the storage efficiency 


Fuel 
$349 
348 
316 
332 
336 
206 
331 
359 = 
323 
286 
25S 
271 
218 
219 
934 
265 
i 20S 4 
165 
202 
166 
164 VOL. 
= 
a as 
359 07 202 0.86 
— 
No 
auxil Col. Cost of 
days ‘te ous fe (a $1 $ ft 
therm 
H + ee 
1.06 
6125 1.059 
561: Or 
1547 1.40 
7535 | .82| ‘ 
1989 1.076 $3.05 
5113. $3.27 
5796 875 $3.84 
3513. 1.112 $3.41 
2826 1.40 $4.31 
2532 1.682 380 $3.43 
4389 1.129 375 $532 $4.42 
2015 1.831 340 $268 $4.22 
1973 | 1.40 300 $3.87 
2641 1.349 280 $4.47 
2361 1.167 265 $5.07 
1698 1.595 230 $5.07 
42 


very high. (See section on “Simplified Calculation,” 
above.) 

Inspection reveals the great sensitivity of the calcu- 
lation to collector efficiency (Fig. 3). For example, 
in Table V the temperature in the 3,000-gal storage 
tank on January Ist is 85.03° when the insulation leaks 
0.1 Btu per sq ft per hour per °F, but 85.13° when the 
insulation leaks heat twice as fast. In both cases, the 
collector area is taken as 1410 sq ft, although it is clear 
that a very slightly smaller collector area, 1408 or 1409 
sq ft, would suffice in the second case. The chief reason 
for this paradox—that lower storage efficiency makes 
the over-all! system more efficient—is indicated in the 
lines immediately below the temperatures, i.e., the 
collector efiiciencies. When the heat leaks out faster, the 
collector temperature is lower; in November and De- 
cember, which are the critical months, the gain of 3 
per cent in collector efficiency more than outweighs the 
greater heat wastage of the less thickly imsulated 
storage tank. 

Another striking example is afforded by Figure 9. As 
shown by Table XXIX, the B, schedule supplies more 
auxiliary heat than the B, schedule in November, 
February and March, and the same auxiliary heat in 
December and January. Yet the B, schedule requires a 
larger collector area, for all storage efficiencies, than 
does the B, schedule. As with the first example, the 
computer calculation (not given for By) was checked 
by desk calculator and the paradox traced to the higher 
storage tank temperature (due to the auxiliary heat 
supplied in November) causing a lower collection effi- 
ciency. These two examples illustrate the value of 
‘arrying out the calculations to four significant figures, 
even though agreement with measurements would not be 
expected to be better than the second significant figure. 

It is clear that auxiliary heat should be furnished with 
‘aution, since it may not give real assistance to the solar 
heating system. To calculate the optimum auxiliary 
heat schedule would be very tedious, but the superiority 
of B, over B, suggests that it should be small. Auxiliary 
heat furnished before the minimum tank temperature 
is reached, i.e., January Ist, lowers the collection effi- 
ciency. Auxiliary heat furnished after January Ist is 
generally unnecessary. However, if no auxiliary heat at 
all is furnished, as in schedule Bo, the result—as can be 
seen from Fig. 9—is a substantial increase in the re- 
quired collector area. It is probable that B, is quite 
close to the optimum. (As pointed out in Principal 
Conclusion No. 13, this does not mean the most eco- 
nomical. ) 

The extreme sensitivity of the system to collection 
efficiency means that the maximum temperature at 
which the water is stored should, in an optimized 
system, be carefully determined. The paradoxes indi- 


cated above strongly imply that 207°F is too high a 


temperature. The optimum maximum storage tempera- 
ture should probably be about 190°F. This would reduce 
the energy storage capacity of the tank, but the gain 
in collection efficiency would (especially in November) 
reduce the required collector area and thus the cost. 

The optimum maximum temperature was not calcu- 
lated, however, because the result would be consider- 
ably in error as long as the calculation is made on a 
month-to-month basis, i.e., the collection efficiency 
does not vary with the day-to-day temperature in the 
storage tank, but is taken as constant over an entire 
month. It is even possible that neither of the above 
paradoxes would appear, i.e., actual cross-overs from 
the expected result, if collection efficiency were not 
quantized monthly. However, the conclusions drawn 
would probably remain valid. 

No calculations were run for collectors with effi- 
ciencies different from those shown in Fig. 3. To a first 
approximation, however, the required collector areas 
can be taken as inversely proportional to the collection 
efficiencies under certain more-or-less average con- 
ditions: 


(17] 
2 Aly 


Equation [17], it will be recalled, holds accurately for 
the simplified calculation. (See Fig. 6.) 

One other paradox deserves mention. In Table V, 
line 1, it can be noted that the insolation for December 
exceeds that for November, although the December sun 
is lower in the sky and the days are shorter. The chief 
reason is indicated in Table IV, i.e., the horizontal-to- 
tilt ratio. A secondary reason is that December has one 
more day than November. 

We turn now to the results of the calculations of 
economic feasibility. The significance of the values com- 
puted for 01//dA depends on the actual present-day 
costs of building and installing collectors. Reasonable 
estimates seem to be: custom-made collectors would 
cost more than $10 per sq ft; if collectors were mass- 
produced, the cost could probably be brought below 
$5, perhaps below $3. In the next decade fuel costs 
will probably rise, but so will materials and labor, so 
that little change in feasibility is to be expected from 
this direction. A breakthrough, on the other hand, that 
would permit collectors to be built and installed for 
$1.00 or $1.50 per sq ft, would make a decisive change 
in feasibility prospects. 

This means that a breakthrough in cost of collectors, 
to make solar house-heating competitive, must accom- 
plish a price improvement by a factor of two or three. 
A breakthrough in storage costs to accomplish the same 
thing would have to be more extreme. In the simplified 
‘alculation above, a factor of 10 was assumed, i.e., 90 
per cent reduction of cost. The fact, mentioned above, 
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that the ‘elbow’? values in Fig. 5 are the most eco- 
nomical on each curve means that 90 per cent reduction 
is somewhat greater than that necessary. In other 
words, a factor of 7 or 8 in storage cost reduction might 
suffice. This assumes, of course, that the collector costs 
indicated in Fig. 2 are realistic. The minimum factor 
which would do the trick depends, however, not only on 
collector costs (0./'0A) but also on the collector 
efficiency, which was assumed in the simplified calcula- 
tion to be a constant 50 per cent. 

A breakthrough in storage costs would almost auto- 
matically make long-term storage more economical 
than short-term storage. As indicated in Fig. 5, large 
reductions in collector area (A) are achieved only by 
using very large storage capacities; such large capacities 
are advantageous only if summer heat is saved for 
winter use, i.e., long-term, high efficiency storage. Thus, 
if storage reduces collector area to where present col- 
lector costs can be met, the storage would be both of 
high efficiency and large capacity. 

The similarity between Fig. 1 and Fig. 2, e.g., the 
relatively favorable areas being in the Southwest and 
the Appalachians, show that the areas suitable for 
solar house-heating are nearly independent of the cost 
of energy storage. In other words, whether or not there 
is a breakthrough in cost of energy storage, the same 
sections of the country remain the most favorable sec- 
tions, and the same sections of the country remain the 
most unfavorable. 

It should be noted that the accuracy of Figs. 1 and 
2 is quite poor; not only were they drawn chiefly on 
the basis of 21 points, but the significant variations in 
climate known to occur within short distances (e.g., 
between Blue Hill and Boston) were ignored. In several 
cases, such as Tables VI and XIV, degree-day data for 
a locality was combined with isolation data for an ad- 
jacent area to complete the set of input data. 

Both Figs. 1 and 2 show the southern tip of Florida 
as an uneconomic place for solar house heating. The 
chief reason is that the annual fuel savings (.V/;) are 
too small. (See Fig. 7.) Since we have stipulated that 
the furnace is to be eliminated, some importance at- 
taches to what might be called the ‘‘furnace line,” 1.e., 
the geographical line south of which the better class 
one-family houses do not have central heating. Inquiry 
as to where this line should be drawn indicates that it 
is moving southward; that is, in new houses along the 
Gulf Coast there is a growing tendency to include cen- 
tral heating. Even in Miami, the electric company is 


yopagandizing for central heating because the elec- 
g 


tric load on cool evenings, when everyone turns on their 
electric room heaters, constitutes a utility problem. The 
line passing through places having an average January 
containing between 300 and 330 degree days hits San 
Diego, Yuma, Corpus Christi, and Jacksonville. New 


a 


Fic. 13—Cost for each square foot of solar collector vs. 
thermal conductivity of storage tank insulation for 
3,000-gal tank costing $350 and insulation costing $7.20/ 
{ 

K/d). 


Orleans, Mobile, and Thomasville (Georgia) have 
slightly more. 

In this connection it may be noted that the widely 
reproduced feasibility map by Siple (Reference 6, page 
15), considers only ‘“‘mechanical feasibility,” 1.e., ‘‘eco- 
nomic considerations were not included.” Some writers 
on solar energy have not been too careful about ob- 
serving the distinction. 

It was shown above that the dividing line between 
long-term and short-term storage is about 0.02 Btu 
per sq ft per hr per °F. Thus, it can be seen from the 
column of insulation factors in Table XX VI that short- 
term storage is cheaper than long-term storage. This 
means that the idea of storing summer heat for winter 
use does not pay at the present cost of energy storage. 

It can be noted that the variation in the economically 
optimized insulation factor for different parts of the 
United States is quite small, averaging about one-fifth 
the insulation required for long-term storage. However, 
the thickness of insulation is not critical as far as over- 
all expense is concerned; using Equation [10], 0.17/04 
goes through a flat maximum when graphed against 
K d. The curve falls more steeply on the side of the 
thicker insulation because the cost of insulation be- 
comes infinite as Ad approaches zero. Two examples, 
for Albuquerque and Lander, are shown in Fig. 13. 

It is difficult to draw general conclusions from Table 
XXVI. For example, there is no correlation bet ween the 
insulation costs (S 0.7 0S) and the optimized collector 
area or the optimized collector cost. There is a slight 
correlation between the latter two; that is, a small col- 
lector area is usually associated with a high permissible 
collector cost, but there are exceptions even to this 
qualitative formulation. Compare, for example, Blue 
Hill and Rhode Island with Nashville, Columbia, and 
l’resno. 

If the storage tank is placed inside the house instead 
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of in the basement, the ambient temperature is higher 
and the heat leakage rate per inch thickness of insula- 
tion is less. Moreover, the heat that leaks out in winter 
is not lost but helps to heat the house, thus providing 
almost 100 per cent efficiency. L6f* and Telkes* have 
ach built a solar house embodying this idea. If there 
is a full basement, extra lolly columns may be needed 
to support the tank. The tank itself could be part of 
the structure of the house, preferably at an inside cor- 
ner, so that only a half or a quarter of the cylindrical 
surface protrudes into a room or closet. However, only 
about 15 per cent saving on insulation is achieved be- 
cause the heat must still be stored for the same length 
of time. If the house is designed for short-term storage, 
the percentage saving in insulation is greater, but there 
is not much insulation cost to begin with. Moreover, 
the heat leakage in summer will probably be objection- 
able. In such cases, cold can be stored in the tank. (See 
the following section.) 

No consideration was given in the calculations to the 
collection of solar energy through south-facing windows. 
As indicated by Anderson® the losses for a sunny winter 
day (in New England) calculated for a 24-hr period, 
will exceed the solar radiation transmitted into the 
house unless the window is double-paned and/or cur- 
tained at night. 


AIR CONDITIONING 


The most expensive part of the solar equipment is 
the collector, and it would be advantageous if the same 
collector could be used to operate a cooling system in 
summer as well as the heating system in winter. Since 
there is considerably more solar energy available to 
the collector in the summer than the heating system can 
make use of, no additional collector area is required if 
a solar-powered air-conditioning unit is added. More- 
over, some parts of the country, particularly the South- 
west, where air conditioning is needed, are also the 
most favored for solar house heating. (These areas 
sometimes have such low humidity that the air can be 
satisfactorily cooled by simple water-spray devices 
rather than refrigeration machinery. ) 

However, in most parts of the United States, the 
annual cost of air conditioning is too small for a solar 
system to pay for itself as a fuel (i.e., electricity) saver. 
Almost any one-family house of 2,000 sq ft or less can 
be cooled with a 2-ton or 3-ton unit. In southernmost 
United States, the cost of electricity for a season of air 
conditioning is generally less than $150, and in most 
parts of the country it is less than $75.’7 '* This is de- 
spite the fact that the cost of removing a Btu may be 
more than double the cost of furnishing a Btu. Fig. 7 
shows that such a small fuel saving can be justified 
only if the cost of the solar equipment does not exceed 


that of the conventional equipment by more than $200 
or $300. 

Secondly, flat-plate collectors which are quite sat- 
isfactory for heating houses do not, under usual weather 
conditions, yield high enough temperatures to be suit- 
able for a refrigeration system. Thus, if collectors are 
to be used for both heating and cooling, the more ex- 
pensive kind needed for cooling must be installed (at 
least, some of the collector area must have the higher 
temperature type.) 

Higher temperatures are needed whether one uses a 
solar-powered compressor with a mechanical air con- 
ditioning system or an absorption-type refrigeration 
cycle. Tabor, who made a theoretical comparison be- 
tween these two alternative methods, has expressed the 
opinion that there is not much to choose between them. 
From the standpoint of considering economic feasibility, 
it is easier to consider the solar compressor. Although in 
practice there is a complicated engineering problem in- 
volved, in principle the solar compressor consists chiefly 
in replacing an electric motor with a cyclical heat en- 
gine, the rest of the air conditioning system being simi- 
lar in the conventional and the solar system. Between 
the difficulties of building a cheap substitute for a 
simple electric motor and the high temperature re- 
quirement on the solar collector, solar air conditioning 
does not appear very promising. 

However, solar house heating can make possible a 
saving if the house is air-conditioned. If the storage 
tank is used to store cold water, a conventional refrig- 
eration unit can be run more or less continually during 
the hot season cooling the water. The required capacity 
of the refrigeration unit is thereby reduced from per- 
haps 3 tons to 1 ton or even } ton. The kilowatt hours 
of electricity consumed are about the same, but the 
difference in initial cost between the large and small 
unit represents a one-shot saving made possible by the 
solar storage tank. The small second tank can be used 
to provide house heat on occasional cool evenings. How- 
ever, this may not be satisfactory in many regions where 
the seasons change irregularly and ‘“‘unseasonable”’ tem- 
peratures are common. 

The question of economic feasibility for solar cooling 
involves two questions: (1) Does the additional capital 
investment involved in adding solar cooling to a solar- 
heated house justify the additional savings achieved? 
As indicated, this question can be answered from I ig. 
7. (2) Is combined solar cooling and heating more or 
less nearly economically feasible than solar house heat- 
ing alone? Even if solar heating and cooling are econom- 
ically infeasible when considered separately, is the com- 
bination any better? 


for each case, certain costs 


aM 
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for the equipment must be assumed. This can be 
avoided, however, by turning the problem around and 
calculating instead the maximum extra capital invest- 
ment (that is, the cost of the solar cooling equipment 
minus the cost of installing equivalent conventional 
air conditioning) that the solar house owner can pay 
without increasing the collector cost per square foot. 
It can be shown from Equations [2] and [10] that 
OM @A remains unchanged if the extra capital invest- 
ment required to add solar air conditioning to a solar- 
heated house equals five times (assuming 10-year amor- 
tization period) the annual savings accomplished by 
the solar air conditioning unit. For example, an annual 
electric bill saving of $50 justifies an extra expenditure 
up to $250; if the solar air conditioner can be made and 
installed for less than $250 extra, then combined solar 
heating and cooling is more nearly economically feasible 
than solar heating alone. Whether or not this figure is 
attainable at today’s prices was not investigated. 
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TaBLe XXVIII 


o = 5.54 V = 3000 gal K/d V = 6000 gal V = 10,000 gal 
L (per min) (Btu/ft?/hr/°F) Kid K/d 


1x 0.499 KX 10°? | 0.628 K 10-? | 0.745 XK 107? 
2 0.997 1.257 1.490 
3 1.496 1.885 2.235 
4 1.995 2.513 2.980 
5 2.493 3.142 3.725 
6 2.992 3.770 4.470 
7 3.491 4.398 5.215 
S 3.989 5.026 5.960 
9 4.488 5.655 6.704 
10 4.987 6.283 7.449 
Taste XXIX—Auxiliary Heat Schedules 
By B; By Bs Bs 
October — .3H - 
November — - 4H OW 
December W .5H W W+ | W+ 
-5H -25H 
January W ~ .5H W W+ W+ 
5H -25H 
February .5H W W+ 
-5H 
March 4H .5W | .bW 
April — .3H 


May thru Sept.—None required 


W = water heat; H = house heat. 


APPENDIX A: THE M.LT. SOLAR HOUSES 


The standard house assumed in this report is similar 
to the solar houses built by Massachusetts Institute 
of Technology at or near Blue Hill, Massachusetts. 
This report is much indebted to the published re- 
ports' * § on these houses. However, several differences 
in point of view may be noted: 

1. The M.1.T. houses are “‘tight,’’ with air change in 
winter reduced to } change per hour, as compared with 
14 or 2 changes per hour in a normal house. The extra 
construction cost involved in sealing all air leaks should 
be added to the solar house heating costs. Not knowing 
what this extra cost is, | have not made the calculation, 
but the B; schedule in Table X XIX and Fig. 9 is caleu- 
lated for a house whose solar heat requirements are 
proportional to those of an M.I.T. house. The B; 
schedule amounts to 34 per cent of the annual heating 
load (house heat plus hot water) being supplied by an 
auxiliary heating unit. 

2. The M.I.T. houses are small, about 1,200 sq. ft. 
in floor area. This is understandable in an expensive 
experiment. I have assumed for calculation, however, 
that since solar heating requires a substantial capital 
investment, a house for an above-average income level 
is indicated. 

3. Heat storage is not considered in the M.I.T. 
studies as a means of reducing the necessary collector 
area because only overnight storage (or storage for two 
or three days) is contemplated. The M.I.T. reports 
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ARY HEAT FOR SOLAR UNIT 
200 OIL AT 18c PER 10° BTU 
+-— GAS AT 22¢ PER 10° BTU. 
~- OFF PEAK ELECTRICITY AT 45¢ PER 10° BTU 
REGULAR ELECTRICITY AT 75¢ PER 10° 


SYSTEMS WITH NO SOLAR HEATING 
A.OiL FOR BOTH HOUSE AND HOT WATER 
B.OIL FOR HOUSE AND HOT WATER IN WINTER. 
100 OFF-PEAK ELECTRICITY FOR HOT WATER ~~ 
> IN SUMMER 
C.GAS FOR BOTH HOUSE AND HOT WATER 
D.OIL FOR HOUSE. OFF PEAK ELECTRICITY 
Naz, FOR HOT WATER. T 


| 
‘ fo $4 cost oF COLLECTOR 
q A PLUS STORAGE, PER 
; SQ.FT. 


> 


° 200 400 600 800 1000 1200 1400 
COLLECTOR AREA SQ.FT. 


Fic. 14—M.1.T. Solar House cost analysis (from Whillier). 


refer to storage capacity as proportional, rather than 
in inverse proportion, to collector area. The houses 
reported on at Phoenix':? had storage tank capacities 
of 1,200 and 1,500 gal. 

4. In the M.I.T. calculations the parameters are 
optimized by calculating the effect of different per- 
centages of auxiliary heat on the dollar value of the 
annual heat collected per square foot of collector. The 
value declines monotonically, (starting with a maxi- 
mum, for collector areas of less than 200 sq ft, of about 
33 cents per sq ft.) Since solar house heating is not eco- 
nomically competitive at an assumed fuel saving of 18 
cents per therm, they could not take for their criterion 
of optimization the maximum dollar value earned 
figured on the total collector area. They used, instead, 
maximum dollar value earned figured on the incremental 
area; that is, a relative law of diminishing returns. They 
found then that between 25 per cent and 50 per cent 
of the heat load should be auxiliary (non-solar) heat. 
This, combined with the smaller house considered, 
explains the fact that the optimum collector areas 
indicated in Fig. 14, taken from Whillier’s paper? are 
only about half that indicated for B; in Fig. 9 (about 
850 sq ft). In other words, the M.I.T. houses, con- 
sidered in the terms of this paper, have auxiliary heating 
of about 65 per cent, 34 per cent being ‘“‘built in” by 
the tight sealing of the house. 

5. The auxiliary heat is furnished by a conventional 
furnace. In the house described above, the furnace is 


eliminated and the auxiliary heat furnished by a hot 
water heater. This changes the economic calculation 
somewhat, since the elimination of the furnace is 
figured, in New England, as a saving of $800 in favor of 
solar heating. 

6. Long-term storage is rejected flatly. Hottel' states 
that it “is undesirable, because of the very poor load 
factor on it.’’ Whillier’ states: ““The load factor on that 
part of the storage system which exceeds 1, 2 or, 3 
days’ capacity is too low.”’ This judgment does not take 
account of the fact that good insulation makes the 
system more efficient all through the year by reducing 
daily heat leakage, not just during exceptionally long 
periods of cold, cloudy weather. However, as the calcu- 
lations show, the collector area reduction achieved by 
efficient storage is quite modest. This is especially true 
in Blue Hill, Massachusetts. (See Fig. 9.) 

A direct comparison between the M.1I.T. results and 
those presented in this paper cannot be fairly made 
because of the above-mentioned differences. In addition, 
there are slight differences in the calculated solar energy 
collected per square foot in each month. However, 
Fig. 14, reproduced from Whillier’s paper,? shows 
satisfactory agreement. The annual fuel oil bill, when 
there is no solar heating (collector area zero) is seen to 
be about $180; at $2 per sq ft, solar heating never re- 
duces the annual heating bill below about $210. 

The break-even point would seem from Fig. 14 to be 
about $1.50 per sq ft, as compared with $1.86 which was 
calculated for B; by the optimizing method explained 
above. 


APPENDIX B: STORAGE TANKS 


We define a dimensionless shape factor, ¢ as the ratio 
of the surface area of the tank to the % power of the 
volume. (See Equation [6|.) A spherical tank, which 


has, of course, the lowest value of ¢ (4.836) costs ap- 
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Fic. 15—Shape factor vs. ratio of height to diameter of right 
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proximately twice as much as a cylindrical tank. A 
“square” cylinder, i.e., one whose diameter equals its 
height, has a o value of 5.536; thus for the same storage 
efficiency a cylindrical tank requires about 15 per cent 
greater thickness of insulation than a sphere. 

To avoid the expense of on-site assembly of the tank, 
the diameter must be kept below the maximum which 
can conveniently be shipped, i.e., about 11 ft. Thus, 
tanks larger than 10,000 gal are usually made longer 
than they are wide. Fig. 15 shows the effect of the shape 
factor of varying the proportions of the cylinder; the 
abscissa is the ratio of height to diameter. On the right 
hand margin the shape factor has been normalized, with 
the sphere taken as unity, so that one can read directly 
the percentage penalty for cylinders of different pro- 
portions. The sphere falls too far below the curve to 
appear, but the hemisphere is indicated. 

The larger the capacity of the tank, the larger the 


surface area, and hence the greater the heat leakage 


per day for the same thickness of insulation. The ques- 


tion arises whether storage efficiency increases or de- 


creases with tank capacity. Storage efficiency increases 
with tank capacity, the period of useful storage being 
longer for a 10,000-gal tank than for a 3,000-gal tank. 
lig. 8 shows that storage efficiency, i.e., percentage 
energy available along any constant time line, increases 
monotonically as L decreases. As can be seen in Equa- 
tion [8|, the L factor varies inversely as the one-third 
power of the volume. The same conclusion can be 
drawn from Table XXVIII where the same insulation 
factor, savy 3 X 10°-*, corresponds to lower L values as 
one moves across the table from left to right. 

Steel tanks, lap-welded and painted with a rust- 
inhibitor-like calcium plumbate, are preferable to 
wooden tanks. Tidewater red cypress is the best wood 
for water tanks; prices are not quite competitive with 
steel tanks. If one insulates the outside of a wooden 
tank so that it cannot “‘breathe,”’ the life of the tank is 
seriously shortened. Moreover, with wooden tanks 


” 


there is a need—from time to time—to tighten the 
hoops which hold it together in order to prevent water 


leakage. 
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Design of a Sea Thermal Energy Power Plant 


By Asa E. Snyder 


A unique 8000-kw gross power and sea water 
conversion plant is described which uses the 
temperature difference of the ocean to produce 
usable energy and fresh water. This plant has 
been designed by the research and development 
department, Clearing Division of U.S. Industries, 
Inc., for installation along coastlines in a semi- 
tropical area. Estimated cost for the sea thermal 
energy plant is 2.5 to 3 million dollars, including 
building and erection costs. This is equal to $310 
to $375 per kw, amounts which are competitive 
with fossil-fuel burning plants in remote areas. 


INTRODUCTION 


The thermal energy of the sea is actually latent solar 
energy which has been stored within the large bodies 
of water upon the earth’s surface. Almost 35,000 times 
the existing annual energy consumption of the world is 
delivered annually to the earth in the form of solar 
radiation. Over seventy per cent of the earth’s surface 
is covered by the oceans; thus, a major percentage of 
the world’s supply of solar energy is deposited in the 
seas. 

Actually, however, only the surface water or top 
lavers of the ocean are warmed. At relatively shallow 
depths in the ocean, no more than 1200 ft, the sunlight 
does not penetrate, and the ocean water temperature 
is 40°F or less, but surface water temperatures of 80° 
to 90°F are quite common in the tropic regions. 

Based upon our existing knowledge of thermal power 
cycles, we know that any such system requires: 

(1) a source of heat, or a hot body, 

(2) a receiver of heat, or a cold body, 

(3) a working medium, 

(4) a means of abstracting or adding heat and a 
method of producing work. 

The temperature difference between the surface and 
the depths can be utilized as the source and receiver 
respectively and is therefore a potential means of ob- 
taining almost unlimited amounts of energy. The water 
itself can be utilized as a working medium, and only an 
economical and practical means for extracting work is 
required. The objective of this paper is to describe the 
fundamentals of a plant designed by the Research and 


USI Clearing Division of U. 8. Industries, Ine. 


Development Department of Clearing Division of 
U.S. Industries, Inc., that will utilize the temperature 
difference of the oceans to produce useable energy and 
fresh water. 


Brief Background of S.T.E. 


The initial impetus to the study of sea thermal 
energy Was given by Monsieur Georges Claude be- 
tween 1925 and 1935 when he carried out tests on the 
installation of a submarine cold water pipe from a 
coastal power station in Cuba, and of a vertical sub- 
marine cold water pipe beneath the floating ‘Tunisia’ 
power station. The accomplishments of his experiments 
were: 

(1) The suecessful operation of a steam turbine with 
a difference of 27°F. 

(2) The design of an efficient air extractor. 

(3) The demonstration of the problems in protection 
from wave action and laying the cold water inlet pipe 
in the open ocean. 

In 1941 the French Government created an inter- 
ministerial organization to continue the research of 
Claude. In 1948 the semi-official company ‘Energie 
des Mers’” was formed by the French Government to 
continue research and study in this field. 

In the United States, the Office of Saline Water, 
U.S. Department of Interior, has sponsored design 
studies of flash-type evaporators which might utilize 
energy in bodies of water. The basic studies for work 
of this nature were conducted under the Sea Water 
Conversion Program at the University of California. 

The Research and Development Department, Clear- 
ing Division of U.S. Industries, Inc., has actively 
entered into a program of design research and study in 
the field of sea thermal energy and desires to stimulate 
interest in this form of energy conversion for use in the 
semitropical coastal areas of the world. 


Distribution of Temperature and Suitable Areas 

for Installation of S.T.E. 

The general distribution of surface temperatures of 
the oceans is given in Table I. 

Analysis of the data given in Table I shows that the 
25°C (77°F) clime falls within the belt of 20° North 
and South latitudes. As stated earlier, at depths of 
1200 ft or more the temperature is sufficiently low, 
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TaBLe I—Surface Temperature of the Ocean! 


" Pacific (°C Indian (°C) 


Latitude Atlantic 


70°-60° N 5.60 

60°-50° N 8.66 5.74 

50°-40° N 13.16 9.99 

40°-30° N 20.40 18.62 

30°-20° N 24.16 23.38 26.14 
20°-10° N 25.81 26.42 27.23 
10°-0° N 26.66 27.20 27.88 
10°-0° S 25.18 26.01 27.41 
20°-10° S 23.16 25.11 25.85 
30°-20° S 21.20 21.538 22.53 
40°-30° S 16.90 16.98 17.00 
50°-40° S 8.68 11.16 8.67 
60°-50° S 1.76 5.00 1.63 
70°-60° S 1.30 1.30 1.50 


with minimal variations, in all oceans. Therefore, the 
appropriate areas for consideration of the feasibility of 
installation of S.T.E. Plants is established by the 
world’s water temperature distribution. 

Countries or areas which can reasonably expect to 
utilize this type of power plant can be listed as follows: 
(1) Central America and the Caribbean Islands 
(2) Northern Coast of East and West South America 

(3) West Coast of Africa—from Mauritania South 
to the Congo 

(4) East Coast of Africa—from Somaliland South to 
Madagascar 

(5) South Coast of Arabia 

(6) Ceylon 

(7) Indonesia, New Guinea, and Borneo 

(8) Phillipine Islands 

(9) Taiwan 

(10) Australia 

More detailed studies of suitable locations have been 
made and are contemplated by Clearing R&D, and are 
available for discussion with interested parties. 


Geological Considerations in Plant Location 
The power plant must be close to deep water which 
has unrestricted communication with the oceans. A 
steep shore, such as that usually found along a voleanie 
or coral island, spits of land, natural peninsulas, or bars 
of river-borne silt are desirable. If at all possible, a 
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Fig. 1—Available energy of sea thermal energy (S. T. E.) 
plant. 
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Fic. 2—Simplified schematic of 8. T. E. plant. 


near-by submarine canyon which is approximately at 
right angles to the shore line is desirable. Excessively 
rugged or rocky profiles along the bottom are to be 
avoided, but are not limiting. Fig. 4 illustrates a plan 
for a typical plant. 

Principle of Sea Thermal Energy 

Fig. 1 illustrates a simple temperature entropy dia- 
gram for water. Observe that cooling the temperature 
of water vapor already heated to 80°F (Point 1) to 
40°F (Point 2) lowers the vapor pressure to that ob- 
tainable at 40°F. The vapor pressure at 80°F is 0.5069 
psi, whereas the vapor pressure at 40°F is 0.1217 psi. 
We thus have a small but decided pressure drop, which 
can be used to extract power from the seas. 

Fig. 2 schematically illustrates a simplified basic 
system for extracting sea thermal energy in the form of 
useful mechanical work. Referring to that figure the 
process is described as follows: 

(1) Heated surface water is lead into a chamber 
where the water is evaporated at the higher tempera- 
ture. 

(2) This chamber is connected to a turbine, which 
extracts mechanical work as the pressure drops through 
the turbine. 

(3) The exhaust from the turbine then enters a con- 
denser where the vapor is condensed into liquid. 


THE CLEARING SEA THERMAL ENERGY 
POWER AND SALT WATER CONVERSION 

PLANT DESIGN 
As part of a long-range program, the Research and 
Development Department, Clearing Division of U.S. 
Industries, Inc., has investigated and designed an 
8000-kw gross power and sea water conversion plant 
suitable for location along a coast line in a semitropical 
area of the world. 

The sea thermal energy power and sea water con- 
version plant described in this paper is a unique type 
of plant; no installation exactly similar to the proposed 
system exists in the world today. The development of 
this type of power plant has awaited the advent of in- 
dustrialization of those tropical land areas where 
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favorable thermal differentials exist in the nearby 
oceans. As mentioned earlier, a sizable backlog of aca- 
demic and industrial programs on components for the 
sea thermal energy (S.T.E.) plant does exist today; and 
one similar installation located on the African Ivory 
Coast was actually carried into a construction phase 
before being discontinued by the French Government 
because of the close proximity of low-cost hydroelectric 
power. No extensive development of component hard- 
ware is required today to put a sea thermal energy 
design into practice. 


Purpose and Scope of Described Design 

A set of hypothetical requirements is assumed for 
this paper. The installation is assumed to be part of a 
process-type facility situated in a desirable location 
(close to deep water) along a shore line on open ocean. 
This installation will be required to produce large 
quantities of power, fresh water, and brine; a 5000-kw 
net output is required as part or all of the output of the 
facility. 

Other sizes of power plants could be studied, pro- 
ducing any preselected combination of power, fresh 
water and brine, depending upon the requirements of 
the specific installation. 


Design Objectives 

The plant is designed for a constant net output of 
5000-kw, and the potable water and brine production 
rate is to be established by this. Design emphasis for 
the entire installation is placed upon high rates of 
availability and reliability. Maintenance costs are to be 
minimal with major overhaul periods of two years. 
Initial cost of the power plant, including all structures, 
power generation, piping, and storage facilities will not 
exceed $375 per kw based on the gross kw output. In 
order to accomplish this, efficiency is to be sacrificed in 
favor of design simplicity and low initial cost. When- 
ever possible in the basic design, consideration is to be 
given to manufacture or subassembly of components 
at erection site to reduce freight differentials and to 
utilize effectively local labor rates in remote areas. 


General Description of System 

A simple schematic of the proposed system is illus- 
trated in Fig. 3. The water is initially gravity fed 
through an inlet open to the ocean, which is essentially 
a ditch cut by bulldozer leading into a solar heating 
flat. The solar heating flat is a scientifically designed 
“black body” heat collector, where the inlet water is 
heated by solar radiation to a maximum temperature. 
The basic design of the system is such, however, that 
under adverse atmospheric conditions it will continue 
to work at reduced efficiencies, i.e., night or cloudy 
days. 
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Fig. 3—Schematie of 8. T. E. installation. 


The water then flows to the deaerator which extracts 
air and other dissolved gases, which would decrease 
the efficiency of the plant, from the sea water. The 
deaerated sea water passes into the evaporator and is 
partly flashed into steam because of a high evaporator 
vacuum. A brine-return pump holds the salinity in the 
evaporator to a level that prevents salt precipita- 
tion. 

The vapor then flows through a reaction turbine, 
which drives a generator, and then passes to the lower 
pressure zone in the condenser, where cold sea water 
circulating through the condenser tubes cools and con- 
denses the steam to potable water. The condensate 
pump distributes the fresh water to suitable storage 
basins. 

The cold sea water is pumped from an approximate 
depth of 2000 ft. A specially designed intake pump is 
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Fic. 4—Plan view of S. T. E. installation. 
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Fig. 5—Cross-section of S. T. E. installation. 


installed from the bay to the inlet depth. The pipe is 
designed in sections and is manufactured from several 
different materials. The sections passing through the 
surf are corrosion-proof steel welded pipe, encased in 
concrete protective jackets, and buried under a surf- 
breaking rock fill. The next sections of the pipe are 
constructed from fiberglass and epoxy resins with inter- 
mittent flexible rubber bellows sections. The inlet to 
the pipe consists of a special water-ballasted tank-like 


device with foul-proof inlet screening. 


Description of the Plant 


A cross-section view of the proposed power plant is 
shown in Fig. 5. Fig. 4 is a plan view of the total instal- 
lation. Table II lists the assumptions used in making 
the caleulations, and Table III shows the ecaleulated 
results. 

The solar heating flats are large black-topped reces- 


sions where water is heated by solar radiation under 


Assumptions for Technical Calculations of 8.T.E. 
Power Plant 


TaABLe II 


1) Net power 5000 kw 
2) Inlet sea water temperature 80°F 
3) Solar heating flat outlet tempera 


ture 90°F 
1) Evaporator steam temperature S4°F. 
5) Cooling water inlet temperature 10°F 
6) Condensate temperature 50°F 
7) Density of fresh water @ 50°F 62.4 lb/cu.ft. 


8) Density of sea water @ 40°F 64.0 lb/eu-ft. 


9) Specific heat of water 1.0 B/Ib.°F. 
10) Deaerator steam loss 1°, of Gross Steam 
Flow 

11) Intake pipe (total length 1000 ft 
inlet depth 2000 ft 

12) Assumed component efficiencies 
AC generator efficiency 95"; 
Turbine mechanical efficiency 90°; 
Turbine stage efficiency 85, 
Coolant pump efficiency 90%; 
Condensate pump efficiency 80%, 
Brine return pump efficiency 90%; 


Tas_e for 5000-kw Net T. E. Power 
Plant 


Fluid qualities 

Inlet flow to evaporator 

Brine return 

Vapor loss in deaerator 

Potable water flow 

Cooling water flow 
Pumping and auxilary power requirements 

Cooling water pump 


192,000 gpm 

| 190,900 gpm 
11 gpm 

1,123 gpm 
140,600 gpm 


1,660 kw 


Condensate pump 22 kw 
Brine return pump 1,300 kw 
Vacuum pump 1S kw 
Total 3,000 kw 


Gross kw output 8,000 kw 


Electric generator 


Speed 450 rpm 
Voltage 4160 
Type 3 Phase—60 cycles 


Steam flow 
Enthalpy drop (Ideal) 
Enthalpy drop (Actual) 
Net flow 


66.0 Btu/Ib 
51.2 Btu/Ib 
156 lb/see 
(9357 lb/min) 
Condenser 


Heat extracted 9.63 & 10° Btu 


min 
Log mean temperature difference 4.97°F 
Friction head loss 56.1 ft 
Area 188,000 sq ft 


3° OD x 049 Gage 
239,000 ft 


Tubing size 
Tubing total length 


Tubing spacing 5’’ centers stag. 
Condenser Dimensions 

Inner diameter 10 ft 

Outer diameter 50 ft 

Height 19 ft 


Condenser weight with baffling 192 tons 


Tu rbi rie 


Speed 450 rpm 
Hub diameter 30.8 ft 
Tip diameter 33 ft 
Mean diameter 31.9 ft 
Bucket height 1.08 ft 


Maximum tip velocity 776 fps 


Inlet Pipe 


Inside diameter 6.30 ft 
Total length 4000 ft 
Inlet depth 2000 ft 


Solar Heating Flat 


Total area Approximately 9 


million ft 
(83000 ft x 3000 
ft) 


favorable conditions. At night the water is received at 
ocean surface temperatures, and the entire plant oper- 
ates at reduced outputs. 

As may be seen from Fig. 5, the heated water from 
the solar flats enters the deaerator chamber through a 
simple vacuum seal, after which it is cascaded over a 
series of vertically spaced, perforated circular trays. 
The water is broken into fine droplets in passing 
through the perforations, which assures complete de- 
aeration because of the large wetted area exposed. The 
pressure in the chamber is maintained by a vacuum 
pump at a value somewhat above the steam vapor 
pressure to minimize the steam loss. This is estimated 
to be 1 per cent of the net steam flow. 

The deaerated water next flows into the evaporator, 
where it partially flashes into steam due to the lower 
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pressure in the evaporator. Mechanical baffling sepa- 
rates the vapor from the liquid and assures that only 
dry steam reaches the turbine. The unflashed liquid is 
returned to the solar flats for reheating by the brine 
return pump, which, of course, works against an 
atmospheric head. 

After expansion in the turbine, which may be con- 
structed of lightweight metal or inexpensive plastic 
materials because of the low temperatures, the steam 
passes into an annular condenser consisting of tubes 
arranged normal to the steam flow and baffled to force 
the steam to take repetitive passes through the tubes. 
The annular condenser represents the most economical 
construction for the large heat-transfer value required. 

The condensate is collected in the bottom of the con- 
denser and removed by the condensate pump, which 
also works against an atmospheric head. 

In the event that potable water is not desired, a mix- 
ing type of condenser can be incorporated in the design 


instead of the surface condenser at some saving in 


initial cost. 
The turbine is vertically mounted, supported within 
simple water-lubricated bearings. The bearing housing 


supports are rigidly mounted to the main shell. The 


turbine shaft connects to the generator, which is inde- 


pendently mounted and isolated on the top of the shell 
and connected to the turbine by means of a flexible 


coupling. A rotary vacuum seal maintains the vacuum 


within the shell. 


Design of Intake Piping 


The inlet piping system is illustrated in Fig. 6. The 
first section of the pipe, extending from the shore to a 
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Fic. 6—Inlet piping system of 8. T. E. installation. 


depth of 50 ft, is constructed of corrosion-resistant 
welded steel pipe encased in a protective concrete 
jacket. This pipe is laid in an excavated trench and 
then covered with stone fill as a surf breaker. The next 
section of the pipe is a flexible pipe designed to follow 
the ocean floor to a depth of 2000 ft. 

The cooling water inlet pump is a Kaplan-type 
pump designed to handle 140,600 gpm. In order to 
operate efficiently this pump must have a positive head 
and is therefore located below the low-tide mark. 

The second section is composed of alternate sections 
of flanged fiberglass pipe and steel and rubber bellows. 
The fiberglass sections are anchored to the floor by 
means of cast concrete ballast blocks, except for the 
initial sections which consist of controlled buoyancy 
steel tanks. The inlet consists of a specially designed 
controlled buoyancy steel tank, containing nonfouling 
screens with self-cleaning features. 

At the time of final construction, the entire assembly 
is floated into position using the buoyant tanks and 
special rafts. The first section is connected to the second 
section or deep part of the piping by skin divers using 
scuba gear, and the assembly is sunk and lowered into 
position. Cables from the raft convoy permit a gradual 
positioning of the final assembly into place. 


Design for Prevention of Marine Fouling 


The growth of marine microérganisms, both plant 
and animal, could be a continual source of trouble 
unless special precautions are taken by the designer to 
prevent this from occurring. Very little is actually 
known about the degree of fouling to be expected in 
deep water inlets. We propose that antifouling measures 
based upon existing power plant practices be in- 
corporated in the design. Two basic control measures 
will be incorporated: 

(1) use of velocity to control growth, 

(2) control by chlorination. 

Actual results of operating steam power plants indi- 
cate these provisions will offer adequate control at a 


minimum of cost. 


Starting and Control of the Plant 

If the proposed plant is to be used as a part of a 
power installation, i.e., an additional steam or diesel 
generation installation is available, then the generator 
of the S8.T.E. power plant may be converted to a 
motor, and at the same time additional power may be 
supplied to the auxiliaries in order to start the plant. 
After the plant becomes self-supporting, the outside 
power supply is cut off, and the plant continues oper- 
ating. 

Should the 8.T.E. plant be the sole power source, 
then a small standby auxiliary plant, either diesel or 
steam, will be required to start and activate the system 


‘ 
4 
4 
| 
159 | 
: 
0) 
40 ft. 
Sa 
fers: 
4 
/ 
/ 
Inlet 
2 
53 
4,4 


and also to act as an emergency energy source during 
maintenance intervals. 

The 8.T.E. power plant, once started, can operate 
at constant maximum output during daylight hours 
and at reduced outputs during evenings, rainfall, or 
other adverse solar conditions. If the power require- 
ments fluctuate and the 8.T.E. is the sole power source 
at the mill, a special governor will be built into the 
turbine to hold speed constant over a varying load. 
Essentially, therefore, the addition of control features 
is dependent upon the over-all requirements of the 


installation. 


ECONOMICS OF SEA THERMAL ENERGY 
POWER AND SEA WATER CONVERSION 


The size of equipment required for a power plant of 
the type described in this paper is relatively large, but 
the low temperatures of the cycle permit the use of 
very low-cost materials The initial capital investment 
of a S.T.E. power plant will not exceed that of 
an equivalent fossil fuel-burning thermal power plant, 
and the operating and maintenance costs of a sea 
thermal energy installation will be minimal. 

In addition to producing electrical power, the S.T.E. 
power plant produces two by-products which make the 
investment additionally attractive. As shown in 
Table III, the potable water by-product will approxi- 
mate 1,123 gpm for an 8000-kw gross installation. 
Furthermore, the installation will produce salts at 
various production rates depending upon the size of 
available evaporating flats. The true value of the 
potable water and salts must be considered in estab- 
lishing the actual operating expense of the plant. 

The estimated cost for the 8000 gross kw Clearing 
S.T.E. installation is 2.5 to 3.0 million dollars, inelud- 
ing building and erection costs. On a dollar per installed 
kilowatt basis this amounts to $310 to $375 per kw, 
which is competitive with installation costs of a com- 
parable fossil fuel-burning thermal power plant  in- 
stalled in a remote area.* 

The Clearing 8.T.E. 8000-kw power and salt water 
conversion plant produces 1.62 X 10° gal of fresh water 
per day. In order to compare the basic cost of the 
S.T.E. plant with a more conventional water conversion 


plant, a simple present-worth engineering economy 
study has been made comparing the 8.T.1. plant with 
a multistage flash plant heated by steam. 

For the Clearing 8.T.E. plant installed at a site in 


an underdeveloped area for a 20-year life: 
First cost (estimated) $2,480,000 
P.W. of annual disbursements based on 40,000 


per vear at 4% (40,000 X 13,590) 543,000 


$3 ,023 , 000 


P.W. of all disbursements for 20 years 


However, the Clearing 8.T.E. plant also produces 
5000-kw of power, which can be sold. Then: 
P.W. of all receipts for electrical sales based 


on §$.05 per kw hour at 4°; ($2,190,000 X 
13,590) $29, 800 , 000 


The present worth of net receipts therefore 
$29,800,000 minus $3,023,000, or $26,777,000. 
teferring to Reference 2, the first cost of a multi- 
stage flash steam heated plant for 1.0 X 10° gpd is 
$1,534,900 with disbursements of $1015 per day. For 
a basis of comparison (1.62 XX 10° gpd plant of the 
same design), the cost and disbursements are assumed 
to be increased by 25 per cent; thus, the first cost is 
$1,920,000, and the disbursements are $1270 per day. 
Then on a present worth basis: 
First cost $1,920,000 
P.W. of annual disbursements at 4°, ($261,000 


13.59) 3,540,000 


P.W. of all disbursements for 20 years $5, 460, 000 


There is therefore a large advantage for the S.T.E. 
installation of ($26,777,000 plus $5,460,000) $32,237,000 
over the more conventional type of distillation process. 
Such economic advantages cannot easily be ignored. 
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All variables influencing the efficiency of a 
flat-plate solar heat collector as a heat exchanger 
can be combined into a single ‘‘efficiency fac- 
tor.”” These efficiency factors are more or less 
design constants of the particular collector de- 
sign, and are only slightly influenced by operating 
conditions. Consequently they are extremely 
convenient for use in accurate design and per- 
formance calculations. The full mathematical 
derivations are presented for several of these 
efficiency factors for various types of collectors, 
together with graphical data and examples of 
their use. 


PRINCIPAL SYMBOLS USED 


A Intensity of atmospheric radiation upon a horizontal 
surface, Btu/sq ft-hr 

Ay Fin cross-section, sq ft 

C.. Specific heat of circulating fluid, Btu/pound-°F 

d (do).. External (internal) tube diameter, ft 

F .Fin efficiency factor, dimensionless 

F’ Collector plate efficiency factor, dimensionless 

Fr Collector plate efficiency factor, dimesionless 

G.. Flow rate of circulating fluid, pounds/hr-sq ft of 
collector 

ha Convection coefficient, collector plate-to-air, Btu/sq 
ft-hr-°F 

h. Convection coefficient, tube-to-circulating fluid. 
Btu/sq ft-hr-°F 

h Simplified radiation coefficient, Btu/sq-ft-hr-°F 

I, .. Intensity of direct beam sunshine upon a surface nor- 
mal to the rays, Btu/sq ft-hr 

k ... Thermal conductivity, Btu/ft-hr-F° 

L. Fin length, ft 

M Fin thickness, ft 

Q Heat transfer rate, Btu/hr 

q Heat transfer rate per sq ft of collector surface, Btu/sq 
ft-hr 
qa total sunshine absorption rate 
ga...-ineident diffuse sunshine upon collector 
qi incident total sunshine upon collector 
qy-...total loss rate of collector 
qu....average useful heat collection rate 


* This paper was supported in part by a grant from the 
National Science Foundation (NSF-G6293). 
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Qu-o.. useful heat collection rate at entering section of 
collector 
Qu-c-.useful heat collection rate at point x 
R Net radiation from a surface to the sky, the surface 
being an exposed thermally black horizontal plane at 
air temperature, Btu/sq ft-hr 
t(T)....Temperature, °F (Rankine) 
T...air temperature 
t plate temperature below bond between plate 
and tube 
t, T>..plate temperature 


tube temperature 
h, fluid temperature at a given point 
: temperature of fluid entering collector 
Us, Over-all heat transfer coefficient for loss through bot- 
tom of collector, Btu/sq ft-hr-°F 
U. Over-all heat transfer coefficient for heat transfer 


through bond between plate and tube, Btu/hr-sq ft 
of bond-°F 


Ur, Over-all loss rate coefficient for collector, Btu/sq 
ft-hr-°F 

w.. Tube spacing, ft 

a Average absorptivity of surface for total solar radia 
tion incident upon a surface, dimensionless 

€ Average emissivity or absorptivity of surface for long- 
wave radiation, dimensionless 

6 Altitude of sun relative to collector plane, degrees 

o Stefan-Boltzmann constant, Btu/sq ft-hr-R 


INTRODUCTION 

In problems involving the design or testing of flat- 
plate solar heat collectors it is very worthwhile to 
calculate the ratio between the actual useful heat 
collection rate of a given design and the useful heat 
collection rate which would have been attained had 
the entire collector surface been at the temperature of 
the circulating fluid. There are two useful forms of 
this ‘‘plate-efficiency factor”’: 


actual useful heat collection rate 


== 
‘ useful heat collection rate attainable with 
entire collector surface at average fluid tempera- 
ture 
a actual useful heat collection rate 


useful heat collection rate attainable with 
entire collector surface at temperature of fluid 
entering collector 
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These factors F’ and Fx, are particularly useful 
because, for many types of solar heat collectors, they 
are essentially design constants determined by the 
construction details of the collector, and do not de- 
pend greatly upon such operating conditions as tem- 
perature, sunshine intensity, or wind velocity. Efh- 
ciency factors of this sort were first used in problems 
of solar heat collector design by Hottel and Woertz,' 
later by Whillier,” and by Hottel and Whillier.’ The 
use of such factors tremendously reduces the empiri- 
cism associated with the design of solar heat collectors. 

For many types of flat-plate collectors, and in par- 
ticular for water-heating collectors formed of tubes 
bonded to a metal sheet, the efhiciency factors for a 
given design can be quickly calculated to a high degree 
of accuracy. Nevertheless, and despite their obvious 
advantages in collector design, these factors are not 
utilized as much as they should be. 

These notes present a full derivation of several of 
these “plate-efficiency factors,” with the thought 
that such a full presentation will be helpful in furthering 
the understanding and use of these basic concepts and 
heat-transfer equations. Our end results will be found 
to be in agreement with those which have been stated 
by Hottel and Whillier.” In general we should note 
that the calculation methods used by the Solar Energy 
Laboratory in attacking problems of collector design 
are closely similar and strongly indebted to those of 
the Massachusetts Institute of Technology group. 
The chief difference between our ‘“‘collector equations” 
and theirs is that we make a small but (we think) 
worthwhile improvement in allowing for the effects of 
atmospheric radiation. 

These notes are divided into four main sections: 

(1) Basic concepts in equation form. 

(2) The temperature distribution and over-all 
heat transfer efficiency of a uniform rectangular fin. 

(3) Collection efficiency factors for collectors having 
uniform flow beneath a flat plate. 

(4) Collection efficiency factors for collectors having 
uniform flow in tubes bonded to a flat plate. 

BASIC CONCEPTS IN EQUATION FORM 

The fundamental steady-state equation of any solar 
heat collector amounts to the simple statement that 
the rate of useful heat collection is equal to the rate 
at which sunshine is being absorbed less the total loss 
rate: 

qu = Ga — (1) 


A good place to start in analyzing the performance 
of a solar heat collector is to consider the heat trans- 
fer from the upper surface of an exposed horizontal 
plate which is imagined to be held at constant tem- 
perature. We can, if we like, look upon this problem 


as that of setting up the performance equations for a 
hypothetical collector having the following charac- 
teristics: 

(a) The collector plate is horizontal, exposed, and 
at constant temperature over its entire surface. 

(b) There are no thermal gradients between collec- 
tor plate and circulating fluid. 

(c) The bottom of the collector is perfectly insu- 
lated. 

We then set up the performance equations for this 
idealized collector, assuming various hypothetical 
cases of increasing complexity. 


Case I: Plate Surface Thermally Black, Plate 
Temperature Equal to Air Temperature 


(a) Operating at night: 


Qu = [2] 
Note: 

R=oT,' —A [4] 

(b) Operating during the day: 
Qu = Ya R [5] 
ga = 4: 
Gi = Ga + sin 6 


Case Il: Real Plate Surface, Plate Temperature 
Equal to Air Temperature 
A real plate surface will have average solar absorp- 
tivity coefficients and average emissivity and ab- 
sorptivity coefficients in the long-wave region (2 to 
50 mu), all of which must be reckoned with. We assume 
that long-wave emissivity and absorptivity can be 
considered equal. The resulting equations: 
(a) Operating at night: 
= [S| 
(b) Operating during the day: 
qu = ag; — R [9] 
Case IIL. Real Plate Surface, Plate Temperature 
Differs from Air Temperature 
Here we must reckon both with absorptivity- 
emissivity effects and with convective heat transfer 
to the surrounding air. 
(a) Operating at night: 
Qu = —eR — (eh, + ha) (tp — ta) {10} 
(b) Operating during the day: 


eR — (eh, + ha) (tp — ta) (11) 


In Equations [10] and [11] above, the long-wave 
radiation loss rate from the plate is seen to be split 
into two components, a procedure which is convenient 
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The heat gain, in Btu per hr, to the fin through the 


/ 
4 area Ay dx on top of the fin is: 


Surrounding dQ d't 
Air Temp, t, = dx = —kMAy dx [16] 


But this same heat gain in Btu per hr is: 


= Gyr Ay-dz 


Substituting between [14], [16], and [17] and letting 


= 9 

b? = U,/kM: 

*j Fig. 1—Heat flow in a uniform rectangular fin. Py R 

2 ¢ — 

b (: ) [18] 
4 in later calculations. The equations do yield the correct = UL 

= long-wave radiation loss, which is: This has the general solution: 

4 

* Qt-radiation = e(oT', A ) y da eR 

= — oT, + R) [12a] Tat 
= eh,(tp — ta) + €R [12b} To determine the value of the constants C; and C, 


in the preceding equation, we apply the boundary 


We find that the “simplified radiation coefficient” 
h, is very convenient for use in calculations. A plot 
of this coefficient is given in Fig. 4. dt 

5 —" and ft, 


THE TEMPERATURE DISTRIBUTION AND 
OVER-ALL HEAT TRANSFER EFFICIENCY 
OF A UNIFORM RECTANGULAR FIN 
_ — eR 


This problem occurs repeatedly in the design of (t, — ta) —— 
solar heat collectors, as well as in the design of many C=C; = [21] 
other types of heat exchangers. The properties of fins a 
are discussed in most of the standard heat transfer 
texts, but, for reasons of brevity, the pertinent equa- cosh be = 4(e* + &”*) [22 
tions are seldom derived in full. Herewith, in “collector 
nomenclature,”’ a full derivation of the case of a rec- 
tangular fin of constant cross-section, with end effects 


conditions: 


[20] 


Application of these boundary conditions will be 
found to yield 


We note that: 


Substituting from Equations [21] and [22] into [19] 
and rearranging, we obtain the solution for the tem- 
perature distribution along the fin: 


ignored. 

Consider a fin of length L as shown in Fig. 1. Sun- a (= i) 

shine is being absorbed along the entire fin at rate cosh bL 23] 
qa Btu per sq ft-hr. Heat is flowing through the fin qg-~ a cosh be 
as shown, and at the base of the fin the temperature is ’ U; cosh bL 


held constant at ¢,;. Surrounding temperature is ¢, , 
and the over-all loss rate coefficient is U;,. Assume 
that the width of the fin (i.e., in the direction per- 
pendicular to the figure) is Ay. Assume also that all 
heat flow within the fin is essentially parallel to the 


The average useful heat collection rate of the entire 
fin can be determined conveniently by noting that 
this average useful collection rate is equal to the total 
heat flow (Btu per hr) at the base of the fin, divided 
by the surface area of the fin: 


X-axis. 
The cross-section of the fin through which heat is er Qr—1 (24) 
flowing is: LAy 
A; = MAy (13] From Equation [15]: 
At a point x, the rate of useful heat collection in Qran _ kM dt, [25] 
Btu per sq ft-hr is: LAy L, Gout 
[14] Performing the indi¢ ated differentiation, (using 
the value of ¢, found in Equation [23] rearranging, 
The heat flow through the fin in Btu per hr is: and remembering that 6° = U,/kM: 


Q = —kAy; dz = [15] qu = bL [da eR l t.)] [26] 


- 
| dx 
At X=0 X= 
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right hand side of the above equation gives the useful 
heat collection rate which would be attainable with 
the entire fin at the temperature of the base of the 
fin. Hence the fin efficiency factor of the actual fin 
is the factor by which the bracketed expression is mul- 
tiplied : 
tanh bL 
bL 


3 


[27] 


It is important to note that the fin efficiency factor 
is not a function of the intensity of absorbed solar 
radiation or of the temperature level of operation, 
except insofar as these variables influence UU’, . 

Two other quantities which are of interest in prob- 
lems involving fins are the maximum fin temperature 
and the average fin temperature. 

The maximum fin temperature is readily obtained 
from Equation [23]: 


cosh bL 


Ga eR l ) 
* (« U, ) ( cosh bL 


The average fin temperature is also obtained from 


[28] 


Equation [23]: 
| al 
-average = dr [29] 


Performing the indicated integration vields: 


tanh bL 


(: #) (1 _ tanh 
bL 


COLLECTION EFFICIENCY FACTORS FOR 
COLLECTORS HAVING UNIFORM FLOW 
BENEATH A FLAT PLATE 


[30] 


Consider the situation illustrated in Fig. 2. The 
temperature of the circulating fluid, as it enters the 
collector, is ¢,... The fluid flows uniformly under the 
collector plate and in thermal contact with it. The 
fluid picks up heat and rises in temperature as it flows 
along under the plate. At a distance x from the enter- 
ing edge of the plate, the plate temperature is ¢, , and 
the temperature of the fluid underneath is ¢,, . 

To avoid complications due to tilt, glass covers, ete., 
all of which are extraneous to the present discussion, 
we assume the collector plate to be exposed and hori- 
zontal. 

Then, at the distance x from the entering end of the 
collector, the local useful heat collection rate may be 


Now we note that the bracketed expression at the 


-Plate Temp.,t, 


Fluid 


At X=0 “Fluid Temp.,ty 
Fluid Temp.,tw-o — Insulation 


Fic. 2—Circulation of fluid in a flat-plate collector. 


expressed, quite closely: 
Qu-z = Qa — ER — U1 (tp — ta) [31] 
in which 
Us [32] 


This useful heat collection rate is transferred from 

the plate to the circulating fluid underneath. Hence: 

Qu-2 = he(tp — te) [33] 

Eliminating ¢, between Equations [32] and [33] vields: 
h. 

h. + Uy 


The average rate of useful heat collection for the 


— R — [BAI 


entire collector (qg, in our nomenclature) is then: 


= ax 
| 
i 


h, 


L 


h, 
h, U, 


The expression in brackets at the right-hand side of 


the above equation is the useful heat collection rate 
attainable with the entire plate at average water 
temperature. Hence, for this type of collector: 
h. 
he + Ui 


We proceed to the determination of Fg. Rewriting 


F’ [38] 
Equation [34] as: 
= F'lqa — — U1 (tw — ta)] [39] 
The differential of the above is 
dqu-2 = —F'U, dt, [40} 


Now the rate, in pounds per hour, at which fluid flows 
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under the area Ay dx is GLAy. ( Note definition of G in calculations the bond width is assumed to be d, the 


the list of symbols). Hence the useful heat gain (Btu same as the tube diameter, and the over-all thermal 

per hr) associated with the area Ay dx may be equated conductivity of this bond between plate and tube is 

to the gain to the fluid as it flows under that area: taken as U.. Btu per sq ft-hr-°F. The fluid flow in each ‘ 
qu-2dy dx = G(LAy dt.) [41] tube is assumed uniform and in the same direction. 


Hence, at any given cross-section along the collector, 


Eliminating dt, between Equations [40] and [41], and 
rearranging: 
=— ‘dx 42 4 
Integrating Equation [42] and applying boundary 
condition at x = 0 yields: 


do 


d wd d d = 


The average useful heat collection rate is: 


= Qu-2 dx = e dx [44] 


Performing the indicated integration yields: 


iC (1 (FUL [45] Fig. 3 arrangement of collector with tubes bonded to 
FU, 
From Equation [39]: 
Qu-o = F'[qa — ER — Un — tal [46] 
Substituting this value into Equation [45]: 
GC —(F’U ,/GC) 
q. = F’ —— -—e 
[47] 
al* 
os In the above equation the expression in brackets is = 
o 
4 seen to be the useful heat collection rate attainable ; 
with the entire plate at entering water temperature. = E 
a Hence, the multiplying factor preceding it is, by 
“ definition, the over-all efficiency factor F» . 
ih U, ° 20 40 60 80 100 (20 140 160 #0 200 220 240 
ty -°F 
flow rate factor” Fic. 4—The “simplified radiation coefficient”’ h,. This coeffi- 
Tha cient is convenient for use in radiative heat transfer 
Phe above expression is gen¢ rally applicable to any equations of the form q = C1h,(ts — t,), where C, is some 
type of collector for which Equations [39] and [41] are constant involving surface emissivities and shape 
factors. 
valid. It may be noted that the expression can be con- -_— 
sidered as being made up of two factors: F’ and a : 
flow rate factor” (the latter factor being also a func 


tion of F’) 
A plot of Equation [48] is presented in Fig. 5. 


COLLECTION EFFICIENCY FACTORS FOR 
COLLECTORS HAVING UNIFORM FLOW 
IN TUBES BONDED TO A FLAT PLATE 


This case, in one variation or another, occurs fre- 
quently in liquid-heating solar heat collectors. The 


general arrangement is shown in Fig. 3. A series of 0 
parallel tubes, of outside diameter d and inside diam- 50 100 
eter d, , are spaced w ft on centers along the collector UL Sqft-br * pound?F * btu 


plate by a continuous bond. For the purpose of these Fic. 5—The over-all efficiency factor F pg. 
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(such as Fig. 3) the temperatures ¢,. , ¢; , and ¢, will be 
the same at every tube. 

We make the following simplifying assumptions: 

(1) The thermal resistance around the cireum- 
ference of a tube is negligible in comparison with all 
other thermal resistances involved. 

(2) The division of incident sunshine between 
tubes and plate is such that the fraction d/w is ab- 
sorbed on the tubes and the remaining fraction | — d/w 
is absorbed by the plate. 

(3) The loss rate factors of the tubes and plate are 
also proportional to the ratios of their projected 
areas d'wand 1 — 

The following five independent equations can then 
be written, in each case expressing a heat transfer rate 
on a “per square foot of collector’’ basis: 

(1) The useful heat gain from sunshine striking the 
tubes: 


d 
Qu-tubes = 
l 


— R — t,)] [49] 


(2) The useful heat gain from the plate: 


(: “) Fla. — R— [30] 
w 


A consideration of the above equation, and of 
Fig. 3, will show that the above F is the fin efficiency 
factor of a fin of length (w — d)/2. 


3) The heat transfer through the bond: 
Qu-plate = — (t. — 


+) The heat transfer from tube to fluid: 


ri, h, 


w t, — ty) 


(5) The total useful heat gain: 
Gu = Qu-tubes 1 Qu-plate 
The preceding set of five equations contains five 
“unknowns” 
Gu Qu-tubes » Gu-plate » fe , bs [54] 


We should like to eliminate all but the first of these 
and, furthermore, to determine if the solution for q, is 
of the form: 

qu = — €R — — tad] [55] 
in which F’ is a function of the following variables 
only: 

: doh U,, and U, [56] 

w w 

The simultaneous solution of Equations [49] through 
(53) involves fifth-order determinants. The solution is 
straightforward, but the necessary calculations ap- 
peared too voluminous even for inclusion in these 
notes. However, if the solution is carried through, it 
will be found that the equation for g, is indeed of the 


form given by Equation [55], and that the value for 


C2U 


Furthermore, it will be found that Equation [57] 
can be manipulated into the following continuing 
fraction: 


CU, 


or substituting values of C; , C2, and C; : 


Ci + 


F’ 
wl, 
rd, h. 


p 
w 


(general case) 
The above equation is the general solution for this 
type of collector. Several particular cases are of in- 
terest. 


Case I. d w = 1.00 
In this case Equation [60] will be found to reduce to: 


l 
F’ = [61] 


wl, 
d, th, 
(d/w = 1.00) 


Note that the above case amounts to assuming 
that the tubes are side by side, with no space between 
them. This is essentially the case of the collector with 
uniform flow under a plate (see Fig. 2) except that 
the tube area in contact with the fluid is x d,/w sq ft per 
sq ft of collector. If rd,/w = 1.00, Equation [61] re- 
duces to Equation [38], as it should. 


Case Il. F = 1.00 


For this case Equation [60] reduces to: 


; 
F’ is: 
(CU. + Ur) + + Ue + 
3 1 
in which 
he (1 “) PC, = 
w w w 
= — [59] 
l 
[60] 
| 
VOL. 
(51) 
| 
| 
1 
w d U. 
w 
(F = 1.00) 
60 


The above equation may be used to estimate the 
efficiency factor of collectors having a relatively low 
thermal conductivity of the ‘‘bond’’ between plate and 
tubes. An example of such a collector is one formed of 
a grid of parallel pipes set inside a glass-covered box, 
with no continuous physical connection between the 
tubes and the inside of the bottom of the box. In 
such a case a reasonable estimate of the over-all con- 
vective plus radiative coefficient of heat transfer be- 
tween box bottom and tubes can usually be made, and 
appropriately substituted for U. in the equation. In 
making such a substitution, it should be noted that 
U’.. as used in these equations is the over-all heat trans- 
fer coefficient between plate and tubes on a “per 
square foot of projected area of tube’ basis. 


Case III. Thermal Resistance of Bond Negligible 


In this case 1/U’. = 0, and Equation [60] becomes: 
= i [63] 
a d 
(1 “) F 
w w 
= 0) 


As has been mentioned previously, the factor F in 
the above equation is the fin-efficiency factor of a fin 
of length (w — d)/2. Hence, from Equation [27]: 


, _ tanh bL 
F= bL [64] 


in which 


b = \ L=" [65] 

Equation [63] represents a case of considerable 
practical importance, as it is applicable both to col- 
lector plates made of tubes soldered to a metal sheet 
and to plates made of ‘‘tube-in-strip”’ or similar integral- 
tube type materials. Most water-heating collector 
plates use one or the other of these constructions. 

We have found that one convenient method of using 
Equation [63] in design calculations is to plot F’ 
against tube spacing for a variety of pertinent condi- 
tions. Examples of such graphs are given in Fig. 6 
through 11. 

Figs. 6, 7, and 8 are plotted for a constant UU, value 
of 1.2 Btu per sq ft-hr-°F. This value is typical of col- 
lector plates insulated on the sunny side with one air- 
spaced layer of glass or transparent plastic, and well 
insulated on the back. Such construction is often used 
in domestic water-heating collectors. The three graphs 
present F’ for copper, aluminum, and steel sheets 
respectively, and in each case for several sheet thick- 
nesses and for two values of the tube-to-liquid film 
coefficient h.. The lower value of h, (60 Btu/sq ft- 
hr-°F) is more or less typical of laminar water flow in 


thermosyphon-type collectors. The higher value of 
h, (300 Btu per sq ft-hr) is in the range of turbulent 
water flow, and may be obtained with pump circula- 
tion. 

Figs. 9, 10, and 11 present the same data for a 
U, value of 2.5 Btu per sq ft-hr-°F which is a value 
typical of a collector with no cover on the sunny side, 
moderately insulated on the back, and operating under 
low wind velocity conditions (3 to 5 mph). This type 
of construction may be used for solar heat collection at 
a temperature level not too much above that of the 
outdoor air, for example in swimming pool heating 
and in some space-heating applications. 

Figs. 6 through 11 are all calculated for 4-in. diameter 
tubes, and it is assumed that the tube walls are fairly 
thin, i.e., d, © d. Under most conditions no large error 
will result from using the graphs for any tube diameter 
between 3 in. and ? in. 

It will be remembered that the various expressions 
for F’ derived in this section are all based on the as- 
sumptions that the flow in the parallel tubes is uniform 
from tube to tube and in the same direction. These 
assumptions are not realized in collectors in which 
the flow is unevenly distributed between tubes, nor 
are they realized when a return bend circuit arrange- 
ment gives flow in opposite directions in adjacent 
tubes. Analytical solutions for these cases are de- 
sirable but do not seem to be available. In general, 
qualitative consideration will show that the F’ ratios 
for collectors with either uneven or opposed direction 
flow will be somewhat less than the F’ 
in this section. 

Turning to a consideration of the factor Ff’, for col- 
lectors of the general type discussed in this section, we 
note that Equations [39] and [41] of the preceding 
section are applicable. Hence the expression for Fx 
(Equation [48]) developed in the preceding section is 
also applicable. The equation is plotted in Fig. 5. 


ratios covered 


ILLUSTRATIVE EXAMPLES 


A few numerical “problems” and their solutions 
may be helpful in illustrating the practical utility of 
the efficiency-factor equations. 


Example | 

A water-heating collector designed for swimming 
pool heating uses an integral-tube type of collector 
plate having 4-in. tubes spaced 4 in. on a 0.025-in. 
copper plate. The collector plate is uncovered, wind 
velocity at the level of the collector is in the range of 
3 to 5 mph, and the loss rate factor U’, is estimated as 
2.5 Btu per sq ft-hr-°F. 

(a) If the flow rate through the collector is 3 |b 
per sq ft-hr, and the flow in the tubes is laminar, what 
is the efficiency factor Fz ? 
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Solution: From Fig. 9, we find F’ = 0.88. Enter- 
ing Fig. 5 with this value of F’ and with 


CG/U, = (3 X 1.0)/2.5 = 1.2, 


we find Fy = 0.62. 

(b) Assume that the absorptivity of the painted 
collector plate for solar radiation is 0.95, and that its 
emissivity and absorptivity in the long-wavelength 
region is 0.90. Assume further that the collector is 
operating under the following conditions: 


Entering water temperature............... 82°F 
Outdoor air temperature.................. 65°F 
Intensity of sunshine incident upon collector 
275 Btu ‘sq ft-hr 
Net outgoing radiation R’ from a black surface 
at air temperature and at collector tilt 
20 Btu sq ft-hr 


What will be the useful heat collection rate of the 
collector in Btu per sq ft-hr? 
Solution: 
Qu = Frilag; — — — ta)| 
= 0.62/0.95 K 275 — 0.90 K 20 — 2.5(82 — 65)| 
125 
(c) If the flow velocity in this collector is such 
that increasing the flow rate to 7 lb per sq ft-hr would 
result in turbulent flow, how much will the useful heat 
collection be increased if the flow rate is increased to 
this value? 
Solution: From Fig. 9 we find F’ = 0.95 for turbu- 
lent flow. Entering Fig. 5 with this new value of 
F’, and with the new value of GC/U, = 2.8, we 
find the new value of F, to be 0.80. Under the 


previous laminar flow conditions we found Fy was 
0.62. Hence the change to increased flow rate and 
turbulent flow will increase useful heat collection 
by the ratio 0.80/0.62, or about 29 per cent. 


Example 2 

A manufacturer of glass-covered domestic water- 
heating collectors produces collector plates by solder- 
ing parallel 3-in. diameter tubes to a metal sheet. 
Several types of sheet metal are available at varying 
costs, and the tubes can be spaced as desired. The 
problem is to determine the particular choice of sheet 
material and tube spacing which will yield the maxi- 
mum ratio of heat collection per unit of initial cost. 

The manufacturer has found that the cost of a com- 
plete collector can be broken into three components, as 
follows: 

(a) A fixed cost of $3.00 per sq ft of collector, which 
covers framing, insulation, glass, supports, painting, 
etc. 

(b) The cost of the sheet metal for the collector 


plate. There are four available sheets: copper sheet in 
three thicknesses, 0.010-in., 0.020-in., and 0.030-in. 
and costing $0.40, $0.80, and $1.20 per sq ft respec- 
tively. Steel sheets are available in 0.040-in. thickness 
costing $0.25 per sq ft. 

(c) The total cost of tubes, soldering, and asso- 
ciated headers. This has been found to be proportional! 
to the length of tube used, and amounts to $0.35 per 
lineal ft of tube. 

The flow in the tubes is laminar. The flow rate, for 
the purpose of simplifying this example, is assumed to 
be the same regardless of tube spacing. 

Solution: The cost of the energy furnished by the 
collector will be directly proportional to the total 
cost per square foot of collector, and inversely 
proportional to the collector efficiency. Hence, if 
we let the collector cost be S dollars per sq ft, 
the problem amounts to determining the lowest 
possible value of the ratio S/F’ for the various 
choices of sheet material and tube spacing. 

From the data given, it is readily seen that the total 
cost S can be expressed: 


S = 3.00 + B + 4.20/w’ dollars/sq ft 


in which B is the cost of the sheet metal used, and w’ 
is the tube spacing in inches. 

The elegant method for determining the minimum 
value of the ratio S/F’ is to find the point at which 
the partial derivatives of S/F’ with respect to sheet 
metal cost and with respect to tube-spacing cost are 
simultaneously equal to zero. However, this problem, 
and many similar problems of this type, can be solved 


Tasie [—Trial and Error Solution of Example 2 


Tube spacing Collector cost’ Efficiency S/F’ 


Type of sheet w’ in. S ($/sq ft factor F’ $/sq ft—F’ 


.00 
10 
13 
.16 
.20 


.010” Copper 
.010” Copper 
.010” Copper 
.010” Copper 
.010” Copper 
.010” Copper 
.010” Copper 
.010” Copper 


™ 


Or or or or 


.020” Copper 
.020” Copper 
.020” Copper 
.020” Copper 
.020” Copper 


.030” Copper 
.030” Copper 
.030” Copper 
.030” Copper 


.040” Steel 
.040” Steel 
.040” Steel 
.040” Steel 
040” Steel 


an 
4 
4 
Le 
- 
0 98 
0 
0 77) 
0 77) 
ae. 0 78 j 
0 7.0 4.40 860 12 
0) 6.0 1.50 (9.07 
0 5.0 1.64 910 
0 4.0 1.85 .933 .20 
Pea. 0 3.0 5.20 954 46 | 
0 7.0 4.80 878 
0 6.0 4.90 901 (9.44 
0 5.0 5.04 .920 47 
0 4.0 5.25 938 “60 
0 7.0 3.85 .725 5.31 
0 6.0 3.95 .780 
0 5.0 835 1.90 
0 4.0 1.30 (4.87) 
0 3.0 4.65 5.00 


equally quickly by a trial-and-error tabulation. For 
this problem we use the values of S given by the above 
equation together with corresponding values of F’ de- 
termined from Figs. 6 and 8. The resulting tabulation 
is shown in Table I. The minimum values of the ratio 
SF’ for the various sheets are those enclosed in 
parentheses. For the particular set of cost and other 
conditions specified in this example, the collector 
plate which ‘delivers the most Btus for the least 
money”’ is seen to be that using 0.010-in. copper sheet 
with the tubes spaced approximately 53 in. 


CONCLUSION 


The efficiency factors derived in these notes are 
applicable to a wide variety of solar heat collectors. 
They are also directly applicable to other types of 
panel heat-exchangers, for example to floor or ceiling 
panels used for interior temperature control. 

The intelligent use of the efficiency factor concept 


can easily eliminate a large area of guesswork in the 
design of a solar heat collector. Although, as has been 
previously mentioned, these efficiency factors are not a 
new concept in solar heat collector design, they do not 
seem to have been used as much as they should be. It 
is hoped that these notes may be helpful in furthering 
the use of these basic heat-transfer concepts and in 
reducing unnecessary empiricism in the design of heat 
exchangers used for collecting solar energy. 
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A Stationary Solar Energy Concentrator 


By Paul B. Archibald 


A parabolic reflector having a short focal length 
will keep the sun’s rays focused within a small 
area for several hours without the employment of 
a heliostat mechanism. Its use as a steam genera- 
tor and as a stove is described. 


Among the more frequently proposed applications 
for solar energy is its employment in cooking. While 
many clever solar stoves have been built, their wide- 
spread use has been restricted mostly to countries poor 
in natural fuels. Americans do not live in such a coun- 
try, and it is difficult to see how solar stoves could ever 
be of much practical importance here. This is not to 
say that a solar stove has not a certain appeal. The 
smug feeling of complacency which the solar broiler 
generates in the backyard chef is to be considered, and 
one is led to surmise that this feeling might be further 
enhanced if the cook had made the stove himself. 

Some experiments directed toward the development 
of a solar steam generator indicated a design for a solar 
stove which would be easy to construct out of material 
normally available in the home. 

The steam generator, from which the solar stove was 
derived, had as its objective a means of concentrating 
sunlight which would not involve the use of any moving 
parts. The parabolic reflector that served as a concen- 
trating device was made in the following manner. A 
parabola having a focal length of 1 in. and legs of about 
16 in. was drawn on a piece of 3-in. plywood, and two 
parabolas of this size were cut out with a band saw. 
The area outside the parabolic curve was used to shape 
and support the reflector. These two pieces were nailed 
to a board 36 in. apart in such a manner that the open 
ends of the parabolas both pointed in the same direc- 
tion. Into these parabolas a piece of 6-mil hardened 
aluminum sheet was draped and then attached with 
tape. 

The reason for the short focal length of the reflector 
was to eliminate the need for a heliostat mechanism. 
The reflector was set level on the ground with its axis 
in an east-west direction and then tilted into the plane 
of the sun. This last adjustment was made by observing 
the shadow of a stick held midway between the two 
surfaces of the parabola and tilting the reflector until 
the shadow fell at the bottom or vertex of the parabola. 


Livermore, California 


Once the adjustment was made the reflector was not 
moved. The focal point of the sun’s rays coming in 
obliquely in the morning and afternoon would be inside 
the l-in. maximum which was obtained when the sun 
was at its zenith. Readjustment to the plane of the sun 
was necessary only every three or four days. 

The boiler tube consisted of a blackened copper tube 
-in. outside diameter within two concentric glass tubes, 
the outside one being 14 in. diameter. The boiler tube 
was placed at the focus of the parabola in such a man- 
ner that the upper side of the copper tube was at the 
1-in. focal point. 

The focal point of the sun’s rays, which in the morn- 
ing and afternoon would strike the reflector obliquely, 
is equal to f = f; cos? a, where f; is equal to 1 in. and 
a is the angle of the sun from the zenith. With such an 
arrangement the 3-in. copper tube is at the focal point 
for a period of about five hours, and the temperature 
‘an be estimated by multiplying the solar constant by 
the ratio of the projected area of the reflector divided 
by the area of the copper tube. With the reflector and 
the boilertube described this ratio was 12, which gives 
a value of about 24 cal per sq em per min. The equi- 
librium temperature for this energy input: is, by the 
Stefan-Boltzmann Law, 463°C. All of the energy is not 
available at this temperature; however, the quantity 
and temperature of the steam which can be obtained 
‘an be readily calculated. 

In the model generator water boiled for seven hours 
and, when empty, the tube reached a temperature of 
350°F during the middle of the day. There were several 
apparent reasons for this low temperature: the boiler 
tube was too long—it should be shorter than the re- 
flector, so as to receive the full 24 cal per sq em per min 
along its entire length when the sun is not at its zenith; 
the glass tube outside the copper tube was not evacu- 
ated, and water vapor would condense on the inside of 
the glass. While none of the mechanical problems would 
appear to be difficult to solve, they would require some 
effort. 

Despite the shortcomings of the first model, its per- 
formance was gratifying, and the idea of a non-moving 
reflector for concentrating solar energy was demon- 
strated. The improvements that would have to be made 
plus the cost of instrumentation before worthwhile 
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data could be collected indicated that as a backyard 
project further development was not feasible. 

However, the use of the reflector as a solar stove 
seemed to be a possibility which could be investigated. 
The shape of the food to be cooked should be long and 
round like the boiler tube. Ears of corn met the require- 
ments, and it was found that corn with the husks re- 
moved and wrapped tightly in cellophane and laid at 
the vertex could be cooked in about forty minutes. The 
flavor seemed much better than boiled corn, but no 
effort was spent on trying to get an impartial opinion. 

One interesting trick was to set the stove with the 
corn in it in the shade so that the sunlight would not 
strike it until about forty minutes before meal time. 

While the reflector described was quite adequate as 
a corn-cooking device, it seemed to be a little too elegant 
for such a prosaic job, and it appeared that a stove of 
more common material and simpler construction would 
be desirable. 

For a reflector surface ordinary aluminum foil glued 
to cardboard was used, and for the parabolic shape stiff 
wire (from a coathanger) was bent into the desired 
shape and held in position by a string across the 
branches of the parabola. The reflector surface was 
bent around the wire parabola and taped in place. The 
stove was then placed ina pasteboard box for support. 


While the amount of solar energy collected is propor- 
tional to the area of the reflector normal to rays of the 
sun, a shallow—say 8-in. focal length—reflector inter- 
cepting the same amount of sunlight did not function 
as a stove. The reason for this seems to be the heat loss 
by convecting air currents. 

The paper and aluminum stove is a very safe device 
even children can operate it—and if of no practical 
significance it makes a rather nice demonstration of 
solar energy. 

The steam generator can be optimized in various 
ways. Thus, to increase the diameter of the copper tube 
from 2? in. to 4 in. will increase the length of time that 
the tube is at the focus from five to seven hours. How- 
ever the equilibrium temperature would be lowered. It 
is not necessary for the boiler tube to be on focus in 
order for it to intercept a considerable amount of solar 
energy. This explains why, in the model, boiling was 
obtained one hour before the focal point moved up to 
the copper tube. To use a shorter-focal-length reflector 
would mean that it would have to be much deeper in 
order to intercept the same amount of sunlight. 

While the stationary solar energy concentrator de- 


scribed is probably not as efficient as one with a helio- 


stat mechanism, its simplicity of construction should 
more than outweigh this disadvantage. 
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Calculation of the Concentration of knergy at Points 
Outside the Focal Spot of a Parabolic Condenser’ 


By Alfred W. Simon 


Scientific and Engineering Staff, Air Force Missile Development Center, Holloman Air Force Base, New Mexico 


An analysis of the distribution of intensity 
and of energy falling outside of the circular pro- 
jection of the cone of reflected radiation coming 
from all infinitesimal areas of the surface of a 
paraboloidal mirror. The results given in a pre- 
vious paper under the same title are incorrect. 
The correct theoretical deductions and the 
correct numerical results are given here. A 
uniform sun is assumed. The original paper 
which is corrected here appeared as Chapter 7 
of Solar Furnace Support Studies, AFMDC-TR- 
58-7, Holloman Air Force Base, New Mexico, 
and on pages 22-24 of Volume II, number 2, 
Solar Energy, April 1958. 


Under the above title in a previous paper I have 
calculated the variation in intensity of the energy 
falling on the focal plane of a parabolic condenser as a 
function of the distance from the center of the spot. 
However, I subsequently discovered that the theoretical 
deductions made in that paper are incorrect and the 
published results, therefore, wrong. In the present paper 
I give the correct theoretical deductions and the cor- 
rect numerical results. 

Consider (Fig. 1) the cone of light of angular width 
a, Where @ is the angular diameter of the sum, re- 
flected at an angle @ to the axis of the parabola. The 
base of this cone of light through the focal point F of 
the parabola makes an angle 6 with the focal plane and 
has a radius R given by the equation: 

R = pa/2 {1] 
where p is the distance PF. This distance is given in 
turn by the equation: 

p = 2f/(1 + cos 8) [2] 
where f is the focal length of the parabola, so that we 
have R given by the equation: 

R = fa/(1 + cos 6) [3] 

The circle constituting the base of this cone of light 
is projected onto the focal plane, thus forming an 


* Originally issued as AFMDC-TN-59-11, ASTIA Document 
No. AD 215462. Released by author May 7, 1959. 


Fic. 1—Geometry of the reflection of a cone of light onto the 
focal plane of a paraboloid. 


ellipse of semiminor axis b given by b = R and semi- 
major axis a given by a = b/cos 6. Hence we have: 
b = fa/(1 + cos 6) [4] 

a = fa/(1 + cos @) cos 6 [5] 


For a focal spot of 5 in. in diameter, that is, for 


fa = 5 in., we have the values of a and b corresponding 


to various values of 6 given in Table I. 

The corresponding ellipses are plotted in Fig. 2 
through 7. 

The amount of energy falling into each ellipse is 
that falling on a ring of radius p sin @ and thickness 
p dé on the face of the paraboloid. In particular, if we 
denote this energy by dE, we have dF given by: 

dE = 2mep sin 6 dé [6] 


where ¢ is the solar constant. 
It will be seen that all the ellipses enclose the focal 
spot of radius r; = 2.5 in.; hence the fractional portion 
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Ellipsoidal spots produced in the foeal plane by cones of light reflected at various angles from the paraboloid. 


of the energy collected at an angle 6 falling into the 
focal spot is given by the ratio: 


xr; mab (7] 
and the corresponding intensity d/; is: 


= rr, sin 6 dé 


rab rr,” 
= sin 6 cos 6 


Consequently, /; is given by: 
Se fe . 
sin 6 cos 6 dé = — sin’, 
vi 
where 4, is the half angle of opening of the paraboloid. 
If now we select a ring of inner radius r; and outer 
radius r,, where r, > r; > r;, this ring will intercept 
a certain area A» of the ellipse corresponding to the 
angle @. 
The energy dF falling into this area will be given by: 


dE = 2rep sind dé 110) 


Trav 


As the point P moves around the face of the para- 
boloid the corresponding ellipse rotates, and this energy 
is scattered over a ring of inner radius 7; and outer 
radius r,. Consequently the corresponding average 
intensity d/, within the ring will be: 


Ag sin 6 dO 


dI, = [11] 


rab x(r,2 — 
and the ratio of this intensity to that in the focal 
spot will be given by: 


dl, Ag sin 6 dé 
I; rab (r,2 — 4e sin? 6; 
Ag 2f'a’ sin 6 d6 
gab (1 + cos 6)2(r,2 — r2) sin’; 


This equation becomes: 


dl Ae 8 5° sin 6 
— = 9 9 9 [13] 
I; rab 3(1 + cos — 


if we substitute for fa its value of 5 in. and for sin” @ 
its value ? corresponding to a diameter of the focal 


spot of 5 in. and a value of 6; equal to 60°, respectively. 
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TABLE I[] 


TABLE III 


T cal sec 


1.000 1017 
0.536 545 
147 149 
3: 46 

0 


Since Ag» is a complicated function of 6, we have to 
resort to approximate integration, rewriting the above 
equation in the form: 


Al, Ae 2000 sin 6 


[14] 


I; rab 3 X 57.3 (1 + cos 0)? (r2 — r?) 
where we have taken A@ equal to 10°. 

The various rings of radii 2.5, 3.5, and 4.5 in. have 
been drawn in Figs. 2 through 7. 

In order to calculate the areas Ag we require the 
intersections of the various circles with the ellipses. 

For these intersections we have, obviously: 


= taV(r2 — / (a — [15] 


The values of the abscissas of the points of intersec- 
tion of the ellipses with the circles of radii r; = 2.5 in., 
ro = 3.5 in., and rs; = 4.5 in., are given in Table II. 

As a typical calculation of the areas A¢ let us take 
the area intercepted between the radii 7; and ry for 
6 = 50° (shaded area in Fig. 6). 

We have Aso given by: 


re 
As = 4 Va — +4 Vr? — 
Jz, a 


hs [16] 
— 4 | Vr? 


5 


Fic. 8—Variation of average intensities within rings of various 
diameters in the focal spot plotted against the diameter 
of the center of the ring. 


where b/a = cos 6. Substituting for x, x2, 71, 72, and 
§ their values and integrating, we find Aso = 9.06 sq in. 
If now we note that for @ = 50°, rab = 45.25 sq in., 
sin 0/(1 + cos @)° = 0.284, and that for this case 
Equation [16] reduces to 

Al, Ag sin 6 


= 1.454 
I; + cos @)? 


[17] 
we obtain A/,/7; = 0.08267. 

If now we calculate the values of A/,//; correspond- 
ing to various values of @ for a given set of radii and 
apply the trapezoidal rule to the results for each ring, 
we have the average intensities within the various rings 
as given in Table III. 

The intensities have been calculated on the basis 
of 1.76 cal per sq cm per minute received by the pa- 
raboloid. In Fig. 8 the calculated values have been 
plotted against the radius of the center of the ring. It 
should be noted that the intensities are average values 
over the area of the ring, that the results hold only 
for a uniform sun, and that no losses have been as- 
sumed. 
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Comments on Papers 


Dear Sir: 

The paper, “Theory and Design of a Mechanical 
Blackbody for Solar Radiation”? by Winskell and Hsu, 
which appeared in the April, 1959, issue of Solar Energy, 
contains a number of errors or inconsistencies. Some 
are important technically; others probably are printing 
errors. 

The first sentence of the Introduction states that 
“one of the current problems in the field of solar energy 
utilization is to obtain a durable, high absorptive, low 
emissive surface to solar radiation.” This is technically 
incorrect, or at least misleading. Besides being durable 
and of high absorptivity the actual requirement is for 
a surface of low emissivity not to solar radiation but 
to the radiation from the collector surface itself which, 
of course, is at a temperature much lower than that of 
the sun. 

The reader is likely to be confused by the expres- 
sions: “angle of incidence,” “angle of radiation,” and 
“angle of surface to radiation,” which apparently are 
used interchangeably. The ‘tangle of incidence,” of 
course, refers to the angle between the incident rays 
and the normal to the surface; but it is not obvious 
which surface is meant. The other two expressions are 
ambiguous until it is determined that the former refers 
to the surface of the sides of the parallelopiped element, 
while the latter refers to the bottom surface of the 
parallelopiped element. 

With respect to the equations above Fig. 12, 


A-(B+C) and H=E — (F+@G) 


the author states that “if C is assumed equal to G, 
D and H are proportional (equal) to (A — B) and 
(E — F), respectively.” This statement is not true ex- 
cept for the trivial case when C and G are both zero. 
Otherwise all that can correctly be said is that, if C and 
G were constants and equal to each other, D would vary 
directly with (A — B) and H would vary directly with 
(E — F). 

In the development of the theory in Phase I the 
assumption was made that all radiation arrived at the 
same angle of incidence. Yet in the description of the 


experimental work in Phase II, the fact that the radia- 
tion arrives at widely varying angles of incidence is not 
mentioned. 

Among the smaller errors noted were: 

(a) In the equation 


Absorptivity = 1 — (1 — 


the exponent n is omitted. Or, perhaps, the error was 
the misplacing of the quantity, (2b tan @/a), which is 
equal to n, and which should have appeared as the ex- 
ponent of (1 — a). 

(b) In the equation immediately following, the ex- 
pression 2b tan @/a in the right hand member is incor- 
rectly shown as a multiplying factor when it should 
have been the power to which the term (1 — a@) is 
raised. 

(c) The authors mention the results of ‘the investi- 
gations of Daniels and Salem” but give no reference as 
to where this work may be found. Also they state that 
“the investigations of Daniels and Salem give values of 
absorptivity to solar radiation and emissivity at some 
of the expected operating temperatures of the surface 
of brass, copper, and copper oxide coatings. This in- 
formation is given in Fig. 14.”’ The only information 
given in Figure 14 is for cupric oxide coatings and none 
of the expected operating temperatures is found. 

There are several annoyances which hardly mislead 
but do detract: in Fig. 12 and 13 it is hard to tell which 
curve is which since it is difficult to distinguish the 
O’s and X’s; the dual ordinate scales, presenting both 
potentiometer readings and temperatures, are quite 
unnecessary; also, for no apparent reason, these dual 
scales are interchanged from one figure to the next; in 
Fig. 12 the intended reference is to the Handbook of 
Chemistry and degrees Fahrenheit has nothing to do 
with it; in Phase I the length-width ratio is referred 
to as b a while in Phase II the symbols B/ A refer to 
the same quantity; the lower abscissa scale of Fig. 7 is 
unnecessary since it repeats information contained in 
Table I on the same page (or vice versa). 
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AKIC-EEI Heat Pump Committee. “First step of 
heat storage research program completed: progress 
report.””. Bull. Edison Elec. Inst. 27(6): 250-51, 
June 1959. 


Summary of reports by Battelle Memorial Institute and by 
Comstock and Wescott, Inc., on first phase of a project to de- 
termine the status of heat storage development and to analyse 
the potential applicability of various heat storage materials 
for use in conjunction with the heat pump and resistance heat 
ing. Battelle’s report indicates that high-temperature storage, 
charged with resistance heating properly controlled, offers the 
most immediate promise as a means of smoothing load curves 
with either heat pump or other electric heating systems. 
Second choice is condenser-side storage charged by the heat 
pump; third choice, evaporator-side storage charged with 
‘free’? heat. Comstock and Wescott’s report recommends that 
future research and development work be directed to study of 
three systems: (1) heat of hydration systems; (2) heat of fusion 
systems; (3) heat of chemical change systems. 
* 


* * 


Battery Research and Development Conference, 12th, 
It. Monmouth, N.J., May 21-22, 1958. Proceedings. 
Power Sources Division, U.S. Army Signal [ngi- 
neering Laboratories, Ft. Monmouth, N.J. 134 p. 
Illus. 


Includes papers on fuel cells, nickel-cadmium batteries, 
special purpose batteries (automatically activated zine-silver 
oxide batteries), new primary battery designs, low temperature 
batteries, and other power sources of the future. Papers of 
particular interest in solar energy utilization are abstracted 
below. 


Hunrath, George. “Solar energy conversion.’ (In: 

Battery Research and Development Conference, 
12th, Ft. Monmouth, N.J., 1958 Proceedings. p. 
111-13. Illus.) 


Aretie regions can benefit from solar energy utilization for 
power requirements of small magnitude because in summer the 
high latitudes receive 24 hours of sunlight per day, even though 
it is at low intensity. Three methods for the conversion of 
solar into electrical energy offer good possibilities: silicon solar 
cells, thermionic emission (Edison effect), and thermocouples. 
Concentrators are required for the last two methods, and 
possibly the first, in northern regions. Conversion efficiencies 
of 10 per cent have been achieved with the first; 20 per cent 
with thermionic emission is possible; present conversion effi- 
ciencies with thermoelectric generators are fairly low. Silicon 
solar cells or a cathode and collector in a thermionic emission 
diode are capable of furnishing power in the Arctic for such 
applications as navigational aides, and unattended com- 
munication or meteorological stations. Energy can be stored 
in storage batteries for peak power demands or to furnish 
power during the winter months. 


* 
Kittl, Emil. “Thermal energy conversion.” /bid. p. 
113-15. Illus. 


The thermal energy conversion aspects of solar and nuclear 
energy are discussed. Additional concentrating equipment is 


Solar Abstracts 


needed for thermoelectric converters utilizing solar energy, and 
only if considerably higher efficiencies are achieved will this 
system be comparable to the solid state solar battery. One of 
the best approaches to meeting military requirements for a 
portable power source is in the field of solid state thermal 
energy conversion. The thermionic diode utilizing the Edison 
effect is discussed in some detail, and a thermionic diode 
powered by concentrated sun rays and another using nuclear 
radiation are described. The load characteristics and output 
voltage vs. output current are compared for the two diodes. 
A large difference of the emitter and collector work functions 
is required for large output voltage of the device, and high 
operating temperature is also necessary to obtain maximum 
output power. Computed values for the theoretical maximum 
power output and maximum efficiency of a thermionic diode 
converter are compared using a group of promising emitter 
materials. Since the maximum efficiency of the thermionic 
converter is about three times the theoretical efficiency of 
the silicon cell, it appears that the thermionic diode offers 
promise of providing a system having a higher conversion effi- 
ciency and higher power output per unit size and weight than 
any other converter. 


Ziegler, Hans K. ‘Energy sources for ‘space-age’ 
activities.”” Jbid. p. 1000-4. Illus. 


Various types of energy sources for artificial satellites and 
space vehicles are compared, taking into account the basic 
requirements and unique operating conditions imposed by the 
outer space environment. Chemical batteries represent the 
best solution for short time and lower power requirements. 
They offer up to 80 watt-hr per watt-hr/lb. For long periods of 
operation in the low, medium, and high-power range, solar 
conversion systems—with 2600 to 3100 watt-hr/lb for non- 
storage type and 1200 to 1600 watt-hr/Ib for storage type—if 
one-year life of solar cells is confirmed, or better, if life is 
longer—appear most favorable. These values apply to space 
vehicles with no orientation control. Nuclear power devices 
with presently 500-2000 watt-hr/lb are expected to improve in 
efficiency and eventually to provide highly reliable sources 
for medium and high-power and long-life requirements 
although competition of solar sources, depending on their 
long-range behavior and practical reliability, may be quite 
close. In this area little field experience is vet available. 


Bliss, Raymond W., Jr. “A general description of the 
solar energy laboratory and its initial research 
program.’ Arizona, University, Institute of At- 
mospheric Physics, May 1, 1959. 45 p. Illus. 


The Solar Energy Laboratory at the Institute is designed 
and built for work in two fields: (1) applied research on the 
development of moderate-temperature thermal processes for 
the direct utilization of sunshine; (2) meteorological research 
on problems involving the measurement of solar, atmospheric, 
and ground radiation. This report provides a general descrip- 
tion of the Laboratory and a preliminary working guide for the 
laboratory personnel. Construction features of the building 
are given and components of the solar-heating-cooling system 
described in detail. A plentiful supply of versatile instrumenta- 
tion has been installed for the experimental system and for 
other projects. For its first two years of operation, a major 
portion of the work will center around the study and analysis 
of the performance of the heating-cooling system. Other 
projects will include maintenance of the local pyrheliometer 
records, atmospheric transmissivity studies, etc. Potential 
graduate student projects are listed. 
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Bliss, Raymond W., Jr. ‘Solar water heating system.” 
U.S. Patent 2,838,043, June 10, 1958. 5 p. Illus. 


In this simple, relatively inexpensive solar water-heating 
system the heat collector is placed in its normal elevated 
position and the storage tank at a low level, and water is cireu- 
lated during sunshine hours by a small pump. Feed water is 
supplied by connecting a source of cold water under pressure 
to the bottom of the tank. An air space is provided, optionally 
at the top of the tank, above the water outlet, which has a 
greater volume than the portion of the system above the air 
space which is to be drained. By providing a bleeder line 
between the air space and one of the lines providing circulation 
between the hot water storage tank and the collector, water 
will automatically flow from the collector to the tank when 
pump action is discontinued, and air will flow through the 
bleeder line from the air space into the collector and into the 
portion of the svstem above the air space 


* 


te- 


‘alifornia, University, Institute of Engineering I 
search, Berkeley. “Sea water conversion program: 
a technical and economic investigation of certain 
schemes for producing potable fresh water from sea 
water. Berkeley progress report for the vear ending 
June 30, 1959.” Series No. 75, Issue No. 15, July 
24, 1959. 28 p. Illus. 


Reports of investigations at the University of California 
at Berkeley and the Richmond Field Station into various dis 
tillation processes, including multiple effect centrifugal evapo- 
rators, vacuum flash distillation, solar distillation, reversed 
osmosis using capillary ultra-filters, ion exchange processes, 
freezing processes, vapor compression tests, magnetic effects 
investigation, biological conversion of saline waters, and 
nuclear reactor studies. The publication includes reports by 
H. W. MeCracken, D. M. Campbell, K. D. Frohman, G. M. 
Voge and J. W. Weldy on recent work at the University in the 
construction of various types of solar stills. A small unit made 
of wood and costing about $60 was the subject of an instruction 
sheet for distribution in response to inquiries received on small, 
low-cost stills; its vield varies from 2? to } gal per day. As the 
result of detailed discussion by a panel of authorities in the 
solar distillation field, some generalizations and recommenda 
tions for specific projects in the Sea Water Conversion Program 
are formulated 


Diamond, P. M. ‘The central problem of large scale 


power generation in space waste heat disposal.” 


Institute of Aeronautical Sciences Paper No. 59-96; 
presented at the IAS meeting, Los Angeles, Calif., 
June 16-19, 1959. 36 cols. Illus. 


\ discussion of certain aspects of the characteristics of 
large power systems designed to operate in a space environ 
ment, with special emphasis on the problem of the weight of the 
unit. The writer concludes that for lower power, very long 
duration applications the solar battery appears to offer the 
greatest potential. The solar turbo-generator systems are 
among the lightest in weight for large scale power generation. 
The limitations of the radiant energy collector, however, may 
ultimately preclude the use of this device for large units 


Khanna, Mohan Lal; Gardner, A. L.; Davey, T. N.; 


Suri, S. P. “Plane-glass mirror solar energy con- 


centrators for concentrating sugarcane & palm 
juices.” J. Sei. Ind. Res. 18A(5): 212-16, 1959. 
Illus. 


The design, construction and working details of a plane- 
glass mirror (9 in. square) solar energy concentrator are 
described. The results of experiments on the evaporation of 
water and cane juice, using reflectors with over-all dimensions 
of 3 by 3 ft and 6 by 3 ft, are presented. The advantages of this 
type of collector over the conventional concentrators, such as 
low cost, simplicity of construction, ease of handling, and 
portability are pointed out, and the probable economy in the 
sugar cane and palm gur manufacturing industries by the use 
of this type of reflector, is assessed. (authors’ abstract) 


Kiepenheuer, Karl. The sun. Ann Arbor, Univ. of 
Michigan Pr., ¢1959. 160 p. Illus. (Trans. by A. J. 
Pomerans from Die Sonne. Berlin ete., Springer- 


Verlag, 1957.) 


Nontechnical account of what we know about the sun, in- 
cluding the nature of sunlight, sunspots, northern lights, solar 
storms, the surface and interior of the sun. Solar energy utili- 
zation is mentioned briefly. 


Telkes, Maria. “Apparatus for storing and _ releasing 
heat.” Patent 2,808,494, Oct. 1, 1957. p. 
Illus. 


Provides a system for storing heat at a relatively high tem- 
perature, in the range 300° to 550°F, with a capacity several 
times that of previous systems. The system makes use of the 
heat of transition of a dimorphic substance, such as sodium 
sulfate, alone or admixed with small quantities of other salts 
as a heat-storing and releasing device. 


Zwick, Eugene B.; Zimmerman, Robert L. “Space 
vehicle power systems.” J. Am. Rocket Soc. 29(8): 
553-64, Aug. 1959. Illus. 


A review of systems which can be utilized for the generation 
of secondary and/or propulsion power for advanced space 
vehicles. Power systems considered include: chemically-fueled 
power supplies, open and closed cycle; solar power supplies, 
closed evele; nuclear power supplies, closed evele; thermionic 
and thermoelectric power supplies; solar photovoltaic and 
solar reeyeling fuel cell power supplies. Also included in the 
review are the results of an initial program relating to the 
optimum temperature ratio, collector and radiator size for 
general power systems when used with nuclear, chemical, 
or solar energy sources. Descriptions are presented of a number 
of typical power systems. An attempt is made to describe the 
general physical construction as well as to indicate the limita 
tions of the system. Analyses and estimates have been prepared 
in order to determine the system weight as a function of power 
level and duration. Comparison of some of the systems studied 
is presented. Weights and sizes of the significant components 
are tabulated for 1- and 100-kw systems. Variation of system 
weight with time and power level is plotted to determine the 
optimum system from a weight standpoint for various applica- 
tions. System selection for various missions is summarized, and 
the role of turbomachinery in space vehicle power systems is 
noted. (authors’ summary ) 
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